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Introduction 



At the beginning of the past century, two basic blocks of the modern physics were established: 
Quantum Theory and the Theory of Relativity. In the following years, Quantum Field 
Theory was developed in order to made both theories compatible (or, to be more precise, 
Quantum Mechanics and the Special Theory of Relativity, as the complete connection with 
General Relativity is still an open problem). One of the surprising points of this theory 
is that some of the calculations involved divergent quantities. At first this was thought to 
be a problem, but soon it was found that it was inherent to any Quantum Field Theory 
calculation, reflecting only the infinite degrees of freedom of these processes. Finally, the 
standard way of treating these divergences was established to be a two-step work: 

• First of all, we have to regularize the divergence. The idea is to introduce an extra 
parameter (the regulator) in terms of which we can rewrite the divergent expression as 
a finite function, being the infinite result a certain limit of this parameter. 

• Once we have parametrized the divergence in terms of this regulator, the next step 
is to drop off the divergent part of each diagram, retaining only the finite part. This 
is what we call renormalization. The easiest way of performing this is to modify 
the coefficients of the terms of the action, making them regulator-dependent, so that 
we have new terms in the calculation (counter-terms) that can be adjusted to can- 
cel the divergent parts. The mass scale at which this procedure is applied is called 
renormalization scale. 

Renormalizable theories are those theories where the previous procedure can be applied 
and only the values of a few parameters are affected. The origin of the harmlessness of the 
quantum fluctuations in these theories can be easily understood with a method developed by 
Wilson jTj. Here, using a functional approach, we begin by considering a theory which has 
an ultraviolet cutoff scale A. Then, integrating over a momentum shell, we define the theory 
e have a new (infinitesimally) lower cutoff scale A'. Although at first this redefinition implies 
that an infinite set of new different terms can appear in the lagrangian, it can be shown that in 
a renormalizable theory only a finite subset of them tend to grow if we iterate this procedure, 
whereas the rest vanishes. The differential equations that govern the flow of the coefficients 
of the lagrangian are called renormalization group equations. Another approach to the 
renormalization group is the Callan-Symanzik equation [2J [3], which is obtained from the 
arbitrariness in the election of the scale that we use to impose the renormalization conditions 
when obtaining the parameters of a renormalized field theory. So, with the Callan-Symanzik 
equation, the renormalization group flows are obtained by looking at how the parameters 
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of the theory depend on the renormalization scale. The shifts of the coupling constants are 
reflected in one special parameter of the equation, which is called the beta function. Hence, 
by studying this parameter, we can obtain relevant physical information, as the validity of 
perturbative approach to obtain the short- or large-distance behaviour of the theory. 

Central to the physics is the idea of symmetry, which is the invariance of a physical 
system under some kind of transformation. When Quantum Field Theory was developed, 
two different types of symmetries were found: the space-time symmetries, generated by the 
Poincare group, and on the other hand internal symmetries. It was shown that both types 
of symmetries can not be non-trivially mixed [5], unless we consider fermionic symmetry 
generators [6](ie., operators that interchange fermions with bosons and vice versa). These 
generators allow us to obtain an extended Poincare algebra, which is called supersymmetry 
algebra. Since its discovery, and although it has not yet been experimentally verified, su- 
persymmetry has become one of the basic elements of modern theoretical physics. Among 
the different reasons for that, we can stand out that is a key ingredient of the theoretical 
efforts for the unification of gravity and the other forces of nature (e.g., supergravity and 
superstring theories), and it provides models that are simpler to study and quantize, as the 
symmetry between fermions and bosons implies that some "miraculous" cancellations occur 
in the calculations. 

As a renormalized quantum theory should have (if possible) the same symmetries as the 
classical one, among the relevant features that we have to maintain in a renormalization 
procedure, we have the invariance with respect to local symmetry transformations, which 
is called gauge invariance. However, the quantization procedures force us to loose gauge 
invariance in the intermediate results (for example, we have to pick up only one represen- 
tative gauge field of each gauge orbit in a functional quantization approach). To maintain 
explicit gauge invariance in every step of the calculations, the background field method [I] 
was developed. Here, the gauge field is split into two parts: quantum and background. We 
quantize the first one, which implies that we have to break the gauge invariance on it. At 
the same time, the second field is treated as a classical one, and therefore gauge invariance 
is retained in terms of it. This has relevant consequences: for example, it imposes a relation 
between the gauge coupling and background field renormalizations and allows us to obtain 
the beta function from a calculation of only the background field two-point function. 

When quantizing a gauge theory, it was found that only some renormalization procedures 
can preserve explicitly the symmetries, except for some exceptional cases called anomalies, 
being the most successful one dimensional renormalization. The key point of this method is 
to rewrite the original divergent integrals in four dimensions as integrals in D dimensions, 
being the regulator the e parameter of this continuous dimension D = 4 — 2s. As we have 
stated before, it preserves explicitly gauge invariance, making also the calculations easy 
to perform (even in the higher-loop cases). However, this procedure has some drawbacks. 
In concrete, due to the fact of changing the space-time dimension, some incongruities are 
expected to appear when applying this method to a dimension-dependent theory as can be 
a supersymmetric one. 

To solve this problem, and offer an alternative renormalization procedure that works only 
in four dimensions, differential renormalization (DiffR) was developed [7]. The basics of 
the method are to work in coordinate space rather in momentum space, rewriting expres- 
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sions that are too singular to have a well defined Fourier transform in terms of derivatives 
of less singular ones. With this prescription, it can be shown that the coefficients of the 
renormalization group equations that are satisfied by the correlators of the theory are easily 
obtained. At the same time, we stay all the time in four dimensions, making this a suit- 
able renormalization procedure when dealing with supersymmetric theories. However, one 
important practical difficult arises. Although gauge invariance is not broken, to recover the 
explicit form of this invariance in the final results we have to fix the ambiguities generated 
by the method. So, we have to impose a posteriori the Ward identities. 

This important point was solved (at one loop) by the introduction of Constrained 
Differential Renormalization (CDR) [21]. The basic idea here is to give a minimal set of 
rules to manipulate singular expressions, so that all the ambiguities of the calculations are 
fixed a priori. At the same time, all of these manipulations are required to be compatible 
with the symmetries that have to be maintained. With this prescription, it can be seen that 
the renormalized expressions directly fulfil the Ward identities without any adjustment. 

The objective of this work is to show that differential renormalization can be easily and 
applied to the renormalization of gauge theories at the two-loop level. In concrete, we will 
show that with this method we can obtain with little effort the two-loop coefficient of the 
expansion of the beta function of these theories. We have to point out that, although it is 
only fully developed for the one-loop case, to perform some of these calculations we will use 
CDR prescriptions. This is due to the fact that when imposing CDR at the one-loop level, 
the coefficients of the logarithms of the mass-scales of the two-loop renormalized expression 
get fixed a priori. No Ward identities are needed to be used. Also, we will show that 
differential renormalization clearly distinguishes between ultraviolet and infrared divergences 
as both are renormalized with different and independent mass-scales. This is not the case 
for dimensional regularization, where both types of divergences get mixed in the results, as 
they are renormalized with the same dimensional parameter e. Hence, this feature allows 
us to revisit one controversial point: the origin (ultraviolet of infrared) of the higher-order 
perturbative contributions to the beta function in supersymmetric gauge theories. Originally, 
Novikov, Shifman, Vainshtein and Zakharov obtained the so-called "exact beta function" of 
N — 1 SYM (Pnsvz) by means of instanton analysis [8], where the origin of the higher- 
order contributions was clearly infrared. However, this was questioned by Arkani-Hamed 
and Murayama [9l [10], as they were able to obtain Pnsvz in a purely wilsonian framework, 
which only depends on the ultraviolet properties of the theory. With our approach, we will 
obtain perturbatively the two- loop coefficient of (3nsvz with the advantage of having the 
UV and IR divergences clearly separated. 

The structure of the work is as follows: In the first chapter, we made a brief presentation 
of DiffR and CDR, showing also how the results of the latter can be used in two-loop calcu- 
lations. In the second chapter, we give a complete treatment of the calculation of the beta 
function of two of the most relevant abelian gauge theories: QED and SuperQED. Although 
these two theories were yet renormalized in the literature with standard DiffR, we will re- 
obtain their two-loop beta functions without imposing Ward identities. The third chapter is 
devoted to the renormalization of non-abelian gauge theories, studying the concrete models 
of Yang- Mills and SuperYang-Mills. Finally, we present our conclusions. In appendices [A] 
and [B] we made a brief presentation of our supersymmetic conventions and the background 
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field method respectively. In appendix [C]> in order to obtain the function that takes into 
account the running of the gauge parameter in the RG equations, we evaluate the one-loop 
RG equations for the quantum gauge field two-point functions of each theory that we treat. 
Finally, in appendix [D] we list some identities and calculations that are used in this work. 



Chapter 1 

Differential Renormalization and 
CDR 

1.1 Differential Renormalization 

Differential Renormalization (DiffR) [7\ is a renormalization method in real space that con- 
sists in replacing coordinate-space expressions that are too singular by derivatives of less 
singular ones. This method does not need cutoff nor explicit counterterms, although they 
are implicitly used when performing formal integration by parts. The basic idea is that 
divergent expressions are well defined for non-coincident points, but at short distances the 
amplitude is too singular and does not have a Fourier transform. Hence, to renormalize 
we are instructed by the method to replace the divergent expression with the derivative of 
a less singular one that has the same values as the original outside the origin, but with a 
well defined Fourier transform (if formal integration by parts is used with the derivatives). 
This method is especially well suited for dimensional dependent theories (such as supersym- 
metric theories), because all the time we stay in four dimensions, which is not the case for 
dimensional regularization or dimensional reduction. 

As an example consider the one-loop contribution of A</> 4 theory. The bare expression is 

A 2 

T(x 1 ,x 2 ,x 3 ,x A ) = — [<5 (4) (xi - X2)5 ( - 4l \x 3 - x 4 )[A(xi - x 4 )] 2 + (2 perms)] , 

(1.1.1) 

where A(x — y) is the massless propagator 

A(«-V) S A, = (S5)(^- 

At short distance \ does not have a well defined Fourier transform, and DiffR proposes to 
replace it for the solution of 

-1 = uG(x 2 ) x^O, (1.1.3) 
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which is 



R 



1 \nx 2 M 2 
— □ 

4 x 2 



(1.1.4) 



Both expressions coincide for i ^ 0, but the new one has a well defined Fourier transform if we 
neglect the divergent surface terms that appears upon integrating by parts the d'alembertian. 
It is in these surface terms where the counterterms hide, and by applying formal integration 
by parts [7J we are implicitly taking them into account, as we will detail later. Thus, with 
the renormalized expression we obtain 
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(1.1.5) 



A constant with mass dimension M has been introduced for dimensional reasons. It parametrizes 
the local ambiguity 

\nx 2 M' 2 \nx 2 M 2 , M' /s . 

A crucial observation is that this shift M — > M' can be absorbed in a rescaling of the coupling 
constant A [7J. This is a hint that renormalized amplitudes satisfy renormalization group 
equations, with M playing the role of the renormalization group scale. 

Let us take a closer look to the implicit counterterms that we are using in all of this 
procedure (this is discussed in [7J and with more detail in [12]). As we have stated previously, 
along with the substitution of the divergent expression with the solution of the differential 
equation, we also have to use the following formal integration by parts prescription 



d 4 x —T(x) = — 

x 4 y ' 4 



l4 . lnx 2 M\ m/ . 1 
afx (a — )T{x) = — 



x" 



, \nx 2 M 2 m/ . 
(fx — DT(x) 



x z 



(1.1.7) 

i.e., we have neglected divergent surface terms. If we made this calculation again, but 
excluding a ball B £ of radius e around the origin and keeping surface terms we have 
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The contribution of the surface integral can be found to be 

In x 2 M 2 



do^ T(x)d tl 



x, 



A7r 2 T(0)(l-\ne 2 M 2 ) + O(e) 



(1.1.8) 



;i.i.9) 



This is divergent as e — ► 0. However, this singular contribution in the 4-point function 
can be cancelled if we add to the action a local counterterm proportional to J d A x 4 (x)(l — 
\ne 2 M 2 ). Hence, as we have seen, the formal integration by parts rule is valid because we 
are implicitly using these counterterms. At the same time, we have to remark that the 
regularization method does not require us to make explicit use of them in any calculation. 
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1.1.1 Higher Loops 

Differential renormalization can be applied not only to one-loop diagrams, but to multi-loop 
expressions. In general, new scales appear corresponding to the renormalization of the dif- 
ferent subdiagrams that form the total expression. From the various types of sub divergences 
that can occur in a typical higher-loop diagram, we will take a closer look to one of them, 
where independent scales are neatly seen to appear at each stage: nested divergences. As an 
example we consider the following amplitude which, in principle, could form part of a bigger 
one: A(x — y)I 1 (x — y), where is 

l\x-y) = J d 4 uA xu A 2 yu . (1.1.10) 

It corresponds to a diagram that looks as follows 



Figure 1.1: Two-loop diagram with nested divergences. 



As can be seen, divergences occur whenever 2 points come together. We can renormalize 
them starting from the most inner one, and proceeding recursively 



A xy J d 4 uA xu [Al u ] R 



R 



1 



\n(y - ufM\ 



4(4vr 2 ) 4 (x - y) 2 
1 ln(a; - y) 2 M\ 
_4(4tt 2 ) 3 (x - y) A 



d A u 



-□ 



(x — u) 2 (y — uY 



R 



J R 

\2 A/r2 



32(4vr 



2^3 



-□ 



ln 2 (x - y) 2 M 2 + 2 \n(x - y) 2 M 2 
(x - y) 2 



where, in going to the second line, we have integrated by parts the d'alembertian and made 
use of □ /^j^p = S(x — u). We observe the appearance of an independent scale associate 
with each renormalization step. 

A systematic implementation of differential renormalization to all orders in perturbation 
theory was presented in [13] . The basics of the method are the separation of the divergences 
in two groups: one corresponds to divergences derived from two points collapsing, and the 
other to three or more points simultaneously closing up. For the first one the singularity 
is replaced with the renormalized form (once the derivatives are pulled in front), whereas 
the other one can be shown to be recursively written as two-point function problems of the 
first type. This procedure follows the BPHZ renormalization program and guarantees that 
differential renormalization maintains unitarity and it can be applied consistently (fulfilling 
locality and Lorentz invariance) to all orders [IB"]. 
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1.1.2 Massive theories 

Differential renormalization of massive theories has been studied in [7J [TT] . The appearance 
of a bare mass does not interfere with the method, since DiffR is related to short- distance 
singularities and masses only change the long-distance behaviour of the correlators. Although 
in this work we will only deal with massless theories, we will give briefly as an example how 
this procedure works with massive A0 4 . The propagator of a particle of mass m is 




Ki(V m 2 x 2 ) 



X1.12) 



where K\ is a modified Bessel function. Let us now consider again the 4-point function 
contribution; in this case is clear that the expression we have to renormalize is 



— —K\{\/ m 2 x 2 ) 



x- 



(1.1.13) 



We have to solve the massive generalization of the differential equation (jl.l.3p . which has 
a solution of the form of 




K\ (V m 2 x 2 ) 



Ft 



'□ - 4m' 



M 2 



Kq(V m 2 x 2 )Ki(Vm 2 x 2 ) 



+7r 2 ln— 6(x) 



1.1.14) 



where M = 2M/7 and 7 is the Euler constant. The general solution has a contact term 
which depends on a new mass parameter M. This guarantees that in the limit where m — > 
the renormalized expressions for fjl.l.4|) and (11.1.141) coincide. 



1.1.3 IR divergences 

DiffR can be also applied to expression with IR divergences [15], i.e. expressions that exhibit 
a divergence for p^ — > 0. The idea is to apply a dual version of Differential Renormalization 
to such quantities 

1 \np 2 /M? R r/ , . 

,,^ = -4^ p 2 (1.1.15) 



We have defined for convenience M/_r = 2Mm/ / yE, where je is Euler's constant, and dis- 
tinguished the IR scale from the UV one. As DiffR is an implementation of Bogoliubov's 
R operation (an operation that yields directly renormalized correlation functions satisfying 
renormalization group equations), in momentum space this is an explicit realization of the 
so-called R operation that subtracts IR divergences. Again, diagrams with IR sub divergences 
are treated according to a recursion formula [161 E] analogous to the UV one. 

Now, when going to higher loops a new effect will show up, namely the coexistence of 
UV and IR scales. Let us start by examining a prototypical diagram where such divergences 
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Figure 1.2: Two-loop diagram with UV and IR divergences. 



arise in a A0 3 type theory (figure [L2|) . The associated amplitude that we have in this case 
is of the form G(x — y) = A(x — y)I°(x — y), where 1° is 



I\x-y) 



d\d 4 v A xu A yv A 2 uv 



;i.l.l6) 



We begin by renormalizing the inner UV divergence as 
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;i.i.i7) 



In order to renormalize the IR divergence we have to pass to momentum space 
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We explicitly observe that the IR singularity at p — > involves an UV scale M. 

Since UV and IR overall divergences are local in coordinate and momentum space, re- 
spectively, the R and R operations commute, and one can define an operation R* = RR to 
renormalize both UV and IR divergences [TBI E] • The fact that the UV and IR renormal- 
izations decouple means that the UV and IR renormalization scales should be independent. 
This is a non-trivial point that in DiffR can be achieved by a careful adjustment of the local 
terms involving both scales^- 

As we have to guarantee that the IR renormalization commute with a rescaling of M, we 
have to fulfill the relation 
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1 IR DiffR, was investigated in [18] where it was concluded that the combination of UV and IR DiffR was 
inconsistent, as the results depended on the order in which integrations were performed. According to |19| . 
however, this corresponds to the natural arbitrariness of the IR renormalization, and this author has actually 
proposed in [5U] a consistent version of DiffR that deals with both UV and IR divergences. Our approach 
here will be closer to the original version of DiffR 
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If we consider the usual expression for the renormalization of (Inp 2 /M 2 )/p 4 



\np 2 /M 2 
p 1 
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1 \n 2 p 2 /M 2 + 2\np 2 /M IR 



we find that the left hand side of (11. 1.191) is 
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whereas the right and side has the form of 
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So, the second expression differs from the first one in a local term in momentum space, that 
will be a non local one in position space. Hence, in order to fulfill (I1.1.19p . we propose the 
following minimal solution that does the job [15] 

\np 2 /M 2 ] 1 -\n 2 p 2 /M? R + 2\np 2 /M? R (l + \np 2 /M 2 ) 

"77 Do 



p l 



J R 



(a m \n— ^ 
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(1.1.23) 

This expression differs from the usual one by scale-dependent local terms proportional to 
In 2 M 2 /M 2 R (apart from the explicit local terms with coefficients am and bm). It should 
be used whenever the "new" scale is to be treated as independent from the "old" one, for 
consistency of the loop expansion. It has to be noted that when we consider a purely UV 
expression as in (11.1. lip , we do not have to take care of all of this, because the extra term is 
cancelled in the RG equation when we derive wrt. the UV scales. With this, the result for 
1° (not taking care of the local terms am and bm) is 



l R 



X 



y) 



32(4tt 2 ) 3 
(x - yf 



32(4tt 



2\3 



d A p \jp- 



d 4 p 



\n 2 p 2 /M 2 R + 2 lnp 2 /M? R (1 + lnp 2 /M 2 



-ip(x-y) 



P 



In 2 p 2 /Mf R + 2 hip 2 /M] R (1 + lnp 2 /M 2 



-ip(x-y) 



P 



That in position space is 



IX 



[In 2 x 2 M} R + 2 \nx 2 M] R (l - \nx 2 M< 



[1.1.24) 
(1.1.25) 



32(4vr 2 ) 2 

Observe that the UV scale M only appears in (11.1.231) . and hence in the above expression for 
J°, in single logarithms. This is fine, for double logarithms of M are expected to appear only 
when the bare expression contains both a UV subdivergence and a UV overall divergence. 
Finally, once we have obtained J°, we can straightforwardly evaluate G to be 

In 2 x 2 M 2 R + 2 \nx 2 M 2 R (l - \nx 2 M 2 ) 



G(x) 



1 

32(4vr 



2\3 



;i.i.26) 
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1.1.4 Symmetries with DiffR 

One of the important properties that is required to every sensible renormalization procedure 
is not to break gauge symmetry when it is applied to a gauge theory. For DiffR it is found 
that gauge symmetry is preserved as long as Ward identities can be always satisfied with 
the renormalized amplitudes (except anomalies). However, at the same time we find that 
we have to make always explicit use of these identities to fix all the ambiguities that have 
appeared in the calculations; in particular with Ward identities we relate the different scales 
that we have to use when renormalizing different amplitudes related by a symmetry(i.e. we 
fix a renormalization scheme). 

As an example of this, consider the case of the one-loop renormalization of the photon 
self-energy in QED [f4j. The bare expression is 



IV(x)U re = -4e 2 [2(^A(x))(^A(x)) - 6^{d p A{x)){dpA{x))] 



'1.1.27) 



and renormalizing, this becomes 

„2 



I\-uM\r 



12tt 4 

e 2 
12tt 4 



R 



l n _ lnx 2 M 2 1 



□ □ 



In x 2 M\ 



- 8(^/! 2 <5( 



x 



1.1.28) 



In this expression we have renormalized with an independent scale the logarithmic diver- 
gence (Mi) and the quadratic one (M2). At the same time, related to the latter, we have 
added a possible local term with a parameter of mass dimensions (/i). The Ward identity 
satisfied by this expression imposes that this has to be transverse, so that we have to fix 
M x = M 2 and /i 2 = 0. 

When going to higher loop computations, the Ward identities play a non-trivial role, as 
they influence part of the divergences that are obtained in the next step. The reason is that 
these identities relate all the relevant mass scales found. So, they allow us to write the one- 
loop renormalized subdiagrams that made up the two-loop expressions in terms of the same 
scale, say M, and fixed local terms, which then are promoted to logarithms of the scale. Going 
back to the example at two loops that was solved in (11 .1.1 II) . suppose that after imposing 
the Ward identities we have the inner divergence written as — |nlnx 2 M 2 /x 2 + aS(x), with 
a a fixed coefficient. Thus, in the two- loop expression we find 



A„, / d 4 u A XU [A 2 } R 



yul 



J R 



1 \n(y - u) 2 M 2 + a 



_4(4vr 2 ) 4 
1 

~32(4tt 2 ) £ 



-□ 



(x-vY ±R 

In 2 x 2 M 2 + 2 In x 2 M' 2 + 2a In x 2 Mf 
x 2 



4.1.29) 



where M' 2 and M' 2 ' are two-loop scales. Hence, as we have anticipated, the one-loop Ward 
identities have fixed the coefficients of the logarithms of the scales in the two-loop final 
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expression. Concerning the new two- loop scales M' 2 and M'^ , it is clear that both can be 
set also equal to M modulo a local ambiguity that will depend on quotients M' 2 /M and 
M'2/M (like in (11.1.61) ). Again, use of the Ward identities would set these quotients to 
certain computable values. In other words, after use has been made of the symmetry, in 
the two-loop expression the only scale that remains can be chosen to be M and sits only 
inside the terms with logarithms whose coefficients were determined from the one-loop Ward 
identities. This observation is at the heart of the present work an permeates implicitly all 
the calculations contained in it. So we repeat it here for full clarity: if one is interested in 
computing a physical amplitude at two loops, a concrete value of the local terms is essential 
and use Ward identities at two loops is unavoidable. If however, as is the case of the present 
work, one is looking for the RG equations, then all the relevant information on the scale M 
resides in the terms with logarithms, whose coefficients only need one-loop Ward identities 
to be fixed. 

1.2 Constrained Differential Renormalization 

Constrained Differential Renormalization (CDR) was developed in [2~T | |2"2" ||2"3~II2~T] to avoid the 
necessity of imposing Ward identities in each calculation to fix the renormalization scheme, 
as we have seen in the previous section. The idea is to give a procedure that allows us to 
fix the scheme a priory. Central to the fulfilment of the Ward identities (and the action 
principle, from which they can be derived) is that the application of the kinetic differential 
operator to some propagator line inside a Feynman graph is equivalent to the contraction of 
the line to a point [2"Tj . This statement is guaranteed to hold if we apply the following set 
of rules 

1. Differential reduction 

• Functions with singular behaviour worse than logarithmic are reduced to deriva- 
tives of (at most) logarithmically divergent functions without introducing extra 
dimensionful constants. 

• Logarithmically divergent expressions are written as derivatives of regular func- 
tions, introducing one single constant M, which has dimensions of mass and plays 
the role of the renormalization group scale. 

2. Formal integration by parts. We do not take care of the divergent surface terms that 
appear when we integrate by parts. Related to this, differentiation and renormalization 
must be two commutative operations: let F an arbitrary function, then = d[F] R . 

3. Renormalization rule of the delta function: 

[F(x, xi, . . . , x n )5{x - y)] R = [F(x, xx, . . . , x n )] R 5(x - y) (1.2.30) 

4. Validity of the propagator equation 

[F(x, x u . . . , x n )(0 - m 2 )A m (x)] R = —[F(x, x lt . . . , x n )6(x)] R (1.2.31) 
where A m is the propagator of a particle of mass m and F an arbitrary function. 
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The upshot is a basic set of renormalized expressions (basic functions) with different 
numbers of propagators and various differential operators acting only on one of them, in- 
volving a single scale M . Therefore the CDR program amounts to the following two step 
operation: 

• Express the Feynman diagram in terms of these basic functions performing all the 
index contractions (this is an important point, because CDR does not commute with 
index contraction) and, by means of the Leibniz rule, moving all the derivatives to 
make them act on one of the propagators. 

• Replace the basic functions with their renormalized version. 

Let us now obtain as an example some of these functions. Consider the one-point basic 
function A(x)5(x) (this corresponds to the one-loop correction to the two-point function in 
A$ 4 theory). Power counting and the locality of the expression implies that the most general 
renormalized value that we have for this is 

[A(x)5(x)] R = (cn + v 2 )5(x), (1.2.32) 

where /i is a mass-dimension constant and c an adimensional constant. However, rule 1 
implies that fi = 0. Now, considering [A(x)5(x)] R 5(y) and using rule 3 we find 

[A(x)5(x)] R 5(y) = [A(x)5(x)5(x + y)] R =[A(x)5(x)] R 5(x + y) 

(1.2.33) 

and integrating over x we arrive to 

6(y) J d 4 x [A(x)6(x)] R = [A(y)S(y)] R . (1.2.34) 

Finally, with this result rule 2 implies that c = 0. Proceeding in a similar way we find that 
all the massless one-point functions in CDR vanish. 

As two-point function examples, we will consider Ad^A and ADA. In the first case we 
have to apply the Leibniz rule to find 

[A(x)d f _ t A(x)) R = d lt [A(x)A(x)] R -[(d ll A(x))A{x)] R 

L rA9/ 1 „ lnx 2 M 2 

= -d,[A \x)} R = -^-^8^—^- , (1-2-35) 

where we have used the result of (jl.l.4p . If we study now ADA, we have only to use rule 4 
to arrive at 

[A(x)DA(x)] R = -[A(x)5(x)] R =Q. (1.2.36) 

Here we present as a summary the most relevant CDR identities that are used in this work. 
We only list the massless examples, although a complete list including massive propagators 
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can be found in 



r a 21 / \ 1 WM 2 

= "4(4^°-^ 

[t _ . . . . 1 _ lnx 2 M 2 

' •,{dpd v -8 ia ,U)8{x) 



288tt 2 

[AnA] fl (x) = 0. (1.2.37) 

CDR can be applied to more than two propagators. In particular, when dealing with 
three propagators, defining T[0] = A AO A, we can find the following relation to hold when 
making a decomposition into trace and traceless parts [211 123] 

T R [d,d v ] = r 8 ^ - i W + J^P] - i2^M(« ■ (1.2.38) 

When using other gauges different from Feynman gauge, some bare expressions are writ- 
ten in terms of a quantity we define as A(x) = 4 ^ n2 ^ lnx 2 s 2 , where s is an irrelevant constant 
with mass dimension. For this structure, CDR prescribes [21] 

r . x ~i / \ 1 lnx 2 M 2 



CDR has been checked in abelian and non-abelian gauge theories [2~Tj [21] and in super- 
symmetric calculations [251 ES]- As an example of its use, we will re-obtain the one- loop 
renormalization of the photon self-energy of QED that we have renormalized with DiffR in 
the previous section. From the bare expression (I1.1.27P we apply rule number 2 to write it 
in terms of the CDR basic functions as 

lWar)|fl = -4e 2 [2^(A^A)-2A^A-^ /3 (A9 /3 A) + ^(ADA)] iJ . 

(1.2.40) 

Now, we have to replace each basic expression with its renormalized value, and straight- 
forwardly we arrive to 



lWa;)|fl = (dpdv - 5^,0) 



e 2 lna; 2 M 2 e 2 r/ s 



_3(4tt 2 ) 2 x 2 36?r 



1.2.41) 



As we have remarked, CDR has fixed all the ambiguities a priori, obtaining a direct final 
result that is transverse, as it has to be to fulfill the Ward identity. 

Finally, it is worth to mention that this method is equivalent to a momentum-space 
regularization method defined also in four dimensions: Constrained Implicit Regularization 
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(CIR). Implicit Regularization [28l [29] is a regularization method based in the assumption 
of a regulating function as part of the integrand of divergent amplitudes, and the extension 
of the properties of regular integrals to regularized ones. As in differential renormalization, 
this procedure generates arbitrary parameters that with CIR are fixed a priori [30J . 

1.3 Two- loop uses of one- loop CDR results 

As we have seen, one of the drawbacks of DiffR is the plethora of scales that pop up at each 
step of the calculation. In symmetric theories, at fixed order in the perturbative expansion, 
these scales should reduce to a single one upon use of the Ward identities. In the previous 
section we have explained how CDR paves the way to this reduction of scales at the one-loop 
level. So far, CDR has not been fully developed at loop-order higher than one, and therefore 
it is not useful for computing, say, scattering amplitudes. However, as is mentioned at the 
end of section II. 1.41 as long as we are interested in the RG equations, all that we need are the 
terms with logarithms, and to obtain them the knowledge of the local terms at one loop-level 
is enough. Hence, we will discuss both the way the logarithms are generated from one loop 
to the next, and the implementation of the CDR rules in such diagrams [3T] . 

1.3.1 Nested divergences 

This case is particularly simple because CDR can be applied in a systematic way. Starting 
from the "inner" divergence, its regularization according to CDR gives an expression with 
logarithms of a single scale (lnx 2 M 2 ) and fixed local terms. The one-loop Ward identities are 
fulfilled. In the next step, when tackling the outer part of the diagram, a simple logarithm 
like the one shown above is promoted to an expression of the form In 2 x 2 M 2 + Clnx 2 M' 2 , 
with C a calculable coefficient and M' a two-loop scale; at the same time, the local terms 
that multiply outer divergences will produce additional logarithms of new scales. CDR does 
not yet prescribe what the different two-loop scales should be; hence, we may take all of 
them the same, and equal to M, at the price of leaving undetermined local terms which are 
irrelevant when obtaining the RG equation. 

This simple scheme has some subtleties when considering diagrams with indices because, 
even at one-loop, index contraction does not commute with CDR. Therefore, the correct or- 
der is to first insert into the outer diagram the non-renormalized expression for the "inner" 
one-loop diagram, perform all the index contractions, and then renormalize. This crucial 
observations looks as the first one in a list of rules that eventually would setup the imple- 
mentation of CDR at higher loops. 

Let us consider now the two-loop example discussed in sections II. 1.11 and 11.1.41 If we 
had imposed CDR in the first step, Mi = M is the only scale generated upon renormalizing 
the most internal divergence, and the local one-loop ambiguity will be fixed to zero, as can 
be seen from f ll.2.37p . We express this by stating that the renormalization of J 1 (a; — y) — 
J d 4 «A IU A^ according to CDR rules is given by 




1 In x 2 M 2 
4(4vr 2 ) 2 x 2 



(1.3.42) 
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Once we have this, to renormalize the complete two-loop expression A/ 1 , we have to apply 
usual differential renormalization and set, modulo local terms, the two-loop mass-scale M' 2 = 
M. So, we arrive to an expression of the form of 



[A/'] 



R 



X = 



32(4tt 2 ) 3 



In 2 x 2 M 2 + 2\nx 2 M 2 

□ 5 + ••• > 

X 2 



1.3.43) 



where . . . stand for the two-loop local terms that we are not taking into account. Notice that 
in the rest of the work (unless explicitly stated otherwise) ... in a two-loop renormalized 
expression like the one shown above will have the same meaning: local terms not considered. 
With this procedure we have renormalized all the different structures made of I 1 that we have 
encountered in our calculations. Apart from the previous one, we have found the following 
relevant expressions 



[Ad,dJ l ] R (x) 



1 



64(4tt 2 V> " 
1 



n \n 2 x 2 M 2 + \nx 2 M 2 
d„u + . . . 



ar 



(1.3.44) 



96(4tt 2 ) j 



In 2 x 2 M 2 + Unx 2 M 2 



1 ln 2 x 2 M 2 + ^lna; 2 M 2 

- -5„.„nn -3 



[An/ 1 ] R (x) 



4^ 
1 

32(4vr 2 )' 



In x 2 M 2 
■□□ \- ... 

x 2 



+ ... 



1.3.45) 
1.3.46) 



To obtain each of these results, we only have to consider the CDR renormalization of I 1 , 
and apply afterwards usual DiffR, setting all the two-loop mass scales equal to M. Let us 
illustrate the simplicity of the procedure with Ad^I 1 . Given the renormalized form of I 1 we 
find 



R 



, s CDR 
[X] = 



1 



4(4tt 2 ) 3 

DiffR 1 



1 . \nx 2 M 2 

—zO, 



X 



2 M 



J R 



16(4tt 2 ) 3 
1 



2\nx 2 M 2 



x 1 



J R 



64(4vr 2 ) 



2^3 V 



n \n 2 x 2 M 2 + \nx 2 M 2 
d„n + . . . 



x^ 



(1.3.47) 



1.3.2 Overlapping divergences 

Diagrams with overlapping divergences are more complex as it is sometimes difficult to 
recognize the one-loop sub divergences that need to be treated with CDR to start with. Our 
approach will be to obtain through different methods (that we will explain later in detail), a 
list of renormalized two-loop integrals with overlapping divergences, where in each calculation 
one-loop CDR rules have been maintained in every step. Although this list is restricted to 
integrals with at most four derivatives acting on the propagators and two free indices, it is 
found to be very useful, as serves as a basis that we can use to express the renormalized 
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overlapping contributions to two-point functions in theories with derivative couplings at two 
loops. As we detail in appendix [Bj these two-point functions are what we need to obtain the 
beta function if we use the background field method. This list will be applied in our work to 
renormalize and obtain the two-loop beta function of (Super) QED and (Super) Yang-Mills. 
We use the conventions of z = x — y and <9 M = We also define H(x — y) = H(z) as 

H[0 1 , 2 ; 3 , O a ] = j d 4 ud 4 v (0$A xu )(OZA xv )(O y 3 A yu )(OiA yv )A uv , (1.3.48) 



with Oi a differential operator. 

1,1] 

H R [d„ 1-1,1} 
H R [l,d x ; l,d x ] 
d x H R [l,d,- l,d x ] 
dxH R [l,l- dxdv,l] 
H R [l,d x] d x d^l) 



H R [d„d x ,d x ■ 1,1] 



d x H R [l,d^; d v d x ,l] 



d x H R [l,d x ; 8^,1) 



d x H R [l,d x ; l,d„d u 



6^(3) 
(4tt 2 ) 4 
3^(3) 



A = a A 



16(4tt 



2^2 
1 



(d,A) = -d,A 

In z 2 M 2 
□ = + . 



16(4tt 2 ) 3 z 2 



32(4tt 



d„a- 



1 \nz 2 M 2 



2\3 t 1 



+ 



32(4vr 



2^3 



d v U- 



±ln 2 z 2 M 2 + f \nz 2 M 2 



+ 



32(4tt 2 ) 



1 \\n 2 z 2 M 2 -l\nz 2 M 2 

d„U 8 -i + 



2^3 ^ 



d„n- 



z* 

±ln 2 z 2 M 2 - In z 2 M 2 



32(1 rr 2 * 
1 



+ . 



32(4vr 2 ) 3 

+S uv an~ 1 



±ln 2 z 2 M 2 + \\nz 2 M 2 
dpd,, □ ^ ^ 

\\nz 2 M 2 



z c 



J fiu L 
1 



+ ... 



-\\nz 2 M 2 



32(4tt 2 ) 3 
1 



+ . . . 



2^3 



9A n — » 

z 2 



±lnz 2 M 2 



32(4tt 

lln 2 ^ 2 M 2 + |lnz 2 M 2 
+ V □ □ J ^ 



2 X 



+ . . . 



2^3 



llnz 2 M 2 



d d □ — 

a; 2 



32(4tt 

±ln 2 2 2 M 2 + |lnz 2 M 2 
s- 5 



2 Z 



+ . . . 



(1.3.49) 
(1.3.50) 
(1.3.51) 
(1.3.52) 
(1.3.53) 
(1.3.54) 
(1.3.55) 



(1.3.56) 
(1.3.57) 



(1.3.58) 



(1.3.59) 
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H R [1,1; 9^,1] 
d x H R [l,l; d x d u ,d,] 
d x H R [l,l; d^dx] 



if*[l,0,A; d v d x ,l] 



Hn 2 z 2 M 2 + %\nz 2 M 2 



32(4vr 2 ) 3 
1 



+ 



32(4vr 2 ) 3 



iln 2 z 2 M 2 + l\nz 2 M 2 



1 



Mn 2 z 2 M 2 + f In z 2 M 



2 



32(4vr 2 ) 3 " v 
1 



+ 



;i.3.60) 
'1.3.61) 
[1.3.62) 



32(4vr 
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±ln 2 z 2 M 2 - 4WM 2 

□ ^ 5-36 



-^lnVM 2 -§lnz 2 M 2 
+<W □ □ s-^ 



+ . 



;i.3.63) 



32(4tt 



2\3 



Mn 2 2 2 M 2 + §ln^M 2 

d»d v U$ ^ 



-Mn 2 z 2 M 2 -^\nz 2 M 2 

+<W D D — 5-^ 



+ . 



(1.3.64) 

These integrals are obtained basically applying two properties: 

• Integral relations presented in appendix [Dl These exact relations allow us to put some 
of the integrals in terms of others that have an explicit d'alembertian acting on one of 
the propagators. Once we have done that, using DA = —5 we can put these integrals in 
terms of the previously defined I 1 . Then, we can straightforwardly apply the procedure 
for nested divergences that we have just presented in the previous section]! 

• The decomposition into trace part, traceless part and fixed local term imposed by CDR 
to T[d„d v ] as flfflgp . 

As in the previous section, let us illustrate the procedure with an explicit example. 
Considering integral (11 .3.531) . this can be evaluated with both methods. First, we will make 
use of integral relation (ID. 1.21) and put this integral as sum of different integrals that have 
the divergences nested 



d*x 



d 4 ud 4 v A xu A xv {d\dlA yu )A yv A l 



+ J d 4 ud 4 v A w4 A OT (DA wu )(^A tfl> )A 14W 
-\dldl J d 4 ud 4 v A.A.^A^A, 



;i.3.65) 



2 Also this is the reason why we have not listed here the cases where the differential operator is a 
d'alembertian. For example, it is obvious that if [□, 1 ; 1, 1] = —A/ 1 . 
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Now, we have to apply as usual OA = —5 and rewrite these integrals in terms of I . 
Note that the third integral can be easily shown to be finite, and its value is obtained in 
appendix ID. II to be a/Ad v (OA) with < 



67T 4 g(3) ' 

(4tt 2 ) 4 - 



01 j A xu A xv{m A M 



-\d v {AI l ) + \iAdJ 1 ) + \d v {OA) . 



;i.3.66) 



Applying the results found in section 11.3.11 for the J 1 expression, the renormalized value 
of this is 



±ln 2 z 2 M 2 + '±\nz 2 M 



32(4vr 



2\3 



nd v O- 



2 /1#2 



+ . . . 



1.3.67) 



where . . . stands for the finite terms that we are not taking into account and z = x — y. 

We can also obtain this integral making use of the CDR relation ( 11.2.381) and perform a 
trace-traceless decomposition of {d v x dlA yu )A yv A uv as 



-dl \ d A ud A v A xu A xv (OA yu )A yv A uv 



+df I d%d\ A XU A XV 



{d\dl - -5 Xu o)A yu 



J d 4 ud 4 v A xu A xv 5(y - u)S(y - v) 



128 

-^d v (AI 1 ) R - ^ 2 d„Al + d\I\ u R 
1 \hi 2 z 2 M 2 + \\iyz 2 M 2 

-d„n* ^ + d x l. 



32(4tt 2 ) 3 



\-L\u R , 



(1.3.68) 
(1.3.69) 



where d\I\ u r is the traceless part, that is finite. As we can see, both results agree. Although 
in this example we can perform the calculation with both methods with the same effort, with 
other integrals the situation is different, and we shall have to study each case in order to 
choose the best one. The explicit evaluation of all the integrals is presented in section ID. II 
of appendix [D] . 



CHAPTER 1. DIFFERENTIAL RENORMALIZATION AND CDR 



Chapter 2 

Abelian QFT applications 



In this chapter we apply the ideas and methods we have just presented to two of the most 
relevant examples of abelian gauge theories: QED and its supersymmetric extension, Su- 
perQED. Although both theories have already been treated in [HI [37] using DiffR, we will 
show that our procedure simplifies the calculations, avoiding the use of Ward identities. 

2.1 QED 

QED is one of the simplest examples of a gauge theory, as the gauge symmetry group is an 
abelian one, U(l). Hence, is a good theory to start with, as we can clearly see all the key 
points of our renormalization procedure. 

2.1.1 The model 

We use the same conventions as [TJ] . The d = 4 massless QED lagrangian is 

C = -F^F^ + ^d^ + ieA^if) , (2.1.1) 

where if) is the fermion field, is the U{1) gauge field and F^ u is the field strength made 
up with as F^ v {x) = d fl A u (x) — d u A fl (x). The 7 matrices satisfy the Clifford algebra 

With w an infinitesimal parameter, the QED action is invariant under the following U(l) 
transformations 

Afj, — > A^ — -d^w 

if) — > if) + iipw 
$ -»• if) - iipw . (2.1.2) 

Hence, when quantizing QED we have to take care of this invariance. We have an infinite 
number of different gauge field configurations (those obtained through gauge transformations 
from a given one) that correspond to the same physical state. In a path integral approach, 
we want to integrate only over the relevant gauge field configurations; hence, we have to pick 
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up only one field from each gauge orbit. To accomplish this there is a well-known procedure 
[32l |4"7] which implies that we have to add to the action a gauge fixing term that depends 
on a new parameter a and possible auxiliary fields (Faddeev-Popov ghosts fields, that in the 
concrete case of QED and SuperQED are not relevant). Different values for a correspond to 
different gauge choices. In particular, in our calculation we will use a = 1 (Feynman gauge). 
The complete lagrangian is then 

C = -F^F^ + ^-^Arf + ^id^ + ieA^ . (2.1.3) 

With this action the gauge field and fermion propagators (in a generic gauge) are 

1 



A^(x-y) = — "—{8 llv U + {a-l)d l ,d v )\n{x-yfs 2 

l07T z 

S(x-y) = -j x d x A(x-y), (2.1.4) 



where s stands for an irrelevant constant with mass dimension. Also, considering the ex- 
pansion of the effective action, we write the terms corresponding to the vacuum polarization 
II Mi , and fermion self-energy £ as 



eff 



x-y) 



AM 



l - J d 4 xd A y A^x) (Yl - ^ d„d v - S^a) 5^(x -y)- n^( 
+ / d 4 xd 4 y 4){x){^d^ {A \x —y) — £(x - y))ip{y) + ... (2.1.5) 



Background field method 

Fixing the gauge is necessary in order to quantize the theory, but has the drawback to make 
us lose explicit gauge invariance in the intermediate results. In order to avoid this, the 
background field method was developed [I] . As the method is detailed in appendix \E[ we 
only briefly outline it here. The key point is the splitting of the gauge field in two parts: the 
quantum and the background fields (A^ and respectively) 

A^A^ + B^. (2.1.6) 

We can use A^ as the integration variable of the partition function, which implies that the 
gauge has to be fixed only for this field. As a result, we retain explicit gauge invariance in B^. 
Along with this, as is shown in appendix [Bj this procedure has other relevant consequences: 
the coupling constant and the background field renormalizations are related, which implies 
that the beta function can be obtained only from the two-point function contribution. Hence, 
we have a background effective action of the form 

Te ff [B] = i I d'xd'y B,(x) [{d,d u - S^U) 5^(x - y) - TL* u B (x - y)] B u (y) + . . . , 

(2.1.7) 
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and our aim is to calculate the two-loop expansion of the two-point 1PI function Il^f. Notice 
also that, as the infinitesimal gauge transformation of the unsplit theory does not depend on 
the quantum gauge field, we can choose in the split theory the usual gauge fixing condition 
G = d^A^. So, we have a split lagrangian of the form of 

C = iF^F^ + lB^B^ + ^id^Arf + tPYid^ + ieA^ + ie^B^ , (2.1.8) 
4 4 la 

with B^ v = d^By — d v B^. From this lagrangian, we have two relevant interaction vertices, 
which are shown in figure 12.11 




Figure 2.1: QED interaction vertices. Thick wavy lines represent background external fields 
and thin wavy lines correspond to quantum gauge field propagators. Solid lines represent 
fermion fields. 



2.1.2 One- loop level 




Figure 2.2: One-loop QED diagrams. 



We first consider the one-loop renormalization of the background photon and fermion 
self-energies. The diagrams that correspond to these contributions are those of figure 12.21 
where thick lines represent B^ fields, and thin ones correspond to A^. The correction to 
the background self-energy is what we need to obtain the one-loop coefficient of the beta 
function, although it will not be used in two-loop calculations. On the other hand, the 
fermion self-energy will be not relevant for the one-loop beta function study, but it will be 
used afterwards as a one-loop insertion in a two-loop diagram. 
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Photon self-energy 

We consider now the one- loop contribution to the photon self-energy H^. The expression 
for this diagram is 

tf!?M(x-y) = -{iefTr [ 7 ,7 A ^A 7 ^^A] 

= -e 2 Tr [ 7m 7 A 7,71 (0 a (A&A) - A<9 A <9<j A) . (2.1.9) 

To proceed with the CDR program, we first have to perform all the index contractions, 
writing this diagram in terms of the basic CDR functions (ll.2.37p . To do that, we have to 
apply in 12.1.91 the following Clifford algebra result 

Tr [7 m 7a7i/7<t] = 4(5 m a<W - <V<$Aa + <W<^0 ■ (2.1.10) 

which allows us to find the fully expanded expression as 

K*d)( x ) = -4e 2 [2^(A^A)-2A^A-^ A (A9 A A)-^ADA] . (2.1.11) 

Now, CDR renormalization of this expression entails only to replace these basic functions 
with their renormalized values. The full one-loop renormalized background self-energy is 
given bv [HI 121] 



127r 2 (4vr 2 ) x 2 36vr 



(2.1.12) 



As was guaranteed by the use of CDR, this result is transverse (as is required by the 
Ward identities), and has the ambiguity (local term) fixed. 

Fermion self-energy 

We will also consider the renormalization of the fermion self-energy E. In a general gauge, 
the bare expression for this diagram is 

S (1) (x) = e 2 7 M A^(x)7 A 9 A A(x) 7 , • (2.1.13) 

By making use of the expression for the photon propagator in a generic gauge (12.1.41) , we 
can replace the basic functions by their CDR values and obtain the renormalized fermion 
self-energy as [HI [2TJ 



, x 1 _ lnx 2 M 2 , 1 ^ lnx 2 M 2 1 or/ 

(2.1.14) 



2.1.3 Two-loop level 

Now we proceed with the two-loop case. There are two relevant graphs with external back- 
ground fields. Diagram (a) has the divergences nested, whereas diagram (b) has overlapping 
divergences. 
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Figure 2.3: Two-loop QED diagrams. 



Diagram (a) 

The expression for this diagram is 

K B (2a)(x-v) = -{ief J d\d\Tr[ 1 , 1 \-dlA xu )^ 1 \u-v)Y(-d:A vy ) 1 ^ 

x(-dy a A yx )]. (2.1.15) 

where Y,^ is the one- loop fermion self-energy. In the following we will restrict ourselves to 
Feynman gauge, as the term that takes care of the running of the gauge parameter in the 
RG equation will be shown not to be relevant for the two-loop beta function [2]. We will 
discuss this later in detail, when we apply the RG equation. So, in this gauge, the bare 
fermion self-energy is 

= -2e V [A<9 A A] (x) . (2.1.16) 

As we stated in section 11.21 CDR imposes a strict order to the operations of index 
contraction and renormalization: first all the indices should be contracted, and only after 
that we can renormalize. Inserting here the bare fermion self-energy we are keeping this 
order. 

Expanding the expression of n^f/ 2 ^ we find 

nSf (2a) (ar-y) = -2e 4 Tr[7^7A7p7 e 7,7 ( x](^A :E3/ )^^| d\d 4 v A xu A yv (A uv d u p A uv ) . 

(2.1.17) 

In order to simplify the notation, we define an integral expression of the form 

I°(x ~y) = J d 4 ud 4 v ^(^AJ . (2.1.18) 
Then, with standard Clifford algebra, we can write ( 12.1. 17ft as 

n£f (2 = e 4 [-32{d,A)d x dJl + m, u {d a A)d x dJl + i6(<9 M A)n/° - %8, v {d p A)ui»] . 

(2.1.19) 

The renormalization of 7° is studied in section ID. 2. 21 of appendix [Dl It is found there 
that this integral expression verifies d^I® R = — and = —^9^1^. Thus, we need 
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only the renormalized value that we have previously obtained for I . We have finally 



p 4 

T1 BB (t\\ 

A V(2a)W|fl 24(47T 2 ) 3 
+ . . . 

Diagram (b) 



ln 2 x 2 M 2 -l\nx 2 M 2 _ In 2 x 2 M 2 + | lnx 2 M 2 



^ ^ + 



(2.1.20) 



This diagram, opposite to the previous one, has overlapping divergences. Following the 
procedure presented in section 11.3.21 we will express this in terms of the integrals listed in 
that section. We begin by considering the bare contribution 

ngWar-y) = -(ie)* J d*ud 4 v Tr [ 7(1 (t^A 11( )/(/^ 
x (7 A ^A^) 7p ( 7 ' J 9 (7 A ra )A ti , ; ] , 



(2.1.21) 



or, written in terms of the expressions we defined as H in (I1.3.48P 

tBB 

V HV (2 b) 



n Jf(2 b)( x ~ V) = ^Tr[-f^ a Ylplulxl P la}H[d a , d a ; dp, d x 



(2.1.22) 



If we use the identity for 7 matrices 7 M 7 I /7 P 7 -7 Ai = —1^ a ^ p ^ v [47] and integrate by parts 
the derivatives acting over A xu and A yv , we find that this diagram can be put as 

nj? (2 b \x - y) = 2e 4 Tr[ llllalxlul p la ] (-^H[l, d a ; d% 1] - d* a H[l, d a ; d p d x , 1] 

+d x x H[l,d a d a ; d/3,1] +H[l,d a d a ; d d x ,l}) . (2.1.23) 

Using the properties of the trace, the Clifford algebra definition {7^ , 7^} = 26^ and 
usual 7 matrices results as Tr [7 M 7 A 72,7/3] = 4(5^ 5^ — 5^6^ + 5^5^), we can write the 
right hand-side of the previous equation as 

e 4 [ -8V°#[1, <9 A ; 9a, 1] + 160£#[1, 0„ ; 1] - 8<^#[1, <9 A ; □, 1] 

-16S£tf[l, <9 A ; <9 A <9„, 1] + 16^tf[l, <9 A ; 1] - 16^F[1, «9 M ; d x d v , 1] 

-16^[1, □ ; d v , 1] + 8<^tf[l, □ ; d x , 1] + 16d*#[l, ; d„, 1] 
-16dJ#[l, d„d„ ; <9 A , 1] + 168%H[1, d x d, ; «9 A , 1] - 16H [1, □ ; c^, 1] 
+8^#[1, □;□,!] + 32tf[l, d^d x ; c\<9 A , 1] - 16F[1, <9,A ; □, 1] ] . 

(2.1.24) 

As we have discussed in section 11.3.21 contributions containing a □ can be reduced by 
means of the propagator equation DA = —5 and expressed in terms of J 1 as given in section 
11.3.11 The remaining contributions can be found in the list of renormalized expressions with 
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overlapping divergences given in section [i.3.2l (or can be easily expressed in terms of integrals 
of the list). Hence, we obtain the renormalized result as 



n^f (2 6) ( x ) I R 



12(4vr 
+ ... 



2^3 



In 2 x 2 M 2 + ^\nx 2 M 2 



ar 



\n 2 x 2 M 2 + ^lnx 2 M 2 



(2.1.25) 



Final expression 

Adding the two previous results, the final two-loop renormalized expression is 



up to local terms. 



R 



2 n^f ( 2 a ) (x 



\R 



+ n 



BB 

ixv (2 6) 



X 



4(47T 



R 

e 4 , Q Q r . lnx 2 M 2 

{d„d v -8„ v U)U — + ... 



2\3 



.1- 



(2.1.26) 



2.1.4 RG equation 

When renormalizing the two-loop diagrams we have restricted ourselves to the Feynman 
gauge. We are allowed to do that as the term in the RG equation that takes into account the 
running of the gauge parameter (jad/da) will be shown not to be relevant when verifying 
the two-loop RG equations. To prove this, in appendix O we evaluate the one-loop RG 
equation for quantum gauge fields. We find there the expansion of 7 a to be 

^ = "3(iV + -" (2 ' L27 » 

Along with this, notice that the tree level background effective action and the one-loop 
correction do not depend on the gauge parameter (in this theory we have no quantum- 
background coupling). Hence, as the first gauge corrections arise at the two- loop level, we 
do not have to take them into account in order to verify the two-loop RG equations {pf a d/da 
acting on them is two orders higher in e). So, we are allowed to perform our calculation in 
the Feynman gauge. 

We define the background field two-point function to be 

r£T(s-y) = {d,d v -5, v U)8^{x-v)-^{x- V ). (2.1.28) 

As is detailed in appendix [HI in the background field method the charge and background 
field renormalizations are related: Z e ^Zs = 1. Hence, we will redefine the background field 
to be Bp = \B' as for this new field the anomalous dimension j^' cancels^- Hence, the 
background two-point function up to two loops is found to be 



(dpd u - 5p V U) 



4(4tt 



2\3 



-□- 



lnx 2 M 2 



X 



+ 



9(4vr 2 ) 



5 (4) (x) - 



3(4tt 2 ^ 



-Cl- 



in x 2 M 2 



(2.1.29) 



1 With our definition B'^ 



eoBpo = Zey/ZseBp = B'^ and then j B , 
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that verifies the following RG equation 

( M s + «s) r "= c - < 2130 > 

where, again, (3(e) is the QED /3-function and we have dropped out the term corresponding 
to the running of the gauge parameter, as it is of higher order in e. Then, we obtain the 
following two-loop expansion for /3(e) 

= W^ e3 + M^f e5 + oie7) - (2 ' 01 > 

These results agree with previous ones found in the literature [HI E31 EH ESI ES] . 



2.1.5 Comparison with DiffR 

To stress the key points of our calculation, let us compare this procedure with usual DiffR 
|14j . With Ms and My the one-loop renormalization scales of the fermion self-energy and 
the three point vertex respectively, the Ward identity 



d 

dz» 



Vn(x — z,y — z) = —ie [5^ (z — x) — 5^ (z — y)] S(x — y) 



(2.1.32) 



imposes that these scales are related as 



h.4 = i 

Ml 2 



(2.1.33) 



When dealing with two-loop contributions, in each case one has to make use of the 
corresponding one-loop scale (Ms or My). After a lengthy calculation we find the final 
values for Il£f (2a) and n£f (26) to be 



(2a) ( X ) | R 



n^f (2b)( x )\ R 



24(4tt 2 ) 3 

„4 



(d^-d^Dp 



12(4vr 2 ) 3 

\nx 2 M 2 



{d^-S^np 



In 2 x 2 M 2 + |lnx 2 M 2N 
;r 2 

ln 2 x 2 M 2 + f lnx 2 M 2 ' 

,r 2 



lnx 2 M 2 



x A 



(2.1.34) 



and to obtain the entire two-loop vacuum polarization, we have to use the mass relation 
(12.1.331) to put one of the scales in terms of the other 



In 2 x 2 Ml = In 2 x 2 M v + In x 2 M v + 



1 



(2.1.35) 
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2.2 SuperQED 

In this section we will deal with the supersymmetric extension of QED, SuperQED. As the 
gauge group is abelian, this is one of the simplest examples of supersymmetric gauge theory 
we can consider. This theory was yet renormalized using standard DiffR in [37], where as 
usual, explicit evaluation of Ward identities played a central role. We will re-obtain those 
results applying our procedure. 

2.2.1 Supersymmetry 

Coleman and Mandula [5J showed that the commutators of the generators of any internal 
bosonic symmetry group and the generators of the Poincare group vanish. Hence, space-time 
and internal symmetry groups can not be mixed in a non-trivial way. However, this no-go 
theorem can be avoided by allowing fermionic symmetry generators [B] , and the algebra that 
we obtain is the so-called supersymmetry algebra. Therefore, supersymmetric transforma- 
tions are generated by traslationally invariant quantum operators which change fermionic 
states into bosonic states and vice versa. Hence, for each particle we have another one with 
the same mass and opposite statistic (of course, as this is not observed in nature we conclude 
that if supersymmetry is a fundamental symmetry of nature it is necessarily broken). Other 
relevant consequences of supersymmetry are the positivity of the energy [38l EU EQl EE] 
and, due to the relations imposed to the coupling constants and the equality of the bosonic 
and fermionic degrees of freedom, some cancellations that occur between different Feynman 
diagrams that make supersymmetric theories to be more convergent [3H1 E01 EE] • 

In order to work efficiently with supersymmetric theories, an extension of the usual 
space-time with additional anticommuting coordinates was developed [72] • With this space 
(called superspace) and the extended fields defined in it (superfields), we can perform all 
the calculations with supersymmetry being manifest. In particular, perturbation theory can 
be extended to superspace, which allow us to simplify the calculations as component graphs 
related by supersymmetry are automatically cancelled out. 

In appendix |A] we give a brief introduction to supersymmetry, superspace and the con- 
ventions that we use. Also, we give there a list of references where the reader can find a 
complete treatment of the subject. 

2.2.2 The SuperQED model 

A supersymmetric abelian gauge theory can be defined in terms of a field strength W a [38] 
which is a chiral superfield (D a W a = 0) that verifies 

D a W a = -D"W a . (2.2.36) 

Hence, this field can be expressed in terms of an unconstrained real scalar superfield V — V 

as 



W a = iD 2 D a V 
W a = -iD 2 D a V 



(2.2.37) 
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From the algebra of covariant derivatives flA.3.9j) we can easily conclude that W a is invariant 
under the transformations 

y> = V + i(A-A), (2.2.38) 

where A is a chiral parameter. 

The relevant action that we find is 

S = J d 4 xd 4 9 W 2 = ^J d 4 xd 4 9 VD a D 2 D a V , (2.2.39) 

which is the supersymmetric gauge invariant generalization of the action for a free vector 
field, as can be seen if we write this expression in component fields [38]. As the matter part 
of the action can be expressed in terms of a chiral field $ as J d s z $e y $ [38| [7TI [69] , the 
supersymmetric extension of massless QED is an action of the form [71] 

S = J d 4 xd 2 9 W 2 + J d 4 xd 4 9 <5 + e 9V $ + + J d 4 xd 4 9 §^e~ 9V §^ . (2.2.40) 
Background field method 

In section IB. 21 of appendix [B] we discuss in detail the application of the background field 
method to supersymmetric gauge theories. It is found there that when dealing with an 
abelian theory, we have a linear quantum-background splitting of the form 

V^V + B, (2.2.41) 

where V and B are the quantum and background gauge fields respectively. The background 
effective action that we have is then of the form 

T[B] = J d 4 xd 4 yd 4 9 B(x, 9) [D a D 2 D a B(y, 9)] T{x-y) + ... (2.2.42) 

Notice that the background field method allows us to obtain the beta function (3(g) only 
from the renormalization of the background field two-point function. Hence, we have to 
obtain the the one- and two-loop coefficients of the expansion of T(x), as with them we can 
obtain (5(g) up to two- loop order. 



2.2.3 One- loop level 

Although we can also consider a diagram which corresponds to a tadpole contribution of a 
loop of <&± fields, this diagram gives no contribution as CDR imposes ADA|_r = 0. Hence, 
the only relevant one-loop contribution from Q± fields to the background vacuum self-energy 
is the diagram shown in 12.41 We will obtain the contribution corresponding to a $ + loop, 
which is denoted cLS r^|_ ^ , clS the other one is exactly the same. With the superspace 
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Figure 2.4: One-loop SuperQED diagrams. Thick lines correspond to external background 
fields and solid lines represent $ + or <£>_ propagators. 



Feynman rules defined in appendix |Al we find the expression of the diagram to be (notice 
that the superspace propagator is P^- = A(xj — Xj)5 4 {9i — 9j) = A XiXj Sij) 



|- / ^rf 8 z 2 C(:,)C(: 2 ! 



^1^12 ^2 



P 2 P 12 D\ 



9_ 
2 



d 8 Zl d 8 z 2 B( Zl )B(z 2 ) [D 2 2 D 2 2 P 12 ] [D 2 2 D 2 2 P 21 ] 



(2.2.43) 



Applying D-algebra we remove all the derivatives from the first propagator and make them 
act over the external fields and the other propagator. So 



£ J d 8 Zl d 8 z 2 B( Zl ) [D 2 D 2 B( Z2 )} P 12 [D\D\P X2 \ 
+^ J d 8 Zl d 8 Z2 B( Zl )B( Z2 )P 12 [uD 2 2 D 2 2 P l2 ] 

-t- J d 8 Zl d 8 Z2 B( Zl ) [D«D a B( Z2 )] P 12 [d 2 aA D 2 2 D 2 2 P 12 ] , (2.2.44) 

where we have used the results of the D-algebra D 2 D 2 D 2 = DD 2 and [D a , P 2 ] = — 2<9 Q(i P Q . 
At this point, we can apply the identity (IA.5.541) for supercovariant derivatives 5i 2 D 2 D 2 5\ 2 = 
5i 2 . This gives us one free 6*-space ^-function that we can use to evaluate one of the 9 integrals. 
With this, we have (using the identifications x\ = x and x 2 = y) 



,(i) 



, 2 



9 — J d 4 xd 4 yd 4 9 B(x,9) [D 2 D 2 B(y,9)] A 

2 f 

+ L I d 4 xd 4 yd 4 9 B(x,9)B(y,9)A xy (nA xy ) 

d 4 xd 4 yd 4 9 B(x,6) [D«D a B(y,9)} A xy d y adi A 



' 2 



(2.2.45) 



This is the bare expression that we have to renormalize applying CDR rules. It is clear 
that the second term does not contribute as CDR imposes AdA|# = 0. The CDR renormal- 
ization of the third term Ad a aA\R = \d a aA 2 R allows us to integrate by parts the space-time 
derivative and make it act over the external background fields. Hence, using the identity 



D 2 D 2 + -d a6l D a D c 



-D a D 2 D r 



(2.2.46) 
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we find the final renormalized expression to be 

,2 



^(1) 



R 



<£_ 

4 



dtxdtydte B(x, 9) [D a D 2 D a B(y, 9)] [A 2 ] R {x-y) 



16(4tt 



2\2 



d 4 xd A yd A 6 B(x,G) [D a D 2 D a B(y,9)] □ 



ln(z - y) 2 M 2 
(x - yf 



(2.2.47) 



As we can see, this term is manifestly gauge invariant, as is guaranteed by the use of CDR. 
The total one-loop effective action is the sum of both contributions corresponding to $+ and 



R 



F (l) , r (l) 
+ R ' — R 



8(4tt 



2^2 



d A xd A yd A 9 B{x,9) [D a D 2 D a B(y, 9)] □ 



ln(x - y) 2 M 2 
(x - yf 



. (2.2.48) 



2.2.4 Two-loop level 







Figure 2.5: Two-loop SuperQED diagrams. Thick wavy lines correspond to external back- 
ground gauge fields, thin wavy lines correspond to quantum gauge field propagators and 
solid lines correspond to $ + or $_ propagators. 



As we can see from figure 12.51 two-loop calculations involve the quantum gauge field 
propagator. Although this propagator depends on the gauge parameter, we can use it eval- 
uated in Feynman gauge, as the term that takes care of the running of the gauge parameter 
in the RG equation (j a (g)d/da), starts its expansion in terms of the coupling constant at 
order g 2 . Hence, 7^ acting over any two-loop diagram will be two orders higher in g, and 
consequently, not relevant when verifying the RG equation. We will detail this later. Also 
notice that, as in the one-loop case, we will obtain only contributions corresponding to $ + 
fields, as the diagrams with $_ fields have the same expression. 
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Before obtaining each diagram, let us discuss an useful identity we will use when we have 
to apply -D-algebra to a supergraph j37j. Let F{ be a function of superspace coordinates 
Zi = (xi, 9i), and suppose that we have an expression where the only dependence in Zi is of 
the form 

J d 8 Zl d 8 Zj d 8 z k ...Fi [B 2 b 2 P^ [D 2 D 2 P ik ] . . . (2.2.49) 

Applying integration by parts rules in superspace (which can be found in section IA.3.21 of 
appendix |A|) and D-algebra, we obtain 

[D 2 F % ] [D 2 P tj ] [D 2 D 2 P ik ] . . . 
Ft [Dfrj] [UD 2 P ih ] . . . 

[DZF,] [DfPij] [d\ a D«D 2 P lk \ . . . (2.2.50) 

If we integrate by parts again to remove all the superspace derivatives from P^-, it is clear 
that we will obtain some contributions that will vanish by the identities (IA.5.541) . when we 
use the #-space 5-functions to set j = k. As an example, consider [-D^Fj] Py [ODfP ik ], which 
is obtained from the second term of (12.2.501) . As can be seen, this contribution will cancel 
when we set j = k and apply 5ijD 2 5ij = 0. So, the final relevant terms of the expansion of 
(12.2.491) are 

j d 8 Zi d 8 Zj d 8 z k . . . [D 2 D 2 F t ] P l3 [D 2 D 2 P lk ] . . . 

+ j d 8 z i d 8 z j d 8 z k . . . FP l3 [uD 2 D 2 P lk \ . . . 

- i J d 8 z i d 8 z j d 8 z k . . . [DfDfFi] P y [d^D 2 D 2 P lk ] . . . 

+ (terms that vanish when j=k) . (2.2.51) 



+ / d 8 Zid 8 Zjd 8 z k . 



I d 8 z i d 8 z j d 8 z k 



Diagram (a) 

The bare expression of this diagram is 

rf } = ~y / d 8 Zl d 8 z 2 d 8 z 3 d 8 z 4 B( Zl )B(z 2 )P 43 [D 2 3 DlP A3 ] [D 2 2 D 2 2 P A2 ] 
x [DlDlP Vi ] [D 2 2 D 2 2 P 21 ] 

4 f 

- 9 - J d 8 Zl d 8 z 2 d 8 z 3 d 8 z 4 B( Zl )B(z 2 )P 43 [D 2 4 D 2 4 P A3 ] [D 2 2 D 2 2 P 2A ] 

x [D?D?P 31 ] [DlD\P l2 ] . (2.2.52) 

We will study the first contribution, as the second one, except for having D and D 
interchanged, is the same. Using the identity P i3 D 2 D 2 P i3 = A| 3 5 43 , integrating by parts 
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and applying P-algebra we find 

rf 7) = -y / d 8 Zl d 8 z 2 d 8 z 3 d 8 z 4 B( Zl )B(z 2 ) (A%5 43 ) [D 2 2 D 2 2 P 42 ] [D 2 D 2 P 13 ] [D 2 2 D 2 2 P 21 ] 

= - y J d 8 z x d 8 z 2 d 8 z 3 d 8 z 4 B( Zl )B(z 2 ) (A 2 3 5 43 ) [D 2 2 D 2 2 P 42 ] (nP 13 ) [D 2 2 D 2 2 P 21 ] . 

(2.2.53) 

Making use of (12.2.511) we can write this a 

r? /} = -y / ^i4<l 8 A%) [£> 2 £> 2 5(z 2 )] (A 2 3 5 43 ) [£ 2 2 £ 2 2 P 42 ] (aPra)^ 

_|. J d 8 z l d 8 z 2 d 8 z 3 d 8 z A B{z 1 )B{z 2 ) (A 2 3 5 43 ) [nP 2 2 P 2 2 P 42 ] (DPia) P21 

+^ y d 8 z 1 d 8 z 2 d 8 z 3 d 8 z 4 B(z 1 ) [D«D a B(z 2 )] (A 2 3 5 43 ) [^D'D^] (DP 13 ) P 21 . 

(2.2.54) 

After using identities (IA.5.541) we can evaluate three of the 9 integrals with the free 5- 
functions. Then, with the identifications x\ = x, x 2 = y, x 3 = u and x 4 = v, the contribution 
becomes 

lf a/) = - 9 - j d A xd A yd A 6 B(x,8) [D 2 D 2 B(y,6)]A xy J A yv (DA xu )A 

- 9 1 J d A xd A yd A 6 B(x,6)B(y,9)A xy J d A ud A v (a A yv ) (OA xu ) A 

+!|_ J d A xd A yd A 6B(x,9) [D«D a B(y,9)] A xy J d A ud A v (c^A^nA^A 

(2.2.55) 

Remembering DA = —5 and the definition of integral expression J 1 , this can be put as 

r (2aJ) 9 



, 2 

uv 



2 

uv 



d A xd A yd A 6 B(x, 6) [D 2 D 2 B(y, 6)] [AI 1 ] (x - y) 
+ t J d A xd A yd A 6 B(x } 6)B(y,6)Al y 

+ !|! J d A xd A yd A 6 B(x, 6) [D dt D a B(y, 9)} [Ad^I 1 ] (x - y) . (2.2.56) 

The second contribution of (12.2.521) only differs from the first one by the interchange of D 
and D. Hence, using P-algebra identities f)A.3.9l) . the total bare expression of diagram (a) 
is found to be 

p (2a) _ r (2a/) v (2aII) 
1+ — 1 + +1 + 

= 9 — j d A xd A yd A 9 B(x, 6)B(y, 6) [□(A/ 1 ) - 2A 3 - d a& {Ad a6l I 1 )] (x - y) 

4 r 

+ 9 — / d A xd A yd A 6 B(x, 6) [D a D 2 D a B(y, &)] [AI 1 ] (x - y) . (2.2.57) 
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Diagram (b) 

The bare expression of this contribution is 

rf> = -g 4 j d 8 Zl d 8 z 2 d 8 z 3 B( Zl )B(z 2 )P 31 [D 2 3 D 2 3 P 31 ] [D 2 2 D 2 2 P 23 ] [D 2 x DIP 12 } 

-g 4 J d 8 z 1 d 8 z 2 d 8 z 3 B(z 1 )B(z 2 )P 31 [D 2 3 D 2 3 P 23 ] [D?D?P 31 ] [D^|P 21 ] . 

(2.2.58) 

As in the previous diagram, the two terms that form (12.2.581) differ only by the interchange 
of D and D, so we will obtain the first one, that we name T+ bI \ With identities (I2.2.5ip 
and (1A.5.54I) we find the relevant expansion of rf° to be 

rf /} = -g 4 J d 8 Zl d 8 z 2 d 8 z 3 B( Zl )B(z 2 ) (A 2 31 5 31 ) [D 2 2 D 2 2 P 32 ] [D 2 2 D 2 2 P 12 ] 
= -g 4 J d 8 z l d 8 z 2 d 8 z 3 B{z x ) [D 2 2 D 2 2 B(z 2 )] {A 2 31 5 31 ) [D 2 2 D 2 2 P 32 ] P 12 
-g 4 J d 8 z l d 8 z 2 d 8 z 3 B( Zl )B(z 2 ) (A 2 31 5 31 ) [UD 2 2 D 2 2 P 32 ] P l2 
+ig 4 [ d 8 z x d 8 z 2 d 8 z 3 B( Zl ) [D a D«B(z 2 )} (A 2 31 5 31 ) [d 2 aA D 2 2 D 2 2 P 32 ] P, 



12 • 

(2.2.59) 

Using identities (I A. 5. 541) we can get rid of the superspace derivatives and obtain free grass- 
manian 5-functions. After evaluating the grassmanian integrals, the identifications X\ = x, 
x 2 = y and x 3 = u allow us to obtain 



r (26/) = _ g 4 J d i xd ± yd ± e B ^ [£> 2 D 2 B(y 1 9)] A xy J d 4 u A yu A 
-g 4 [ d 4 xd 4 yd 4 9B(x,9)B(y,9)A xy [ d 4 u (UA yu )A 2 xl 



2 

xu 



, 2 

xu i 



+ig 4 j d 4 xd 4 yd 4 9 B(x,9) [D a D 6l B(y 1 6)] A xy j d 4 u (d y a< ,A yu )A 2 a 

(2.2.60) 

or, in terms of the J 1 integral expression 

T (2bi) = _ g4 J d A xd 4 yd 4 B ^ ^ \p^ 2 B(y, 9)} [A/ 1 ] (x - y) 
+g 4 J d 4 xd 4 yd 4 9 B(x, 9)B(y, 9)A% 

-ig 4 J d 4 xd 4 yd 4 9 B(x, 9) [D a D«B(y, 9)] [AO^I 1 ] (x - y) . (2.2.61) 
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Finally, adding up the other contribution and using D-algebra relations flA.3.91) we find the 
total bare contribution to to be 



,(26) 



//rfW* B (xM) [-□(AT) + 2A3 + fl- (Ada^I 1 )] (x - „) 
-# 4 / d 4 xd 4 yd 4 6 B(x, 0) [D a D 2 D a B(y, 9)] [A/ 1 ] (x - y) . (2.2.62) 

Diagram (c) 

This diagram is given by 

lt c) = -y / d s z 1 d 8 z 2 d s z 3 d 8 z i B( Zl )B(z 2 )P u [DlDlP lA ] [5?D?P U ] 

X [Z^Pjb] [^2^2^24] • (2.2.63) 

Applying identity fl 2 . 2 . 5 1 [) we can split this expression in three contributions as 

rf) = r (^) +r jc/i) + r (2c//i )) (2 2 64) 



with 








r (2c/) 


9 2 


1 + - 


2 




x[ 






r (2cI7) 

1 + ~~ 


9 4 


2 






v {2cIII) 




1 + — 


2 




x [ 



d 8 z l d 8 z 2 d 8 z 3 d 8 z 4 B{ Zl ) [D 2 D 2 B(z 2 )} P 23 P U [D 2 2 DlP 24 \ [D 2 D 2 P 14 ] 



d 8 z l d 8 z 2 d 8 z 3 d 8 z 4 B(zi) [D«D a B(z 2 )} P 23 P 34 [d 2 aA D 2 2 D 2 2 P 24 ] 
[D 2 D 2 P U ] [DfflPu] . (2.2.65) 

We will evaluate each contribution separately. Starting by r+ c/ ' ) we can apply again fl 2 . 2 . 5 1 1) 
and write this as 

lt c7) = -y / d 8 Zl d 8 z 2 d 8 z 3 d 8 z 4 [D 2 D 2 B( Zl )] [D 2 D 2 B(z 2 )] P 34 P 23 P 24 P 14 [D 2 D 2 P 13 ] 
_|. y d 8 z l d 8 z 2 d 8 z 3 d 4 z 4 B{z x ) [D 2 D 2 B(z 2 )} P 34 P 23 P 24 P l4 [udId{P 13 ] 

■ 4 /■ 

+ l JL y d 8 z l d 8 z 2 d 8 z 3 d 8 z 4 [D«D a B( Zl )} [D 2 D 2 B(z 2 )} P 34 P 23 P 24 P l4 

(2.2.66) 

Integrating by parts and using the anticommutative nature of the superspace derivatives, we 
find that an expression of the form J dxdyd6[D a A(x , 9)][D 2 B(y, 9)]f(x — y) vanishes. Hence, 
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the first and third expressions automatically cancel. Using the identifications x\ = x, x 2 = y, 
x 3 = u and X4 = v we have for this contribution 



, 2 

XV 



rJ c/) = ^ J d 4 xd 4 yd 4 9 B(x,9) [D 2 D 2 B(y,6)] A xy J d 4 v A yv Al 

J d 4 xd 4 yd 4 9 B(x, 9) [D 2 D 2 B(y, 9)] [A/ 1 ] (x - y) . (2.2.67) 



9 4 



We continue now evaluating T+ cII \ As in this case we have the product of the superpropa- 
gators P23-P34) we can use one °f the free grassmanian 5-functions and evaluate the integral 
over 9 3 . With this, we can write this expression as 

Tf U) = - 9 - J d 8 z 1 d 8 z 2 d 4 x 3 d 8 z i B(z 1 )B(z 2 )A 23 A u (nA 24 )5 2i [D 2 D 2 1 P u ] [D 2 DlA 13 6 12 ] 

= 9 - J d 8 Zl d 6 z 2 d 4 x 3 B{z x )B{z 2 )A 2 3 [D 2 D 2 P 12 ] [D 2 D 2 A 13 S 12 ] . (2.2.68) 

We have performed the integral over z 4 applying OA(x 2 — x&) = —6(x 2 — £4) and the 
grassmanian 5-function 5 2 4- Introducing again the grassmanian coordinate 9 3 with a 5- 
function and integrating by parts the superspace derivatives that act over A 13 5i 3 (= P 13 ) we 
find 

rf JJ) = y / d 8 Zl d 8 z 2 d 8 z 3 B{ Zl )B{z 2 ) [D 2 2 D 2 2 A 2 23 6 23 ] [D 2 2 D 2 2 P 12 ] P 13 . (2.2.69) 

At this point, applying (12.2.511) . we can integrate by parts the superspace derivatives acting 
over P\ 2 and arrive to 



t H) = 7T / d 8 z x d 8 z 2 d 8 z 3 B{ Zl ) [D 2 D 2 B(z 2 )} [D 2 2 D 2 A 2 3 6 23 ] P 12 P, 



13 



4 f 

+ 9 - J d 8 z l d 8 z 2 d 8 z 3 B( Zl )B(z 2 ) [UD 2 D 2 A 2 3 5 23 ] P 12 P 13 

- x g l J d 8 z 1 d 8 z 2 d 8 z 3 B( Zl ) [D a D«B(z 2 )] [d 2 a „D 2 2 D 2 2 A 2 23 5 23 ] P 12 P 13 . 



(2.2.70) 



These expressions can be evaluated straightforwardly with the identities ( 1 A. 5. 541) . Using also 
the coordinate identifications x\ = x, x 2 = y, x 3 = u and X4 = v, we have this contribution 
written in terms of the integral expression J 1 as 



,(2di) _ g 



d 4 xd 4 yd 4 9 B(x, 9) [D 2 D 2 B(y, 9)] [AI 1 ] (x - y) 
-t. J d 4 xd 4 yd e B(x,9)B(y,9)Al y 

+i t J d 4 xd 4 yd 4 9 B(x, 9) [D a D«B(y, 9)] [Ad^I 1 ] (x - y) . (2.2.71) 
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Finally, we take care of y+ cIII \ With (12.2.511) . taking into account that the superspace 
derivatives anticommute and making the usual identifications X\ = x, x 2 = y, x 3 = u and 
x 4 = v, we find 

Vf UI) = i 9 — j ' d A xd A yd 4 6 B(x,6)[D«D a B(y,6)] 
x/«,A„4, Ki A 1 ,)A.,pA„) 

4 i> 

+ L / d 4 xd 4 yd A 9 \D $ D p B{x,6)j [D & D a B(y,e)} 

x J d\d 4 v A yu A uv {d y ^A yv )A xv {^A xu ) , (2.2.72) 

or which is the same, integrating by parts the superspace derivatives of the last integrall and 
using the integral expression H defined in (11.3.481) 

Tf IH) = / d 4 xd 4 yd 4 9 B(x,6) [D A D a B(y,6)] [Ad a6l I l ] (x - y) 

+ L J fatfytfQ B (x, 6) [D l5 D 2 D a B{y 1 6)] C^H[d^, 1 ; 1, d a& ] . 



[2.2.73) 



Adding up the three contributions, we find the final bare expression of rj c) to be 
lt c) = Y J d * xdA y d4e B ( x ' 6 ) [D a D 2 D a B(y,6)] [A/ 1 ] (x - y) 

+|. j d A xd 4 yd 4 6 B(x, 9)B(y, 9) [□(A/ 1 ) - A 3 - d a& (Ad^I 1 )] (x - y) 



+ >L J d^xd%d'BB{x,B) [D 13 D 2 D a B (y , 9)] C p& H[d pfi , 1 ; 1, d a& ) , 

(2.2.74) 

with C@ a given in section IA.1I of appendix [Al 
Diagram (d) 

The bare contribution of this diagram is 

rf ) = -y / d 8 Zxd 8 z 2 B( Zl )B(z 2 ) [DlDlP 12 ] P 12 [D 2 2 D 2 2 P l2 ] . (2.2.75) 

Applying the identity 5\ 2 D 2 D 2 5i 2 = 5\ 2 it is clear that, with the identifications X\ — x and 
x 2 = y, the non-renormalized expression of this diagram is 

4 p 

= - 9 — d 4 xd 4 yd A 6 B(x,6)B(y,6)Al y . (2.2.76) 



2 J[D^D^B(z 1 )][D^B(z 2 )]f a&t ^(z 1 -z 2 ) = -J B(z 1 )[Df i Df ! D-D a B(z 2 )}f a . M (z 1 - z 2 ). Also it has 
to be noted that due to the anticommutative nature of the superspace derivatives D^D? = -D 2 C &13 
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Renormalization 

As our renormalization procedure verifies CDR at one loop, in order to obtain the final renor- 
malized result we can replace each bare expression for its renormalized value and simply add 
them. Even more, as with our procedure each expression always has the same renormal- 
ized value, we can first add all the bare expressions, and then perform the renormalization. 
This is forbidden when we use DifTR, as each expression has to be renormalized with its 
corresponding scale. From the explicit form of the different contributions, it is clear that all 
the terms cancel exactly, except the last part of diagram (c) (12.2.741) . Hence, the two-loop 
renormalized contribution to the vacuum self-energy is (multiplying by two as we consider 
both contributions from the chiral matter superfields $ + and $_) 



= 2 (rf ) + rf ) + rf ) + rf ) 



R 



g 4 J d 4 xd 4 yd 4 9 B(x, 9) [D /3 D 2 D a B(y, 9)] C^H R [d^, 1 ; 1, d a& ] 

d 4 xd 4 yd 4 9 B(x, 9) [D a D 2 D a B(y, 9)] □ ^ V) ™ +... 



16(4tt 2 ) 3 J * " v *' /J (x- y y 

(2.2.77) 

where, when obtaining the final result, we have directly applied the corresponding identity 
from the list of integrals with overlapping divergences of section 11.3.21 Let us remark again 
that, as was guaranteed by fulfilling CDR rules at one loop, this expression is directly gauge 
invariant. 

Before dealing with the background RG equation, we have to justify the use of the 
Feynman gauge in the calculations. As for the QED case, in appendix [UJ by means of 
the evaluation of the one-loop RG equation for the quantum gauge fields, we obtain the 
expansion of the term that takes into account the running of the gauge parameter in the RG 
equation {pf a d/da). This is of the form 

7. = ~(£^ 2 + ... (2.2.78) 

At this point, notice that the first gauge corrections to the background effective action arise 
at the two-loop level. Hence, when verifying the two-loop RG equation, we may not take 
into account ^y a d/ da acting on them, as this is two orders higher in g. This is the reason why 
we are allowed to use the Feynman gauge in our calculations in both QED and SuperQED 
models. 

Let us now proceed to the evaluation of the RG equation for the background gauge 
field self-energy. As in the QED case, if we make the redefinition B — > ^B, the anomalous 
dimension term of the renormalization group equation cancels (remember that the coupling 
constant and the background field renormalizations are related: Z g \fZ~B = 1 )• So, with this 
definition, the background field two-point function to two-loop order is 

1 s 1 lnx 2 M 2 q 2 \nx 2 M 2 

r,(x) = -6* {x) -—0— -^a—^ + ..., (2.2.79) 
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and it fulfills the following RG equation 



T R (x) 



(2.2.80) 



where we do not consider the term that takes care of the running of the gauge parameter. 
By solving this equation order by order in g it is clear that the beta function is of the form 



Pi 
ft 



/3 1 g 3 + f3 2 g 5 + 0(g 1 
1 



167T 2 
1 



5(4tt 



2\2 



(2.2.81) 



However, as with our supersymmetry conventions the gauge coupling constant differs from 
the usual one gsQED (g = V^gsQED) [38j , the expansion of the beta function to two-loop 
order in terms of the usual coupling constant is 



P(9sqed) - 2U^t\ 9 sqed + 2 ^ 7C 2y g SQED + OWsqed) 



(2.2.82) 



This agrees with previous results found in the literature [39j HOJ SI]- 





(d) 



Figure 2.6: One-loop SQED Ward identities. 



Let us compare this procedure with the steps that we have to take when using usual DiffR. 
We have to consider first the Ward identities, that can be shown to relate the 3-point 1PI 
Green's function < T £>$ + <l> + >ipj and the 2-point 1PI Green's function < T $ + <l> + >ipj 
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as [37] 



D 2 ( Zl ) < T B( Zl )^{z 2 )^{z 3 ) > 1PI = -g<T $_(^)$_(^ 3 ) >ipi D 2 ( Zl )S 8 ( Zl - z*) 
D\ Zl ) < T B( Zl )^( Z2 )^(z 3 ) > 1PI = -g<T $_(z 2 )$_(^) > 1PI D 2 ( Zl )5 8 ( Zl - z 3 ) 

D\ Zl ) < T B( Zl )$ + (z 2 )$ + (z 3 ) > 1PI = g<T <Z> + ( Z2 )$ + ( Zl ) >ipi D 2 ( Zl )6 8 ( Zl - z 3 ) . 



With these identities, we can obtain the one-loop relation between the scales that renor- 
malize the B — $ — $ vertex functions (diagrams (a)-(c) of figure with scales My a , My b 
and My c respectively) and the self-energy corrections (diagram (d) of figure 12.61 with 
scale My s ). Performing the explicit renormalization and imposing identities (I2.2.83p . this 
relation is found to be [37] 



Hence, when renormalizing each of the two-loop diagrams, we have to use the corre- 
sponding one- loop scale, add up all the results and apply (12.2.84[) . In [37] is shown that 
this relation cancels contributions that came from different diagrams and are grouped in an 
expression multiplied by hi[(My My ) /(My MyJ]. As can be seen from our procedure, these 
cancellations take place automatically once we have renormalized the one-loop divergences 
with the rules of CDR. 



(2.2.83) 



M 2 M Vs = 



Ml M 2 r 



(2.2.84) 
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Chapter 3 

Non-abelian QFT applications 



3.1 Yang-Mills 

3.1.1 Conventions and definitions 

Relevant group theory definitions 

Let G be a continuous symmetry group with generators T a . We can define an associated Lie 
algebra through the commutation relation 

[T a ,T b ]=if abc T c , (3.1.1) 

where f abc are the structure constants of the Lie algebra, which obey the Jacobi identity 
jadejbcd _|_ jbdejcad _|_ jcde jabd _ q ^y e can ^ severa ] representations of this Lie algebra 

in terms of matrices t": one of them is the adjoint representation, denoted by r = G, 
where the representation matrices are given by the structure constants {t b G ) ac = if abc . These 
representation matrices are found to satisfy 

tr[t a r t b ] = C(r)S ab 

Y,tX = C 2 (r)l, (3.1.2) 

a 

where C(r) and C*2(r) are constants, being the latter called the quadratic Casimir operator. 
For the concrete case of the adjoint representation, we write the relation for the Casimir 
operator as 

^acdjbcd = CAS ab _ (3.1.3) 

where we define Ca = 02(G). 
Yang- Mills model 

Yang-Mills theory is one of the simplest examples of non-abelian gauge theory [32]. It is 
obtained by imposing invariance under a continuous symmetry group. We start by consid- 
ering V(x) to be an unitary n x n matrix representing one of the elements of a gauge group. 



39 



40 



CHAPTER 3. NON-ABELIAN QFT APPLICATIONS 



Then, the fields transform according to 

il>(x) — > V(x)ijj(x) 

= (l + iw a (x)t a + 0(w 2 ))tfj(x) 



(3.1.4) 



where we have considered an infinitesimal parameter w a , which has allowed us to expand V(x) 
in terms of the generators of the symmetry group. Now we have to construct a covariant 
derivative that, when acting over ip(x), has the same transformation as the field. This 
derivative, expressed in terms of a connection (gauge potential) A^^ 

□3 



A°;i%, is found to be 



D 



WAX 



(3.1.5) 

If we choose the adjoint representation, this becomes D ab = <9 M 5 afe + gf abc A b . As this 
derivative has to transform covariantly under the gauge group, the infinitesimal gauge trans- 
formation of A® is found to be [H] 



4" 



1 



A a a + -d„w a - f abc w b A c u + 0(w 2 ) 
9 

A; + - g (D,w) a + 0(w 2 ). 



(3.1.6) 



By considering the commutator of covariant derivatives, we can define a field strength as 
—igF® u t a = [Dp , D v ] that in terms of the gauge potential has the form 

F; v = 8^1-8^ + gj ahc A\Al. (3.1.7) 

With this field strength is straightforward to define an gauge invariant quantity that is 
the Yang-Mills action: 



S 



(3.1.8) 



As in the abelian case, when quantizing the action in a path integral approach, we have 
to fix the gauge in order to suppress all the equivalent field configurations obtained from 
a given one through gauge transformations. The result is that the gauge-fixed partition 
function Z is 



Z[J]= / [dA] det 



6G a (A v 
5w b 



exp 



J w=0 



-S(A) - — I d*x G a G a + J a u A 
2a 



a 



(3.1.9) 



where G a is the gauge-fixing function. Writing the determinant in terms of anticommuting 
ghost fields^ and choosing for the gauge- fixing function G a = d^A^ we can find the complete 
Yang-Mills lagrangian to be 

c = ^F^ + ^d^rid^r + ^nriD^r. (3.1.10) 

This implies that we have the following gauge field and ghost propagators 

<Al{x)A h M> = 5,J ab A(x-y) 
<ri a (x)ff(y)> = 5 ab A{x-y) (3.1.11) 
1 Being 8 an anticommutative variable J HdOdOe^ 6 ^^ 6 ! = det[a] 
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Background field method 

As is detailed in appendix [Bj with the standard quantum-background splitting A® — > A^+B^ 
we can define two gauge covariant derivatives D^ c = d^5 ac + gf abc B^ and V> a ^ = d^5 ac + 
gf abc (B^ + A b ^). Using them, and with a background covariant gauge-fixing function as 
G a = (jy^A^Y, we find the split lagrangian to be written as 

£ = \f; v f; v + ^-(B,A,y(B u A u y + (n^nv^y , (3.1.12) 

with the field strength depending in both quantum and background fields F£ v = F^ U (A + B). 
As in the previous abelian examples, we will perform the calculations in Feynman gauge; 
however, there is one important difference. As can be seen from the lagrangian (13.1. 12|) . we 
have an interaction term of the form 

gf* [2(d lt B;)A i lt At - B^A^Al] , (3.1.13) 

which implies that the one-loop background self-energy depends on the gauge parameter, as 
we have a loop with quantum gauge fields propagators. Hence, although the term in the RG 
equation that takes care of the running of the gauge parameter, jad/da, will be shown to be 
of order g 2 (like in QED and SQED), in this case we can not leave it out, as when acting on 
the one-loop contribution it will affect the verification of the two-loop RG equation. Then, 
our procedure will be as follows: first of all, the standard gauge fixing parameter a will be 
redefined here as - = 1 + £, so that usual Feynman gauge (a = 1) will correspond to £ = 0. 
We will obtain the one-loop contribution to the background self-energy in this gauge. Then, 
by means of functional methods, we will expand the complete effective action at one loop at 
second order in the background fields and retain only the linear dependence in £, that we 
term 1^ . We apply this procedure as in the renormalization group equation we will first 
take derivatives with respect to this parameter and after that impose Feynman gauge (£=0). 

Hence, we have a background effective action of the form 

r.//[B] = \J d 4 xd 4 yB;(x)T^ ab (x - y)B b u (y) + ... 

= i y d 4 xd*yB;(x) [5 ab (d,d u - 5, u U)5^\x - y) - U B ^ ab (x - y)] B b v {y) + ... 
= S [B] + Tf d'xd'y B;(x)n*r b (x - y)B b v (y) + . . . , (3.1.14) 

where So is the tree-level background two-point function. 

In figure 13. 1[ the interaction vertices derived from the quantum-background split action 
which are relevant to this work are shown. We have the following corresponding Feynman 
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Figure 3.1: Relevant interaction vertices of the Yang-Mills quantum-background split action. 
Thick lines represent external background fields, thin lines are quantum gauge propagators 
and dashed lines correspond to ghost propagators. 



rules (evaluated in Feynman gauge) 



(vl) 
(v2) 

(v3) 
(v4) 
(v5) 
(v6) 



gf 



aba 



+ / J ( VV ~ pp 6, 

gf abc 



A " da") + 25^8, 



PI 

p,v w p 



-9f abc dT 



-gf abc {dl .„ 

2 radx rxbcc 

y J J 



(3.1.15) 
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3.1.2 One-loop level 

Although the background self-energy is all that we need to find the one-loop beta function, 
we will also obtain the linear dependence in the gauge parameter £ of the effective action 
calculated in a generic gauge and expanded to second order in background fields. We need 
this contribution in order to take care of the running of the gauge parameter in the RG 
equation. 



Correction to the propagator 

This is the sum of two different diagrams: one with a loop of quantum gauge fields, and 
another of ghost fields, as can be seen in figure 13.21 In these diagrams we have to apply the 
CDR procedure: First we write the expressions in terms of the basic functions defined in 
(11. 2. 371) . and after that we replace them with their renormalized values. 

Here and in the rest of the diagrams of the Yang- Mills theory, D^' y denotes a space-time 
derivative acting over one external field. Applying the Leibniz rule, D^ ,y becomes a minus 
derivative acting over the propagators. The bare expression for both contributions is 

• Gauge loop 



g2 jacd jbdc 



-A 



-26^D* + 6 pa (&°-&°)+26 tip D* 



x 



-26 vp D* + 5 pa {dl- dl) + 25 ua D\ 



A 



9 2 C A 5 



ab 



[8d p d u A 2 - 8^DA 2 + 8d p {Ad u A) - WAd^A] 



Ghost loop 



(3.1.16) 



- g 2 f abc f bcd A xy (d; -d${ft- dl)A xy = -g 2 C A 6 ab [2^(A9,A) - AAd^A] . 

(3.1.17) 

Adding the two previous results we find the total non-renormalized contribution to be 
: /C^ afe [4^A 2 -4^nA 2 + 2^(A^A)-4A^A] . 



n^f (i) (%) 



(3.1.18) 
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It is worth to mention here again that we are allowed to do the previous step (adding 
up the expression even before renormalizing) because we are using CDR, as we pointed out 
previously. With CDR the basic functions are renormalized always with the same expression, 
despite their origin. In an usual DiffR procedure, we first have to renormalize each diagram 
in a separate way, relate the different scales that appeared via the Ward identities, and only 
after that we can add up the results. 

Replacing the values of CDR for A 2 and Ad^d^A, the renormalized one-loop contribution 
to the B a propagator is obtained as 



R 



g 2 C A 5 a \d,d v -5, u u) 
g 2 C A 5 a \d^ u -5^n) 



11 



A 2 R (x) 



11 



lnx 2 M 2 



487T 2 (47T 2 



72tt 



■Ax) 



(3.1.19) 



As a check, the result found here is automatically transverse, fulfilling the corresponding 
Ward identity. 



Effective action in a generic gauge 

As we have discussed previously, in order to take care of the running of the gauge parameter 
in the RG equation, we will obtain the linear dependence in £ of the one-loop background 
effective action expanded to second order in the background fields. To perform this calcula- 
tion we consider a functional approach: to obtain the exact one-loop effective action it is well 
known that we have to consider only the part of the lagrangian quadratic in the quantum 
fields H61I17]. This part is 



r(2) 

gauge 



gf^B^Al + ~(D M A,r(D M A,r + |(D^r(D^r 
[5^ ab - 2gr b B; u + £(D M D„r fe ] A\ , 

abc r>b rjc 



(3.1.20) 



where D ab = (D^D^ 6 and B% v = d^-duB^+gf^B^. Then, the generating functional 
for connected Green functions can be put as 



W 



■-irln [5^ b - 2gf cab B; u + £(D M D,) ab ] . 



(3.1.21) 



At first order in £ and second order in B fields, this is expressed as 

1 



W 



W + ZC A g z tr 



-AB; V AB; V - 2AB; v AB a vX Ad x d, 



+ 0(^,3"), (3.1.22) 



where as usual □ = d^d^ and A = — □ 1 . We can write the renormalized expression of the 
first term of 03.1.221) as 



(A) 



~ J d'xd'y B« u {x)BlM^\R 



(3.1.23) 
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whereas the second one is of the following form 

(B) = -2 J d 4 xd 4 yd 4 u (d u x d;A ux )B; u (x)B a uX (y)A xy A yu 

= -2 J d 4 xd 4 y B; u (x)B a uX (y)A xy d x x d; j d A u A xu A yu . (3.1.24) 

In order to evaluate the latter expression, we must apply CDR in momentum space 

If 1 1 

d 4 u A„,Am, = , . / d 4 u 



(47r 2 ) 2 J (x — u) 2 (u — y) 2 
' ^ d\d A pd A q -l- e - ip{x - u) e- 



(4tt 2 ) 4 J * * p 2 q 2 

1 f d A p L e -M*-v) 



(4tt 2 ) 2 J * p* 

4(4tt 2 ) 2 J P p 2 
ln(x — y) 2 m 2 



4(4vr 2 ) 

= -A(x-y). (3.1.25) 

(B) = 2jd 4 xd A yB; u (x)B: x (y)(A xy cf x d;A xy )\ R . (3.1.26) 

Hence, remembering the CDR renormalization of Adxd^A (11.2.391) we can obtain the 
renormalized expression for (B) as 

(B) = ~\j d'xd'y B; u (x)B; u (y)A 2 \ R - JL J d^yB^B^d^A^ . 

(3.1.27) 

Adding up the two contributions we have 

(A) + (B) = J #x*v BZ u (x)B a uX (y)dZd* x A xy , (3.1.28) 

which can be written in a more familiar form at explicit second order in the B fields as 
d A xd*yBl v {x)B a vX {y)d^dlA xy = J d 4 xd 4 y(d,B a u {x) - d u B;(x))(d u B a x (y) - d x B a u (y)) 

xd;0 x x A xy + O(B 3 ) 
= - J d 4 xd% Bl(x)BZ(y){dffi - 5^DpA xy + 0(B 3 ) . 

(3.1.29) 

With this result we obtain the previously defined if to be 

if = -^0, f d A xd 4 y B*{x)B*{y){&pZ - 5, v U){UA{x - y)) . (3.1.30) 





Figure 3.3: Two-loop YM diagrams (a)-(e). 



3.1.3 Two-loop level 

Now we follow with the two-loop contribution to the background field self-energy. The 
relevant diagrams are those of figures [3731 and [3^1 ((a) to (k)). Diagrams (a) to (h) have 
nested divergences, whereas diagrams (i), (j) and (k) have overlapping divergences. 



Diagram (a) 

This diagram has the following bare expression 



-y) = -2 fl < r/ W / 8 */ ,, 7 ^ ^ ^ 
x(d v x A uv )A uv (d u x A xu ) , 



(3.1.31) 

which can be rearranged in terms of the integral I 1 

K B (U X ) = ~9 4 C 2 A 5 ab [Ad v {Ad^I x ) - d^AI 1 ) - AAd.dJ 1 ] . (3.1.32) 



In order to renormalize, we have to replace these expressions with their renormalized 
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i 

(f) 











values, arriving to 



Figure 3.4: Two- loop YM diagrams (f)-(k). 



nfif (2a) ( X ) | K 



32(4tt 



2\3 



-±ln 2 a; 2 M 2 - |lna; 2 M 2 
^□-^ ^ 



±ln 2 x 2 M 2 + ±Mnx 2 M 2 
+<W □ □ ^ 5-5 



+ . . . 



(3.1.33) 



Diagram (b) 

This diagram is of the form 



n 



BB ab 
fii/ (2b) 



(x-y) = g 2 f ace f bed J d 4 ud 4 v A 2 



-28» x D* p + 8 Xp {dl-dl) + 25» p Di 



xUf a A ^(u-v)A 



-28 ua Dl + 8 aX {d; 



26 vX Dv 



A 



xy ; 



where Yl^^(x — y) = 5 ab Ii^^{x — y) is the one-loop correction to the quantum gauge field 
propagator. Its bare and renormalized expressions are found in section IC.2.11 of appendix 
[C| where it is used to obtain the leading term of the expansion of the function that takes 
care of the running of the gauge parameter in the RG equation. It has to be noted again 
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that, in contrast with dimensional regularization, the renormalized one-loop expression for 
the quantum gauge field propagator flC.2.17j) can not be used in the two-loop diagram. The 
reason is that the indices of the one-loop insertion the will be contracted in a second step, 
and one of the rules of CDR is to make first all the index contractions before performing 
the renormalization. Hence, only the bare one-loop contribution (IC 2.161) can be inserted. 
Therefore, expanding (13.1.341) we find 



n^f (26) ( x ~~ y) 



-g 2 C A 5 ab J d^v -Adffi [A XU UX (1) (u - v)A vy A xy ] 



tAA 



u 



V^)A V yA X y 



tAA 

+4D [A xu U^ {1) (u - v)A vy A a:y \ 

tAA 



v)A vy {dl- dl)A 



ftcr (1) 



, U 



V 



A V yA X y^ . 



If we use the bare result (1C.2.16I) for straightforward operations lead us to write 

this in terms of the previously defined 7° and a new integral expression of the form 

I°^(x -V)=f d A ud\ A xu A yv (A uv dffiA uv ) . (3.1.34) 

So, we have 



n 



BB ab ( } 
pv (2b)\ x ) 



R 



g 4 C 2 A 6 ba [-24d,d a (Ad u d a I°) + Ud^d v (ADl°) + 32d^(AI c v 
+ 125^8^8^1°) - 16^n(AD/°) - 165^8 p d a (AI° pa ) 
+ 208^A8 u ai°) - 20A8 IM 8 u ai° + l2n(A0 M $,J°) - 16D(A/J 1/yj ,, 

(3.1.35) 

It is clear that, as Dl° = —J 1 , with the renormalized values found for I 1 we can obtain 
all the expressions made up with ni°. In appendix ID. 21 we study the renormalization of the 
rest of the expressions made up with J° and I® that appear in this diagram. It is found 
there for Ad^d^I and AI® U the following renormalized values 

\nx 2 M 2 



[Ad,d u I°] R (x) 



32(4tt 2 ) 3 
1 

32(4tt 2 ) 3 



8pd u - 



x A 



±ln 2 x 2 M 2 + \lnx 2 M 2 



+ 



Unx 2 M 2 -±lnx 2 M 2 
dpd^ = + 5^a- " 



x. 



+ 



(3.1.36) 



With this, it is easy to arrive at the following renormalized expression 



n^f (2b)( x ) 



2 sab 



R 



9'C 2 A 5 
32(4tt 2 ) 3 



■f ln 2 x 2 M 2 



f lno; 2 M 2 



x z 



f ln 2 x 2 M 2 + £lnx 2 M 2 

+ <V nnJ 5-5 



+ ... 



(3.1.37) 
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Diagram (c) 

This diagram is easily renormalized as its expression is 

n2f(£)(x) = -g^f^ r d f My fy ec 5 ^ A 3 

= —g^S^A 3 



R Q A C 2 A 5 a \ ±lnx 2 M 2 

4 k^r^ aa ^^- (3L38) 



Diagram (d) 

This diagram is similar to the previous one, and we find 

TrBBabf \| ^ n 4 r 2 X ab X A 3 

1 V (2d) \ X ) I _ 2 9 A °^ iX R 

g 4 C 2 A 5 a \ -Unx 2 M 2 

Diagram (e) 

The bare expression of this diagram is 

n2T ( £)(s) = -\gY cd f b9e Jd 4 uA x ^-2S^D x p + S pa (^-d;) + 2S w D^ 
xA xu [r x f^ d (5 pX 5 £a - 5 p J eX ) + f C9X f xde (5 p J £X - 5 p£ 5 aX ) 

+f cdx f xe9 (S p£ 5 Xf7 - 5 pX 5 £fT )]A yu [-25 U£ Dl + 5 £X (d«- k) 
+26 vX Dy\A yu , (3.1.40) 

which, making all the index contractions can be written as 

n2T ( £)(*-y) = - s^n) J d 4 u A 2 xu A 2 yu . 

The renormalized expression of the integral is easily obtained as 

^ \n 2 x 2 M 2 



f ,4 11 7T 
/ d 4 u -r— -> □- 

7 X-M 4 M 4 4 



(x-w) 4 w 4 4 x 2 

so that 

n,%t)(-)U = 8^53^^(9A-^°)°— ^— + ■■■ (3-i-4i) 
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Diagram (f) 

This diagram is of the following form 

nSf ( g)(x-i/) = 2g i rv xdf f^f be} J d\/^ u {d;Ky){di-di)^ 

+2 g*f°f°f«*f»f*f*°* J d 4 u A xu (d; -d;)A xy (d:A yu )A 



*xy 



uy ■ 



Operating, this can be written in terms of I 1 , which allows us to write 



yiBBab . . 



)\ R = -g A C 2 A 5 ab [-2d,(AdJ l )+AAd,dJ 1 ] R 



g A C\5 ab 
32(4tt : 



2^3 



\\n 2 x 2 M 2 -\\nx 2 M 2 
d u d v n* 9 



-\\n 2 x 2 M 2 - 1 ±\nx 2 M 2 

+ S^nn-z 



+ . . . 



(3.1.42) 



Diagram (g) 

Contracting the indices of the bare expression 



n2f ( g) (x - y) = -2g*r*f*f'f«*ff< I d'u 5, a A xu {d u x A xu )A uy [-28 vX D* 



*xy 



(3.1.43) 



it is easy to write this diagram in terms of I 1 , which implies that the renormalized form is 



,4>nr2 tab 



g A C 2 A 5« b 



32(4tt 2 ) 3 



-d^AdJ 1 ) - Ad.dJ 1 - S^dxiAdxI 1 ] 

fUnVilf 2 + ±flnx 2 M 2 
d„d v U^ ^ 



R 



X^ 



-^\n 2 x 2 M 2 -^\nx 2 M 2 
+ <V DD — ^ 



x" 



+ ... 



(3.1.44) 
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Diagram (h) 

From the bare expression 



n^f (2h) ( x ~ y) 



-g* J d% - SM + f adx f xfc 



cdc 



A xuf 



A( e ) fbfe 



-25 vX DV + 6 eX { 



26 ve D y x 



A 



xy , 



where A^ A^ = (d^A^^A^, evaluating all the index contractions we can also express 
this contribution in terms of I 1 . Hence, the renormalized result is 



(2h) ( x ) I R 



-g 4 C 2 A 8 a < 

g 4 C 2 A 8 ab 
32(4tt 



45 

-d^AdJ 1 ) - 27Ad,8J 1 - dd.MAdJ 1 ] 
%ln 2 x 2 M 2 + 2 ±\nx 2 M 2 

f) f) n- - 



J R 



-%\n 2 x 2 M 2 + ^\nx 2 M 2 

+ S^nn-* ^ 



x- 



+ 



(3.1.45) 



Diagram (i) 

This diagram and the two following ones have overlapping divergences. In order to renor- 
malize them, we will make use of the list of integrals obtained in section 11.3.21 
The bare expression for this diagram which is 



n 



BB ab 
[XV (2i) 



X 



-y) = -g 4 f aSc f 9 ^f bde f 9ef d A ud A v AUdl-dl)(d u x A 



*uy) 



x {di- dl)A yv {dlA vx )A 



(3.1.46) 



Expanding this expression, we can write it in terms of the H integrals defined in (11.3.481) 
The bare contribution is then found to be 



n 



BB ab . 
[iv (2i) 



(x-y) = --g^C A b ab [ d£d»H[l, d x ; d x , 1] - 2d x H[l, d x ; d x d„, 1] 

-2d y v H[l, d x d p ; d x , 1] + 4ff[l, d»d x ; d v d x , 1] 



(3.1.47) 



At this point, we have to straightforwardly use the list of overlapping divergences of 
section 11.3.21 In concrete, with (11.3.511) . (I1.3.54p . and fll.3.64[) we find the renormalized 
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result to be 



R 



32(4tt 



2\3 



-^ln 2 x 2 M 2 -§lnx 2 M 2 

x 2 



±\n 2 x 2 M 2 + 2 %\nx 2 M 2 

+ V DD - ^ 



ar 



+ ... 



(3.1.48) 



-25^ + 5^-81) 



(3.1.49) 



Diagram (j) 

The basic form of this diagram is 

K B 3)( x ~y) = ig 4 r s r 9d f hde f Se9 f d\d A v& 

+ 25, p Dl\ {d u p K y ){dl- $)A^A U „)A 

Evaluating the index contractions this becomes 

nSf ( g)(x -y) = -g 4 C\5 ab [ -4dZd£H[l, d,d x ; 1, 1] - 4%dZH[l, d x ; 1, 0„] 

-4d'tH[l,d,d x ; l,d v ]-4d%H[l,d x ; l,d,d u ] 
-Ad x u H[l,d^d x ; l,d x ] + dffiH[l,d x ; l,d x ] 
-4if[l,fyd A ; l,cU> A ]+4D#[l,c> M ; 1,9,] 
+4D#[1, 1 ; 1, d^} - 2d x u H[l, dp ; 1, □] 

l,n]-4ff[l,ft,D; 1,^] 
-2^//[l,^D; 1,1]] . (3.1.50) 

From this list, the expressions that have a □ (remember □A(x) = — can be written 

in terms of the integral J 1 , and their renormalization is straightforward. The rest of the 

integrals can be found in the list of integrals with overlapping divergences. So, we arrive at 
the following renormalized result 



U fiv(2j)( X )\ R ~ 



g A C 2 A 5 ab 
32(4tt 2 ) 3 



n n ±ln 2 x 2 M 2 + ±lnx 2 M 2 
f) f) □— - 

M x 2 

±ln 2 x 2 M 2 - §lnx 2 M 2 



x A 



+ 



(3.1.51) 



Diagram (k) 

In order to obtain all the contributions that form this diagram, the Mathematica package 
'FeynCalc' was used, so that all the index contractions were performed by the computer. 
The output of this process are the final relevant expressions that need to be renormalized. 
The contributions shown here are those that have a divergent part, omitting those terms 
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that are finite. 



n^f (2k) ( x - y) 



~-Cl5 ab 



[ +168 tMl/ d x x d*H[l, d x d a ; 1, 1] - 20c^tf[l, d,d x ; 1, 1] 

-124cP9*tf[l, <9 A ; 1, 9J + 72d x x H[l, d,d x ; 1, 9„] 
+565, u d x x d x a H[l, d x ; 1, <9 CT ] - 72d x u H[l, d,a ; 1, 1] 
+20^ff[l,D; 1,1]-144//[1,9 M D; 1,9,] 
+729*#[1, □ ; 1, 9,] - 72%H[1, d,d x ; 1, 9 A ] 
+Md?dZH[l,d x ; l,d x ]-72H[l,d,d x ; l,d v d x ] 
+40Off[l, 9 M d„ ; 1, 1] + 32D#[1, d M ; 1,9,] 
+I6^n#[l, □ ; 1, 1] - I6^n//[l, <9 A ; 1, 9 A ] ] 
+ (finite terms) . (3.1.52) 



Proceeding in the same way as in the two previous diagrams, we can easily found the 
renormalized form of this contribution to be 



tiBB ab ( 

1 V (2k) \ x 



g A C 2 A 5 ab 
R ~ 32(4tt 2 ) j 



-f ln 2 x 2 M 2 -f lnx 2 M 2 
d&U—* 



ar 



27 ln 2 a; 2 M 2 + f lnx 2 M 2 



+ ... 



(3.1.53) 



Two-loop final results 

In order to obtain the total two-loop renormalized contribution to the background gauge 
field self-energy, we only have to add all the renormalized expressions for the diagrams that 
we have obtained. So, we have 

K B 3(*)\n = + (3-1-54) 

By fulfilling one-loop CDR rules, we have fixed the renormalization scheme a priori, 
which implies that one-loop local terms have defined values. Hence, as is imposed by gauge 
invariance, we have obtained directly a transverse result. 



3.1.4 RG equation 

With the previously obtained expressions for the one- and two-loop corrections of the back- 
ground gauge field propagator, we can obtain the first two coefficients of the expansion of 
the beta function of this theory. If we define 

T* u Bab (x) = (d,d u - 8, u D)5 ab ^ 2 \x) , (3.1.55) 
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then, with the one- loop contribution (I3.1.19p . the gauge fixing renormalization (I3.1.30P and 
the two-loop contribution (13.1.541) . the effective action for the background gauge fields is 



r (2), 



X 



1 p . , 11C A \nx 2 M 2 
—5{x) + ^——r-D- 



g* ■ ' 487r 2 (47r 2 



x- 



72tt 



8tt 2 



2 ^ 2 lnx 2 M 2 

□ s— + ... 



x- 



2(4tt 2 ) 3 

With this definition, the equation we need to consider is 

r (2) 



(3.1.56) 



5=0 



(3.1.57) 



Notice that 7^ is the coefficient that takes care of the running of the gauge parameter. 
The expansion of this function to order g 2 is obtained in appendix O To do so, we consider 
there the one-loop RG equation for quantum gauge fields. With this, we find for 7^ 



7£ 



5C A 2 ^ 
;9 + 



24vr 2 - 



(3.1.58) 



Also notice that if the background gauge field is redefined as B' = gB, this implies 



7s = (the charge and background field renormalizations are related: Z g 



-l/2x 



' B 



So, 



with (I3.1.56p . (13.1.571) and (13.1.581) . we evaluate the first two coefficients of the expansion of 
the beta function to be 



Pi 



(3i9 3 + (329 5 + 0(g 7 ) 
IIC A 
48vr 2 
17C A 



24(4vr 2 ) 2 ' 

These results agree with those previously obtained in the literature HU S3 EE]. 



(3.1.59) 



3.2 N = 1 Super Yang-Mills 

In this section we consider the two-loop differential renormalization of the supersymmetric 
extension of the previous model, N = 1 Super Yang-Mills |49, 50J. With this calculation 
we revisit an old controversy: the origin of higher-order perturbative contributions to the 
beta function in supersymmetric gauge theories [H [56J, [631 Ml [65]. Differential renormal- 
ization has one important advantage over usual renormalization methods (as dimensional 
reduction), as in this case we have UV and IR divergences. With dimensional methods 
both renormalizations mix (we need to subtract the IR part in the final result), but dif- 
ferential renormalization clearly distinguishes between UV and IR divergences as both are 
renormalized with independent scales. 
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3.2.1 N = 1 Super Yang-Mills model 



As is detailed in appendix |Aj in order to formulate a supersymmetric gauge theory we can 
follow two different approaches. In the first one, named chiral representation, we begin 
by considering a multiplet of unconstrained gauge superfields (V = V a Ta, with Ta the 
group generators), which are a generalization of the results found after studying the off-shell 
representations of the linear free theory. We use these gauge superfields to construct covariant 
derivatives = (e~ 9V D a e 9 V, D a , — i {V^ , V£j) that allow us to obtain gauge invariant 
expressions. On the other hand, with the second approach, called vector representation, 
we begin by considering covariant derivatives (termed = (V^, V^, V£ d )) that, after 
imposing covariant constraints on them, can be expressed in terms of prepotentials. With 
both approaches we find field strengths defined in terms of an spinorial field. In particular, 
for chiral representation we have 



W a = iD 2 (e~ 9V D a e 9 V) 
Wa = e~ 9V W a e 9 V = e~ 9V { 

which allow us to define a gauge invariant action as 



■W a ) 



o9V 



(3.2.60) 



So 



1 

7 



tr / <Fxd 2 9 W 



— \tr [ d 4 xd 4 9 (e- 9V D a e 9 V)D 2 (e- 9V D a e 9V ) (3.2.6V 
2g 2 J 



For vector representation, the field strength is defined as 

with 

>p_ ■ 

and with this we write the gauge action as 



V" iX7 v ■ 
v a > L v pp 



v a > v p 



ap 



(3.2.62) 



On 



9 Z 



tr / d A xd 2 6 W A 



29 



-tr J d*xd 2 6 Q [V v& , {VI,V:}]j . (3.2.63) 



Quantization of a supersymmetric gauge theory is also discussed in appendix [A] We have 
to add to the action a gauge-fixing term that depends on a gauge parameter a (Sgf) and 
anticommuting chiral ghost fields c, d [Spp). Their explicit expressions are 



Sgf 
Sfp 



tr 



-tr f d 8 z (D 2 V)(D 2 V) 
a J 

d 4 xd 4 e (d + d)Li 



±9V 



[c + c) + cothLi gV (c 



L X Y 



[X,Y] . 

(3.2.64) 



One of the relevant features of the supergraph techniques applied to this theory is the 
appearance, along with the usual on-shell infrared divergences of Yang-Mills theory, of ad- 
ditional infrared divergences due to the form of the gauge propagator in a general covariant 
gauge [HU [52], [531 EH [55]. This can be clearly seen if we consider the expression for this 
propagator in momentum space [55] 



A(k) 



1 + (1 - a)(D 2 D 2 + D 2 D 2 )k~ 2 



-8\9 - 6') 



(3.2.65) 
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As the leading term in D 2 D 2 + D 2 D 2 is constant when k 2 — > 0, this propagator goes as 
\jk A at small k, which is the origin of the infrared divergences. Although Feynman gauge 
(a = 1) seems to be the solution, this is not the case: the one-loop correction of the gauge 
propagator takes us out of the Feynman gauge; hence, when we consider two-loop diagrams 
that have as an insertion the one-loop corrected propagator, the infrared divergences reappear 
|55j . With dimensional methods this situation represents a severe problem, as both UV and 
IR divergences are renormalized with the same dimensional parameter e (d = 4 — 2e). So, 
we have to subtract the IR contribution^]. However, life gets simpler if we use differential 
renormalization. As we have seen in section ll.l.3[ with differential renormalization UV and 
IR divergences are renormalized with different and independent scales. 



Background field method 

We discuss the application of the background field method to supersymmetric gauge theories 
in appendix [B] As is explained there, due to the non-linear gauge transformations of SYM, 
a linear quantum-background splitting is unsuitable. The accurate spitting is achieved if we 
replace e gV with 

e gv ( s pm) = e n e 9V e fi ; (3.2.66) 

where V is the quantum gauge superfield and f2 is the background prepotential. It is 
worthwhile to mention that we have redefined the usual background field B (17 = Cl = ~B) 
as gB — ► B. This splitting implies that we write the covariant derivatives in a quantum-chiral 
but background-vector representation as 

V Q = e- 9V V a e 9V , Va = V A , V Q « = -i {V Q , V^} , (3.2.67) 

where V a and V„ are the background covariant derivatives. Relative to the quantization 
of the theory, we remark one of the particular features of Super YM: the appearance of the 
Nielsen-Kallosh ghost [58]. In the usual quantization procedure we gauge-average with a 
simple exponential factor; however, if we use a more complicated function e.g., exp J f'Mf 
with M an operator, we have to normalize the procedure dividing by detM. As we have 
this situation when we use background chiral superfields, we introduce a new ghost field 
b (Nielsen-Kallosh ghost), which have opposite statistics to / and allow us to properly 
normalize the gauge-averaging procedure. As this new field does not interact with the 
quantum fields and enters quadratically in the action, we find that it only contributes at the 



2 In [SS] this was achieved at two-loops by choosing a non-local gauge-fixing term that cancels exactly the 
contributions that takes us out of Feynman gauge. In [56] the procedure was to define a R— operation [57] . 
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one-loop level. Hence, all the relevant terms that form the split Super Yang-Mills action are 



Sym 
Sgf 
Sfp 
Snk 



1 



tr / d 4 xd A 6 (e- 9V V a e gV W 2 (e- gV V n e 9V ) 



[1+Otr j M(W)(VV) 



tr / d A xd 4 6 



cc - cc + -(c + c) [gV , c + c] + . . . 



1 + £)tr / d^xd^O bb . 



(3.2.68) 



Notice that we have redefined the usual gauge parameter a as - = 1 + £. Hence, as is 
discussed in appendix |Bl we have a background effective action of the form 



TIB] 



S [B]+T^ + 

"1 



+cxpS. 



int 



' 5 


5 


5' 




Jj' 






exp 



-In '•/ ./::.•/ 



(3.2.69) 



J=j=j=0 



where Sq[B] is the "free" part of the background action, stands for the one-loop contri- 
bution in the gauge £ and J, j and j are the sources. □+ and □ are operators defined in 
appendix 03] as 



□ 



□_ 



□ 

□ 



iW a V c 
iW a V r 



-(V a W r 



(3.2.70) 



To apply supergraph techniques to this theory, we can follow two procedures: In the first 
one we expand all the background covariant derivatives in terms of the explicit background 
connections V a = D a — iT a , so that we can employ usual D-algebra. This procedure is 
applied in [SS], and we will use it in the one-loop level. However, another approach was 
considered in [5^1 [60] . In this case the spinorial background connection is not explicitly 
extracted from the background covariant derivative, and covariant D-algebra is used in the 
diagrams. At the end, we have all the diagrams expressed in terms of the background 
space-time connections T a6l or field strengths W a . Thus, as we do not have explicit spinor 
connections F a (which are of lower dimension), the diagrams are more convergent and fewer 
in number. For the two-loop calculations we will follow this procedure. 



3.2.2 One-loop level 

We will proceed with the one-loop case. As the previous examples considered, we will obtain 
the one- and two- loop corrections evaluated in Feynman gauge (with our conventions, £ = 0). 
Therefore, at this level we have not only to consider the background gauge field self-energy 
correction, but as we will have to take care of the running of the gauge parameter in the RG 
equation (j^d/d^), we will obtain the linear contribution in £ of the one-loop background 
gauge field two-point function. 
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Background gauge field self-energy 



/TOP 



Figure 3.5: One-loop background gauge field two-point function contribution. Thick lines 
correspond to external background fields and thin lines represent ghost propagators. 

To obtain the one-loop contribution, we begin expressing the covariantly chiral ghost 
fields in terms of ordinary chiral fields by 



c ^ e B/2 ce -B/2 



(3.2.71) 



where B is the background gauge field. Hence, the one-loop relevant interaction terms are 
(see appendix [Bj 



tr / <Fx<P0 



-~V{n - W"V„ - W"V n ) r + c'Bc + c'Bc + bBb 



(3.2.72) 



Note that in the previous expression, from the part which corresponds to the interaction 
between the quantum and background gauge superfields, it is clear that we do not have 
enough covariant derivatives to obtain a non- vanishing contribution to the background gauge 
field two-point function (remember that at least we need two D 2 and two D 2 ). Hence, we 
conclude that the only non-vanishing contribution comes from ghost superfields (as they are 
chiral superfields, we have at the vertices additional superspace covariant derivatives [3"8]). 
With the definition of the superspace propagator P^ = d^A^, the total contribution of the 
ghost superfields is straightforwardly evaluated as [371 EH] 



-pi loop 



3C A S ab 



2 

3C A 



d 8 z x d 8 z 2 B a (z 1 )B b (z 2 ) 



D\P 12 Dl 



D 2 2 P l2 D\ 



J d 8 Zl d 8 z 2 B a ( Zl )B a (z 2 ) [D 2 2 D 2 2 P 12 ] [D\D\P X2 \ . (3.2.73) 



Applying the identity f)2.2.5ip . we can write this expression as 



■pl loop 



2 

i3C A 



d 8 Zl d s z 2 B a ( Zl ) [D 2 D 2 B a ( Z2 )] P 12 [D 2 2 D 2 2 P 12 ] 
d 8 Zl d 8 z 2 B a ( Zl )B a ( Z2 )P 12 [UD 2 2 D 2 2 P 12 ] 
d 8 Zl d 8 z 2 B{z x ) [D«D a B a (z 2 )] P 12 [dlJ)\D\P X2 \ 



(3.2.74) 
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And finally, with the usual superspace 5-function property 6uDlD 2 5i2 = £12 and the 
identifications x\ = x, X2 = y we arrive at 



pi loop 



2 
2 

i3C A 



, 2 

xy 



d A xd A yd A 9 B a (x,6) [D 2 D 2 B a (y,6)] A] 
d 4 xd A yd A 6 B a (x, 6)B a (y, 9)A xy UA xy 

d A xd A yd A 9 B a (x, 6) [D«D a B a (y, 9)] A xy d v adl A xy . (3.2.75) 



Applying CDR rules, we find the renormalized form of this expression to be 



pi loop 



3( ' 1 f d A xd A yd A 6 B a (x,6) [D a D 2 D a B a (y,6)} n Hx 



16(4vr 2 ) 2 7 » " v *' J1 (x-y) 2 ' 

(3.2.76) 

where we have used the superspace derivatives identity D 2 D 2 -\-(i/2)d a aD a D a = l/2D a D 2 D a . 



Effective action in a generic gauge 

Now we proceed with the additional result that we have to obtain in order to deal with 
the running of the gauge parameter: the contribution to the background effective action of 
quantum gauge fields evaluated in a generic gauge. We have to follow a procedure similar 
to that used in the Yang-Mills case. First, we have to consider the one-loop effective action 
contribution at second order in background gauge fields evaluated in a generic gauge. Then, 
we have to expand this in terms of the gauge parameter £, retaining the linear part. The 
reason for this is the same as in the non-supersymmetric case: in the background gauge field 
RG equation, after considering the term that takes care of the running of the gauge parameter 
7£(g)<9/<9£, we will impose Feynman gauge (£ = 0); hence, the only relevant term for us in 
the ^-expansion is the linear one. As in the Yang-Mills case, to perform this calculation 
we consider a functional approach. The quadratic action that we have in this case implies 
that we find the following contributions from V fields and Nielsen-Kallosh ghosts b to the 
one-loop effective action [55] 



eff 



— tr In 
2 



□ 



v 2 v 2 



tr In 



□_ +£V 2 V 2 



(3.2.77) 



which is written in terms of the operators defined in appendix [B] as 



□ 



□ 



□ 



□ 



iW a V c 

lW a V r 



iW ' V 

i 



V a W a ) 



n-iW V* - -(V W, 



(3.2.78) 



Expanding (13.2.771) in £ we find 



\ ff = — -tr Inn + tr In □_ + T ? 



1 £ 
— tr\nn + tr\nn_ + -Tf + 0(£ 



(3.2.79) 
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and the linear part is 



.(i) 



tr 



, tr [i(V 2 V 2 + V 2 V 2 ) 
_□ 



+ tr 
+ tr 



l(n-n + )— V 2 V 2 



□ 



□_ 



tr 



□ 



- I -iW a V a + -(V~W A ) I — v 2 v 2 



□_ 



+tr 



□ 



^ -WT, + -(V a W Q ) — VT 



1 



□ 



(3.2.80) 



where in the second step we have used the property irCII 1 V 2 V 2 = trD^ 1 V 2 V 2 [59J. At this 
point, applying □I 1 V 2 V 2 = V 2 D^ 1 V 2 [59J, the anticonmutative nature of the covariant 
derivatives (V Q V 2 = 0) and the Bianchi identity ^ a W a = -V"^ we arrive at 



% 

—tr 
2 



I 

-tr 
2 



(3.2.81) 



So, considering the inverse of the operators, 



1 


l 


% 




= □ - 


-□ 


i 


l 


% 


□7 




□ 



-(V a W a ) I = 



w 6 v d + i(v Q ^ 



1 

□ 
1 



□ 



+ 



(3.2.82) 



where □ = l/2V aa V Q , (i , the contribution at second order in the background gauge fields is 



1 

— tr 
2 

1 

+ -tr 
2 



□ 



o 



□ 







1 ,^ 



72^2 



□ 







□ 







W Q V a + -(V a W c 



□c 

-v 2 v 2 

□ n 



C(fi 3 



(3.2.83) 



with D = l/2d aa d a a being the usual d'alembertian. Hence, as the terms that corresponds 
to ^ a D 1 V Q V 2 V and r^'VaV V 2 have not enough covariant derivatives to give a 
non-vanishing result at second order in the background gauge fields (remember V Q , V = 

-iVaaV" + iW a ), the contribution is found to be (once we have imposed the Bianchi 
identities) 



1 

--tr 
2 



u u u 



0(B 3 ) 



-tr I d 8 Zl d s z 2 d 8 z 3 [D a W a {z 2 )\[DPW p {z^)\ [D\D{P 12 ] P 23 P 13 + 0(B 3 ) . 



(3.2.84) 
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After simplifying this result with the usual superspace 5-function identity (IA.5.541) and 
using the identifications X2 = x, £3 = y and x± = u we have 

lf } = l -tr J d i xd 4 yd i 9[D a W a (x,e)][D' 3 W0(y,9)}A xy J d\ A xu A yu . (3.2.85) 

This expression is IR divergent, and it has been evaluated in (13.1.251) . The renormalized 
result found there is 



1 ^ rxr«t„ mw MnM 1 ~ V) M IR 



tr I d 4 xd 4 yd' 2 9 W a (x,9)W a (y,9)n^ — y \^ lu +0{B 3 ) , 



(4tt 2 ) 2 J * y ' ' ' (x-y) 



(3.2.86) 



where in the last step we have applied the identitjH 

d 4 xd 4 yd 4 9 [D a W a (x, 9)} [D^Wpiy, 9))f{x - y) 

= J d*xd 4 yd 2 9 W a (x, 9)W a (y, 9)Df(x - y) + C(5 3 ) . (3.2.87) 

Hence, the linear term of the expansion in the gauge parameter of the one-loop effective 
action evaluated at second order in the background gauge fields is 

r € = -j^ytr J #x*y*0 W"(x,9)W a (y,9)n ^ X -y^ . 

(3.2.88) 

3.2.3 Two-loop level 

We proceed now with the calculation of the two-loop contribution to the background gauge 
field self-energy. As we have stated previously, we will use here covariant D-algebra, which 
simplifies and reduces the number of the diagrams we have to consider. We begin with the 
pure contribution of quantum gauge fields, and leave ghost contributions for a later section. 

Quantum gauge field contribution 

In order to obtain these contributions, we have to expand the gauge action So of (13. 2.68ft 
and obtain the different interaction terms. With them we are instructed by the background 
field method and covariant Feynman rules to consider diagrams with external background 



3 J d 4 xd 4 yd 4 9(D a W x a )(D W y p)f{x - y) = - J d 4 xd 4 yd 4 6W x a (S p a D 2 W%)f(x - y) 
= - J d 4 xd 4 yd 2 6W xa {[D 2 ,D 2 ] W°)f(x~y)+(D(B 3 ) = J d 4 xd 4 yd 2 9W^(aW y a )f(x - y) + 0(B 3 
= J d 4 xd 4 yd 2 eWZ t W ya nf(x-y) + 0(B 3 ) 
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Figure 3.6: Two-loop contribution to the background effective action. Wavy lines correspond 
to quantum gauge field propagators. 



fields and covariant propagators □ 1 inside loops. From this expansion, the requirement of 
having at least four V and four V to get a non- vanishing contribution imposes that the only 
relevant interaction term in our problem is 



9 + 
—tr 

2 



v { v a v, v 2 v„y 



(3.2.89) 



which allow us to construct a vacuum diagram like the one shown in figure 13. 6[ where the 
wavy lines correspond to covariant propagators □ 1 . Note that we have implicit interactions 
with the background field in the covariant derivatives and the covariant propagators. 

To evaluate this diagram, we have to rearrange the covariant derivatives at the vertices, 
where their explicit expression is {ig/2)f abc V a V a V b V V a V c . At the left vertex we choose 
a given configuration of derivatives which, after using the commutation relations of the 
covariant derivatives and integration by parts, can be rewritten as 

l -jLf abc V a V a V h V 2 V a V c = % j-f abc (-2V a V 2 V h V 2 V c + V a V a V b iV a& V & V c 

-V a V a V b [iW a , V] c ) . (3.2.90) 

Once we have this fixed arrangement, at the right vertex we have to choose the six 
permutations that are possible. Then, we integrate by parts in each of them so that one 
specific line if free of any operators. This implies that this vertex is written as [59J 

l -lf abc V a (-2V 2 V b V 2 V c + 2V a V b V 2 V a V c - 2vV 6 v 2 v d v c - v a v b zV ad v d V c 

+V"V b iV a6 ,V a V c - iV a6t V h VaV a V c + iV a V b [W a , V} c - 2zVV [Wa , v] c ) . 

(3.2.91) 



It can be shown that most of the different combinations of terms at each vertex either 
cancel by not having enough covariant derivatives, produce pairs of divergent contributions 
that cancel each other after using the Bianchi identity V a W a = -V W A or give finite 
Feynman integrals [59]- The only non- vanishing divergent contribution that we find is given 
by the first terms of (13 . 2. 90[) and (13.2.911) . At this stage, we can obtain explicit background 
gauge fields by expanding each of the covariant propagators to second order in background 
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gauge fields as 



1 1 
□ ~ □ 



A ( W a -V c 



□ 





v«, - 


+ w & 


v d , — 


) 




' □ 




' □ 





i 



(3.2.92) 



Notice that the first term of the expansion gives no contribution, as when combined with a 
similar term from another line is finite, and when combined with a space-time background 

connection T a J^ it gives a contribution of the form of r(V a W Q + W^Vg) that clearly 
vanishes by the Bianchi identities [59] . 




Figure 3.7: Diagrams corresponding to the expansion of □ . Thick lines correspond to 
external background fields. 



Let us consider first the third term of (13.2. 92p . which generates a diagram of the form of 
diagram (a) of figure 13.71 Notice that all of these contributions have a common symmetry 

As we have two explicit background 

= Dn , with 



factor of g that we will take into account at the end 
field strengths, at second order in the background gauge fields we have 



□o = (l/2)d aa d a a the usual d'alembertian. Therefore, we have an explicit expression of the 



form 



2 loop 



3g 2 C 2 A tr j d 8 Zl d 8 z 2 d 8 z 3 d 8 z 4 [VjP 12 ] [w a (z 2 ) V 2 a + W*(z 2 ) V 2 , 



V3P34 



P 



v 2 A v\p A1 



p. 



23 



(3.2.93) 



4 Recall that there is another implicit dependence in background fields, as we have □ = l/2V QQ V i: 
l/2(d aa — iT aa )(d a a — ir a &), where is the background space-time connection. 
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Due to the anticommutative nature of the covariant derivatives V 3 ^ V 3 = 0. Hence, we have 



T 2ioo P = _3 g 2 C ^ tr J d 8 z 1 d 8 z 2 d 8 z 3 d 8 z 4: \ylPi2]W 6l {z 2 )V 2 aP2zW l3 (z z ) 



V^ 3 V 3 P 3 4 Pu 



X 



vlvlPu 



3g 2 C 2 A tr / d 8 Zl d 8 z 2 d 8 z 3 d 8 z 4 [V 2 P 12 ] W a (z 2 )V 2 a P 23 W?(z 3 ) 



V 3/3 V3P 3 4 



14 



X 



v 2 vjp 4i 



(3.2.94) 



As can be seen, we have divided this expression in two contributions. When dealing with 
the first one, as W 13 is a covariantly chiral superfield (V^W 13 = 0) and taking into account 
the basic relation between the covariant derivatives {V^ , V Q } = i'Vaa, we find 



2 loop 



1 1.1 



3ig 2 C 2 A tr / d 8 Zl d 8 z 2 d 8 z 3 d 8 z A [V\P l2 ] W & (z 2 )P 23 W a (z 3 ) 



VLV3P34 



x 



v 2 v\p 41 



Pu 

(3.2.95) 



As we need four V and four V to get a non- vanishing result, and realizing that at sec- 
ond order in the background fields in this expression V Q and V Q q, commute (W a = 

— \ V° , V a(i ), it is clear that integrating by parts we obtain the non- vanishing contri- 
bution to be 

1*5"* = 3tg 2 C 2 A tr J d 8 z 1 d 8 z 2 d 8 z 3 d 8 z 4 P 12 W'"(z 2 )P 23 W a (z 3 ) 



x 



v 2 v\p lA 



+ £>(P 3 ) 



Pu 
(3.2.96) 



At this point, replacing the covariant derivatives by the usual ones (as we have already two 
explicit background field strengths), using the usual 5-function superspace identity and the 
identifications x\ = u, x 2 = x, x 3 = y, x A = v, we have 



, 2 

uv 



V\ l r = 3ig 2 C 2 A tr J d 4 xd%d 4 6 W a (x,6)W a (y,9)A xy J d 4 ud 4 v A xu {d y a6l A yv )Al 

+0(B 3 ) . (3.2.97) 

To obtain the second contribution to T 2 loop , we consider the covariant derivatives that 
are acting over P\ 2 and, integrating by parts, make them act over P 34 



T\T P = -3g 2 C 2 A tr d 4 xd 4 yd 4 6P 12 W a (z 2 )V 2a P 23 W^z 3 ) 



v 3 pvlv 2 3 p 3A 



x 



V^vJPu 



P41 

(3.2.98) 



Now, with the relation 



-V a( iV a + iW a , we conclude that this contribution 



vanishes, as we do not have enough covariant derivatives. 
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We now consider the contributions that come from the second term of (13.2. 92p . In this 
case we study the commutator between D _1 and the covariant derivatives, finding out 



V a , — 



— [V a , □]- 



(3.2.99) 



and 







Vd , — 







W a V a& + -(V aa ,W a ] 



1 

□ 



(3.2.100) 



Hence, for this case we also find diagrams of the form of diagram (a) of figure 13.71 As we 
have two contributions that are identical, except for having W a and interchanged, we 
detail the calculation of only one of them. 



It'r p = 3ig 2 C 2 A tr J d 8 Zl d 8 z 2 d 8 z 3 d 8 z 4 [V?P 12 ] W a (z 2 )P 23 



W a (z 3 )Vl a + -(V a ,W a (z 3 )) 



VjP 3 4 



P 



ii 



vlvlp u 



.(3.2.101) 



As when obtaining T 2 l 2 oop , we can integrate by parts the covariant derivatives acting over P 12 
to find 



-p2 loop 
1 2.1 



3ig 2 C 2 A tr \ d 8 z x d 8 z 2 d 8 z 3 d 8 z 4 Pi 2 W a (z 2 )P 23 



v\vIp u 



(3.2.102) 



that is an expression that can be directly simplify with the application the 5-function identity 
(IB.2.32P obtained in appendix [B] With the same identifications as in the T 2 loop contributions, 
we find the total T 2 loop contribution to be 



-p2 loop 
1 2 



-3ig 2 C 2 A tr J cPxdtycPe W a (x, 9)W & (y, 9) + W"(x, 9) W a (y, 6 
xA xy J d\d A v A xu d y a& A yv A 2 uv + 0(B 3 ) 

-J9 2 C> I d 4 xd 4 yd*9 [w a (x,9)d y a& W"(y,9) + W*(x,e)& cA W a (y t 



xA xy J d\d 4 v A xu A yv A 2 uv + 0(B 3 



(3.2.103) 
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Hence, the sum of the contributions T 2 loop and T\ loop (written in terms of the 1° integral 
expression defined in section H.1.3P is 

r 2 loop + r 2 loop = 3ig 2 C 2 Atr j # x tf y # e w <*^ 9)W & (y, 9) [Ad a J°] (x - y) 

W a (x, e)d v a6t W* + W"{x, 9)d y a . W a (y, 6 

x [AI°] (x-y) + 0{B 3 ) . (3.2.104) 

Once we have finished the study of the expansion of (13.2.921) . we proceed now to consider 
diagrams with background space-time connections I\ We begin by considering the expansion 
of the inverse of the □ operator to second order in T, which is obtained as 



.^^Atr I d'xd^yd'e 



□ Ur 



□ r 



□r 



■(d aA T adl ) + T a "d aA 
-(d a «T aQi ) + T a °d Q 



□ 



1 1 1 r-T,,,,; 1 



2 D f 



□ r 



□ r 



PP 



□r 



+ . . . 



(3.2.105) 



If we start with the third term of (13.2.1051) . what we have is a contribution of the form 
of diagram (a) of figure 13.71 So, this is written as 



t2 loop 



3g 2 C 2 A tr / d 8 z l d 8 z 2 d 8 z 3 d 8 z A [V?P i2 ] 



T a& {z 2 )^ + -{^T a ^){z 2 ) 



P 



2.3 



vlp 



3^34 



P 



11 



v\v\p lA 



3g 2 C 2 A tr I d* z x d* z 2 d* z 3 d 8 z± P 12 



8„ j8„ j8. 



^K + ^d P %^(z 3 ) 



r^(z 2 )d 2 Qa + -(d aa r aa )(z 2 ) 



p 



23 



V*V§i>34 



Pi 



14 



V4V4P14 



(3.2.106) 



where in the last step we have applied the same procedure as with T^l and T 2 , loop , inte- 
grating by parts the covariant derivatives that are acting over P\ 2 . Hence, after the usual 
identifications and some simple algebra, we find for this contribution 



-\2 loop 



■-g 2 c 2 A tr j d i xd i yd i er^(x,e)r^( y ,e)d:^ [ai°] ( x - y ) 



+3g 2 C 2 A tr I d 4 xd 4 yd 4 9T aa (x,e)T f3/3 (y,9)d^ 



Ad* 7° 



3g 2 C 2 A tr / d 4 xd 4 yd 4 9r a ^(x,9)T^(y, 



Acf d x 1° 



{x-y) 

(x-y) + 0(B 3 ) . 

(3.2.107) 

Let us now consider the first term of the D _1 expansion (13.2.1051) . At second order in 
the background fields, we find the relevant contribution to be obtained when we expand two 
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different lines, obtaining a diagram of the form of diagram (b) of 13.71 Explicitly, we have 



-p2 loop 
1 4 



-g 2 C\tr / d 8 z l d 8 z 2 d 8 z 3 d 8 z i [V?Pi 2 ] 



r-(z 2 )^ d + -(^r Qd )(^ 2 ) 



r,3/? (^)^ + -(^r /3/3 )(^ 4 ) 



vjp 41 



Pl3- 



^23] 

(3.2.108) 



In this expression we have only four V and four V, so it is obvious how we can obtain 
the non- vanishing contribution: we integrate by parts the covariant derivatives to move four 
of them (for example those that are acting over P 12 and -P34), and made them act on the 

— 2 — 2 

other four (V3P23 and V4P41 in our example). After some simple algebra, we find for this 
contribution 



-\2 loop 



-g 2 C\tr J d 4 xd 4 yd 4 9T^(x,6)T^(y,9)H[l,d a( ,; l,d p$ ] 

+ \g 2 C 2 A tr j d 4 xd 4 yd 4 eT a& (x,e)T^(y,e)^&; $ H[l,l; 1,1] ,(3.2.109) 



where we have used the H integrals defined in (11. 3. 481) . 

Finally, we consider the second term of the D _1 expansion. It is clear that this generates 
a contribution of the form of diagram (c) of figure 13.71 which is written as 



-p2 loop 
1 5 



--g 2 C\tr j d 8 z 1 d 8 z 2 d 8 z 3 d 8 z i [V?P 12 ] r aA (z 2 )5 8 3 r aA (z 3 ) 



VJP34 



xP 



11 



(3.2.110) 



After freeing Pi 2 of covariant derivatives and make them act over P34, the expression can 
be treated with the usual steps (5-function identity and point identifications). Hence, we 
find the bare expression to be 



-p2 loop 
1 5 



-g 2 C 2 A tr j d 4 xd 4 yd 4 6 T a& (x,9)T^(y,e)(2C aP C^) [□(AJ ) - d^Ad^I ) 
+AUI ] (x - y) + 0{B 3 ) . (3.2.111) 



Ghost contribution 

As the Nielsen-Kallosh ghosts only interact with the background gauge field and they enter 
quadratically in the action, it is clear that they only contribute at the one-loop level. Hence at 
two-loops we have contributions from c and d ghosts, which are proportional to the difference 
of the two graphs shown in figure 13.81 Expanding the propagators as for the quantum gauge 
field, we easily find that terms with two Y cancel [59]. For Vt^-terms, we obtain that the 
only divergent contributions come from factors acting on the same line (either V^Vd -1 V 
or W [V , D -1 ]). However, all of these terms can be listed and showed to cancel each other 
or produce a combination that vanishes once we impose the Bianchi identities [59]. Hence, 
we conclude that the two-loop ghost contribution vanishes [551 EH] . 
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Figure 3.8: Two-loop ghost contribution to the background effective action 



Total renormalized contribution 



As we have seen, all the divergent contributions have been written in terms of the integral 
expression we have defined as 1° and two of the H overlapping integrals listed in section 
II .3.21 Hence, we only have to directly use the renormalized results previously found and 
replace the bare expressions with the renormalized ones. Thus, we find for the first two 
contributions 



2 loop\ 



3ig 2 C 2 



tr / d*xd 4 yd*9 W a (x, 6)d y adl W & + W & (x, 6)d y aA W a {y, 9) 



x 



[A/°U*-y) 



\R 



^£^tr I d'xd'yd 2 Q W a (x,6)W"(y,6)d: o 



ln(x - y) 2 M 2 
(x - y)< 



+ 0(B 3 ) . 
(3.2.112) 



Appliying the Bianchi identities we find an useful relation to simplify this result and express 
it in a explicit gauge invariant form. Let us consider an expression of the form 

J d 4 xd i yd i eW a (x,e)(d y a .w\y,e))f(x-y) . (3.2.113) 

Then, at second order in the background gauge fields, we can replace W^d^Wy with 
W'^.'V^Wy . If we write the space-time covariant derivative as the anti-commutator of the 
spinorial covariant derivatives, taking into account that W a is a covariantly chiral super- 
field and using the Bianchi identity Vl^ + VH^ = 0, we find that the expression we are 
considering can be written as 



w a x di & w; = iW a x V a V P W% + 0(B 3 ) 

= %W a x C afi V 2 Wl + 0(B 3 ) 

= iW a x V 2 W ya + 0(B i ) 

= iW«D 2 W ya + 0{B 3 ) . (3.2.114) 
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Hence, the integral expression (13 .2. 1 13[) can be put as 

J d 4 xd 4 yd 4 9W a x {d y a& W\f{x-y) = i j d 4 xd A yd 4 9 W a x {D 2 W a ) y f(x - y) + 0{B 3 ) 

(3.2.115) 

Therefore, we have a gauge invariant expression for the ^-contributions of the form 

2 r 

Y,^V° OP \r = 3g 2 C 2 A tr j d 4 xd 4 yd 2 9W a (x,9)W a (y,9)n[AI°] R (x-y) 

ki(x - y) 2 M 2 



i=l 



3g2 ° 2A ;tr I d 4 xd 4 yd 2 9 W a (x,6)W a (y,9)n'- 



2\3 



16(4tt 



The r contributions are also added up and renormalized as 



(x - yY 



0{B 3 ) . 
(3.2.116) 



i=3 



x 



1 ln(a; - y) 2 M 21 
32(4vr 2 ) 3 {x - yf 



-[AI°] R (x-y) 



+ 0(B 3 ). (3.2.117) 



As this expression is transverse we can rewrite it in terms of the background field strength. 
To do so, we have to use the following property: Let / be a generic function, then, to second 
order in the background gauge fields we have 



tr j d 4 xd 4 yd 4 9 T a «(x, 9)T^(y, 9) - 2C aP C.p) f(x - y) 

= -3tr J d 4 xd 4 yd 4 9 [D a B(x,9)] [D 2 D a B(y, 9)] Uf(x-y) + 0{B 3 ) 

= 3tr J d 4 xd 4 yd 2 9 W a (x, 9)W a (y, 9)Df(x - y) + 0{B 3 ) . (3.2.118) 

To prove this relation we have only to write the connection in terms of the background 
gauge field as T a a = [DaD a — (i/2)d a a) B + 0(B 2 ). As a consequence of having a transverse 
structure, the integral is then written as 



tr d 4 xd 4 yd 4 9 [D & D a B(x,9)} D^D^B(y,9) (d x a6l d x b ,-2C a pC &$ n)nf{x-y) 



tr d 4 xd 4 yd 4 9 [D a B(x,9)} [D 2 D^B(y,9)]C^[d x a& d^ $ -2C a pC & p)nf(x-y) 



(3.2.119) 
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where we have integrated by parts the superspace derivative D a . Then, as C^ a O a aOpp = 

6j>n, C^CcpC^ = -2C a(3 and W a = (l/i)D 2 D a B + 0(B 3 ) we straightforwardly obtain 
the relation (I3.2.118p . With this, we can express the space-time connection contribution 
(that is transverse as is required) in terms of the background field strength. Adding up 
all the results and taking into account the symmetry factor, we find the total two-loop 
renormalized contribution to the background gauge field self-energy to be 



5 EH 



2loop - trjd i xd 4 yd 2 eW a (x,e)W a (y,e)Tf 2loop (x-y), (3.2.120) 



with 



r(2 ) 2 i ooP( , 3g 2 C 2 A \ In 2 x 2 M 2 R + \ \nx 2 M] R (l - In x 2 M 2 ) + lnx 2 M 2 

B [X) ~ 64(4vr 2 )3 x 2 + '" ' 

(3.2.121) 

where Mm corresponds to the IR divergence and M to the UV one. As is expected, we 
have found a result that has both types of divergences. However, opposite to the case of 
dimensional regularization, we are able to clearly distinguish them. 



3.2.4 RG equation 

With our conventions, remembering that in the background field method the background 
field and the coupling constant renormalizations are related {Z g \/ Zg = 1), the anomalous 
dimension of the background field cancels. Hence, the RG equation for this fields is^| 





M- h 0(a)— + 7^(9)7- 

QM Og ^ yy! 0i 



r 



(3.2.122) 



C=o 



with containing the free action, the one- and two- loop contributions ( (13.2. 76[) and 
(13.2. 121 j) respectively) and the linear expansion in £ of the effective action in a generic 
gauge (1372788]) : 

In x 2 M 2 r 



1 3C A D lnx 2 M 2 

V W + 16(47r 2 ) 2D— x 2 



3g 2 c 2 



256(4vr 



2\3 



-□ 



16(4vr 

In 2 x 2 M 2 i j + 21nx 2 M 2 p (l 



2^2 



-□- 



L IR 



l IR\ 



lnx 2 M 2 ) +41nx 2 M 2 



+ ... 
(3.2.123) 



As in all the previously considered models, the gauge-running term 7^ is evaluated in 
appendix [C] with the one- loop RG equation for the quantum gauge field self-energy. We find 
there that this function has an expansion to second order in g 2 of the form 

75 = 



4(4tt 2 



(3.2.124) 



5 Note also that in the RG equation we have not included a term — Mnid/dMm because in the 
renormalization we have required that both IR and UV scales were independent, which implies that — 
M/M IR dM IR /dM = 0. 
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Hence, inserting (13.2.1231) in the RG equation, and using the result obtained for 7^, we 
find the following value for the one- and two-loop coefficients of the expansion of the beta 
function 



M = - K H 9 5 - (3-2.125) 



4 V8tt 2 J y 4 V8vr 2 

Let us remark again that as our supersymmetric conventions differ in a \/2 factor in the 
coupling constant wrt. the usual ones [M] (g = V^gsYivi), this results matches the stan- 
dard beta function expansion (3(g SY M) = -(3/2)[C a /(8tt *)]gly M ~ (3/2)[C A / (8vr 2 )] 2 c/f yM + 
0(9sym)- 

Discussion of the result 

There has been some controversy about the origin (UV or IR) of the higher-order perturba- 
tive contributions to the beta function in supersymmetric gauge theories. The "exact beta 
function" of N = 1 SYM was discovered by Novikov, Shifman, Vainshtein and Zakharov 
(NSVZ) in [8] (although the expression was first derived in [C32J ) , and it is of the form (in 
this discussion, we will use the usual coupling constant, thus implicitly we have g = Qsym) 

Mg) = ^_ . (3.2.126) 

HKyj 16tt 2 i_^m! v ; 

8tt 2 

In the NSVZ derivation of the function, instanton analysis was used, showing that the 
higher-loop corrections to the one-loop result were due an imbalance in the number of 
fermionic and bosonic zero modes. Thus, these correction had and IR origin. Afterwards, 
the two-loop coefficient of the expansion of the beta function was obtained using dimensional 
reduction [551 EE EH], matching (13.2. 126[) . However, as we have pointed out previously, in 
this calculation IR divergences appear and dimensional reduction regularizes both types of 
divergences with the same parameter [52] • I n [56J this is solved by subtracting the IR diver- 
gences by a R operation [16] , obtaining that the two- loop correction of the beta function has 
its origin in an specific operator of dimensional reduction that is not available in a renormal- 
ization procedure that stays in four dimensional, which seems to imply, as is pointed out in 
[55] . that no divergence should occur beyond one- loop. 

This situation is cleared if we distinguish between the running of the physical coupling 
constant (the constant that is used in perturbative calculations when obtaining the 1PI 
effective action), and the running of the coupling constant in the Wilson's effective action 
approach^- In the latter case, using different arguments as the holomorphic dependence on 
the complexified coupling constant or the fact that the relevant domain of the non-local 



6 This operator is written in terms of the background space-time connection T aa and the Kronecker delta 
functions in four (6^£) and n {Sj^) dimensions as tr f d 4 xd 4 9 r QQ T gb{5^P — 5j&), which, by means of a 
Bianchi identity can be put in terms of the classical action as —etr J d 4 xd 2 8 W a W a 

7 In this method initially we consider a theory defined with a cutoff scale M, and study how the different 
terms of the lagrangian flow when we integrate over momentum slices down to another scale M'. This flow 
implies that the coupling constants verify RG equations of the form Md/dMX = /3(A) 
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operators responsible of the higher-loop corrections to be virtual momenta of the order of 
external momenta (and therefore excluded by definition in the Wilson action) [HI |63] , we 
can conclude that the flow of the coupling constant is exhausted at one loop. However, in the 
case of the physical coupling constant, we have higher-order contributions that appear when 
we take the expectation value (in a external field) of the operators in the Wilson action, 
being the relevant IR pole related to an anomaly [63J. 

This IR origin of the higher-order corrections is questioned in P, |T0] . In these works, in 
a purely wilsonian framework, an NSVZ flow is obtained. The key idea is the differentiation 
of the flow of an "holomorphic" coupling constant and a "canonical" coupling constant. The 
first "holomorphic" coupling constant corresponds to a lagrangian which is normalized at a 
scale M as 



with \jg\ = l/g 2 + i0/(8ir 2 ) being the complexified coupling constant, whereas the "canon- 
ical" coupling constant corresponds to a lagrangian of the form 



Although the g^ running can be shown to be one-loop, when we consider the running of 
the coupling constant g c we found that obeys a NSVZ flow. The reason is that in order to 
maintain "canonical" normalization in the lagrangian, we are forced to perform a rescaling of 
the fields, being this rescaling anomalous and the origin of the higher-order terms. In these 
works it is also argued that the flow of the 1PI coupling constant is closely related to this 
"canonical" coupling constant flow, being this confirmed in [M], where it is found that the 
first two coefficients of the 1PI beta function coincide (in any mass-independent scheme) with 
the first two coefficients of the "canonical" wilsonian beta function [15] . As the construction 
we have just described is made a la Wilson, it is claimed in [9j [10] that it only depends 
on the UV properties of the theory. However, in [SS], this interpretation is criticized, as 
it is pointed out that the IR degrees of freedom must be included in the derivation of the 
anomaly, if we want to maintain low-energy physics unchanged under rescaling [65, 66J. 
We proceed to the discussion of our result. We have taken the following steps 

1. We first renormalize the one-loop UV sub divergences with a scale M. 

2. If we have had an overall UV divergence, we could have renormalized it with a different 
scale, say M 1 . By power counting only IT diagrams could have overall UV divergences; 
nevertheless, as the traceless parts multiplying IT are finite, they can only depend on 
the one-loop scale M, which implies that we can not have any M' dependence, as gauge 
invariance imposes transversality in the IT expressions. 

3. Although we could have when taking the derivative wrt. the UV scale in the RG 
equation a non-local dependence in M (see (13.2.1121) and (I3.2.117P ). after integration 
over half of the supercoordinates these contributions become local. 




(3.2.127) 




(3.2.128) 
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4. The one-loop scale is cancelled by the two-loop coefficient of the beta function in the 
RG equation, playing the off-shell IR scale a passive role, as it is exactly cancelled in 
this equation. 

Hence, we conclude that the scale associated to the one-loop renormalization of the quan- 
tum superfield is the one that gives rise to the two-loop coefficient of the beta function. The 
fact that no overall UV scale appears, is directly related to the conclusion obtained in [55] 
that in a four- dimensional regularization method there are no superficial divergences. How- 
ever, as we have seen, this not implies that the two-loop coefficient of the beta function 
cancels. The mechanism presented here agrees with previous calculations in which the cor- 
rections to the one- loop result arise from a one- loop anomaly [63j [9j [TOj EH [68] . In our case, 
the anomaly is to be associated with the external loop, and is responsible of the promotion 
of the M dependence into a non- vanishing non-local structure that eventually generates the 
two-loop coefficient of the beta function. So, we have found that the two-loop coefficient of 
the beta function arises from a one-loop UV scale which survives at two loops when IR effects 
are included. 
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Chapter 4 
Conclusions 



We conclude that differential renormalization is a useful method we can employ when renor- 
malizing a gauge theory. Although the strengths of the method are well-known (e.g., gauge 
invariance is not broken, and we stay all the time in four dimensions, which is crutial when 
studying supersymmetric theories), in its original formulation it has, at least for us, one 
important drawback: the necessity of imposing explicitly the Ward identities in each cal- 
culation with gauge theories. However, this point was solved for the one-loop case by the 
introduction of Constrained Differential Renormalization. 

We have shown that we can made fruitful use of the one-loop CDR results in two-loop 
calculations, due to the fact that CDR fixes all the ambiguities related to the logarithms of 
the scales at the two-loop order. We have distinguished two cases: in the first one we have 
diagrams where the divergences are "nested" . This implies that we can directly apply CDR 
to the "inner" divergence, which straightforwardly fixes all the coefficients of the logarithms 
of the scales in the total two-loop expression. The second case corresponds to diagrams with 
overlapping divergences. To deal with them we have obtained a list of renormalized two-loop 
integrals where in each calculation one-loop CDR rules have been maintained in every step. 
Although the problem is not solved for the general the list is restricted to two-point 

functions with at most four derivatives acting on the propagators and two free indices, these 
integrals are the expression we have to deal with when we use the background field method. 

We have discussed also the application of DiffR to IR divergent expressions. Although 
the renormalization procedure in momentum space resembles the usual one, we have one 
subtle point: the co-existence in the same expression of UV and IR divergences. In this case, 
as both renormalizations should decouple (UV and IR divergences are local in position and 
momentum space respectively, so that Bogoliubov's UV R and IR R operations commute), 
the scales related to each type of divergence must be independent. In order to guarantee 
this, we have found that we have to modify the usual renormalization relations by means of 
an adjustment of the local terms involving both scales. 

In this work we have re-obtained two results that where previously derived with usual 
differential renormalization and Ward identities: the two-loop beta function of QED and its 
supersymmetric extension, SuperQED. In both cases, we have shown that using one-loop 
CDR simplifies the calculations, as expressions that vanish by symmetry automatically can- 
cel, and we do not have to relate the different scales via Ward identities in the renormalized 
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results. 

However, we have not only re-obtained previous results, but we have also performed 
two of the relevant calculations that were pending with differential renormalization: the 
two-loop renormalization of Yang-Mills and N — 1 Super Yang-Mills. With the first one we 
have found no difficult when performing the calculation, as our use of one-loop CDR results 
allow us to perform the renormalization with the same ease as with standard dimensional 
methods. With Super Yang-Mills theory, the use of differential renormalization has one clear 
advantage over dimensional reduction (which is the usual regularization method employed 
with supersymmetric theories): in this case we have both UV and IR divergences, and with 
dimensional reduction they become mixed (both are renormalized with the same infinitesimal 
dimensional parameter), being necessary to subtract the IR contribution in the final result. 
Differential renormalization, however, clearly distinguishes between both divergences as they 
are renormalized with different scales. This feature allows us to give new insight on the 
origin (UV or IR) of the higher-loop contributions to the beta function, which has been a 
controversial point. We have found that higher order corrections to the beta function come 
from the one-loop UV scale, which survives in the higher-loop expression by the presence of 
the IR divergences. 

Finally, among the different open problems that we have, it is clear that the principal one 
is the extension to higher-loop order of CDR. To achieve this, we have to obtain a complete 
set of rules that fix the local ambiguity of the higher-order expressions as it has been derived 
for the one-loop case. We think that the results we have found are a step ahead in this 
direction, as the complete CDR renormalized expressions must coincide, at least in parts 
corresponding to the logarithms of the scales, with the renormalized results that we have 
presented here. 



Appendix A 

Conventions for supersymmetric 
calculations 

In this section we will briefly review the most relevant results and conventions of supersym- 
metry and superspace that are used in this work. Although this topic is covered with great 
detail in various references (e.g. [SH EH 171]). we will follow closely [38], where the reader 
can found a complete treatment of the subject. 

A.l Notation 

Setting up the notation, we express the vectors (representations (5,5) of the Lorentz group) 
with two pairs of two- valued spinor indices, one dotted and another undotted V aa . To relate 
a vector in an arbitrary basis with index a we use the Pauli matrices: 

Fields: V a « = ^af'V^ V h - = ^o~ 6 y adl 

Derivatives: d aa = a\ & db d k = |cr£ d 9 ad 
Coordinates: x aa = \a^ a x- x- = cr^x aa 
Pauli matrices satisfy 

<^ Q = 2^ , a b -^ = 26 a %? (A. 1.1) 

With this conventions, the Super Yang-Mills coupling constant (g) that we use is related to 
the usual one {qsym) by g = \f2g S YM [31] ■ 

A graded commutator [f^, Q B } = Q A Q B — (—) ab Q b Qa is defined as the anticommutator 
{Qa, &>b} when both Qa and Q B are fermionic operators and the commutator [Qa, &b] 
otherwise. Symmetrization and antisymmetrization (sum over all permutation of indices 
with the corresponding sign in the case of antisymmetrization) are indicated by ( ) and [ 
respectively. Indices between vertical lines | | are not taken into account in the previous 
operations. We also define a graded antisymmetrization [ ) so that [Qa,^b} = ^[a^b)- 
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For raising and lowering spinor indexes, we use the matrices C a p that have the following 
properties 

a/3 - ~ % 6 P] 

C aP = cfe C a p = C^ (A. 1.2) 

With these matrices, if ip a denotes a spinor we define 

ip 2 = \c aP r^ = \r^ a (a.i.3) 

i[>a = ^Cp a , r = C a ^ p (A.1.4) 

Similar relations hold for the hermitian conjugate ip a . In the case of a vector V aa , we define 
the square to be 

V 2 = \v™V a «. (A.1.5) 

A. 2 Supersymmetric algebra 

Coleman and Mandula [5] showed that any group of bosonic symmetries of the S-matrix is 
the direct product of the Poincare group with an internal symmetry group. This implies that 
the commutator of the bosonic generators of the internal symmetry group and the generators 
of the Poincare group (J a p, J a biP a p) van ish. In [6], Haag, Lopuszanski and Sohnius avoided 
this no-go theorem by allowing fermionic symmetry generators Q aa (where a = 1, . . . , N is 
an isospin index). They found the most general super- Poincare algebra to be 



^ Qaa 


,Q b p} 


— () b P ■ 

— °a r aP 




{Qaa ; 


Qbp) 


— CapZab 




[Qaa 


' p pp] 


= [Paa , Ppp] 


— [-Jap > Qn] — 


[Jap 


, Qcy\ 


— 7^~i{aQcP) 




[Ja/3 


, P 77 ] 


= -C 7 (a-P/3) 7 




[Ja/3 


, J^] 


= 

2 ( a P) 




[Jap 


' Ja$] 


— [Zab > Z cc [] = 


[Z ab , Z cd ] =0 



(A.2.6) 

where Z a \, are |iV(jV— 1) complex central charges. The N = 1 case is called simple supersym- 
metry, whereas N > 1 is called extended supersymmetry. We consider only the N = 1 case. 
For theories satisfying super-Poincare algebra some remarkable properties can be derived 
[38]: 



Equality of bosonic and fermionic freedom degrees. 

Energy is positive, and in particular, vacuum energy can be shown to vanish. 



A.3. SUPERSPACE AND SUPERFIELDS 
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A. 3 Superspace and superfields 

A compact technique for working with supersymmetric theories is what is called superspace 
|72j . By means of anticommuting parameters we can integrate the super- Poincare algebra 
and obtain a group, the super- Poincare group|38, 69J. As usual spacetime can be defined as 
the coset space Poincare group/Lorentz group, superspace is defined to be the coset space 
super-Poincare group/Lorentz group. Hence, superspace is a space spanned by the usual real 
commuting space time coordinates and new anticommuting coordinates z A = (x aa , 9 a , 9 a ). 
Supersymmetry generators are realized as coordinate transformations in superspace [38J . As 
the usual generators of the Poincare algebra can be represented by differential operators, Q a 
and Qa are found to have the following expression [38] 

Qa i i^da ^aa 

Q a = i(d a - l -6 a d aa ^ (A.3.7) 

It has to be noted that neither usual coordinate derivative d aa nor fermionic coordinate 
derivatives d a , d a are invariant under supertranslations (those that are generated by Q a , 
Qa). However, supersymmetric covariant derivatives that are invariant under these trans- 
formations can be defined as [38] : 

% — • 

D a = a + -0 a d aa 

D a = d a + l -e a daa (A.3.8) 

We list here some relevant algebraic relations of these derivatives 

^D a , -DqJ id a a 

[D a , Dp} = {D a ,Dp}=0 

[D a , D 2 ] = id a «D a 

{D 2 ,D 2 } = a + D a D 2 D a 

D 2 D 2 D 2 = UD 2 

D 2 9 2 = -1 . (A.3.9) 

Also, with these derivatives we can define two projection operators as 

III = -D a D 2 D a /0 

2 ' 

n = (D 2 D 2 + D 2 D 2 )/U , (A.3.10) 

which verify that IIi + IIo = 1. 
A. 3.1 Superfields 

Superfields are defined to be multispinor functions over the superspace $ = $(a;, 9, 9). It is 
clear that Taylor-expansion of these functions in terms of 9 a and 9 a breaks off at order 9 2 9 2 , 
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due to the anticommuting nature of these parameters. The different terms of the expansion 
are called the component fields. We can impose constraints into a superfield and reduce 
the number of independent component fields. The most relevant ones for our work are the 
following two: 

Chiral superfields With the aid of the covariant derivative ()A.3.8j) we can impose the 
constraint 

D & $ = 0. (A.3.11) 

If <3> verifies the previous relation it is called a chiral superfield. This constraint implies 
that this superfield has the following component expansion 



$ = A + d a i) a -6 2 F + l -6 a 6*d a aA+ l -d 2 6*d a ai) a + ^6 2 6 2 UA . (A.3.12) 

Real superfields A superfield V is called a real superfield if it verifies the constraint 
V = V + . With this constraint, the component expansion is found to be [38] 



V = C + 9 a \ a + O^Xa ~ 2 M — 6 2 M + 9 a 9 & A a & — 9 2 9 a \ 



a 



-e A 9 a Xa + 9 A 9 A D . (A.3.13) 

A. 3. 2 Superspace integration and superfunctional derivation 

In order to obtain supersymmetric invariant actions, we have to define the integration over 
the anticommuting coordinates. The basic properties of these integrals are 



d9 9 = 1 

d9 1 = (A.3.14) 

With this definition, the delta function of the anticommuting variables is found to be 

S(6 - 9') = (9 - 9') . (A.3.15) 

These properties imply that the integration over 9 is identical to differentiation, or, which 
is the same, that inside a d A x integration we have D a = J d9 a . Also, since supersymmetric 
variations are total derivatives, if we consider \I/ a general superfield and $ a chiral superfield, 
the following quantities are supersymmetric invariants 



S* = J d 4 xd 4 9 * 

S* = [ d 4 xd 2 9 $ (A.3.16) 



As it will be necessary when studying perturbation theory in superspace, we have to 
consider the superfield extension of the functional derivative. It is found [38] for a general 
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superfield \I> and a chiral superfield $ that the superfunctional derivation is 



S^(z') 



5\ z - z') = 5\x - x')5\9 - 9') 
D 2 5\z-z') (A.3.17) 



Finally, this integral definition implies that we have the following integration by part 
rules for the superspace derivatives D a 

J d 8 z A(D a B)C = J d s z \-(D a A)BC - AB(D a C)} 

J d 8 z A{D a E p )C = J d 8 z \-{D a A)E p C + AE p (D a C)] 

J d 8 z Ep(D a A)C = J d 8 z [{D a Ep)AC - EpA(D a C)] 

J d 8 zE p {D a F 1 )C = J d 8 z [{D a E p )F 1 C + EpF^D^)] , (A.3.18) 

where A, B and C are bosonic (conmmuting) superfields and Ep, F 1 ferminic (anticommut- 
ing) superfields. 

A. 4 Superspace formulation of supersymmetric gauge 
theories 

Among the different models we can consider in superspace (Wess-Zumino, nonlinear a— models, 
etc.), we will restrict ourselves to gauge theories. To deal with them there are two different 
approaches that we detail here, as both will be useful when studying the supersymmetric 
extension of the background field method. 

A. 4.1 Chiral representation. Prepotentials 

Starting with the linear case, by studying the field content of the N — 1 vector multiplet [38J, 
we can obtain that the corresponding irreducible off-shell field strength is a chiral superfield 
W a that satisfies 

D a W a = -D^Wa . (A.4.19) 
Hence, W a can be expressed in terms of an unconstrained real scalar superfield V by 

W a = iD 2 D a V 

W A = -iD 2 D a V (A.4.20) 

This definition is clearly invariant under gauge transformations with a chiral parameter A 
of the form of 

V = V + i(A-A) . (A.4.21) 
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This linear study has to be generalized to the non-abelian case. To do so, we consider a 
multiplet of chiral scalars fields transforming according to some representation of a group with 
generators T A - These fields are postulated to transform with chiral parameters A = A a Ta 
as 

$' = e l9A $ , (A.4.22) 

with g a coupling constant. Considering an antichiral field $, we found that it transforms 
with the complex conjugate representation with antichiral parameter A. To obtain an in- 
variant expression we introduce a multiplet of real superfields V A that transforms as 

(e 9V )' = e l9A e 9V e~ i9A , (A.4.23) 

which implies that $e 9y $ is invariant under gauge transformations. With the prepotential 
V we can construct derivatives gauge covariant with respect to A transformations as (38] 

v A = D A -ir A = {v a ,v & ,v a6l ) 

= (e- 9V D a e 9V , D a , -i { V a , V d }) , (A.4.24) 

which satisfy the requirement 

(V A $)' = e l9A {V A <5>) V A = e l9A V A e~ i9A (A.4.25) 

These derivatives are called gauge chiral representation covariant derivatives. Their conju- 
gates V, that we call gauge antichiral representation covariant derivatives, are covariant with 
respect to A transformations 

V A = (D a ,e? v D & e-* v ,-i{V a , V d }) . (A.4.26) 

Both representations are related by a nonunitary similarity transformation 

V A = e 9V V A e- 9V • (A.4.27) 

Field strengths defined by commutation of the covariant derivatives can be expressed in 
terms of the following fields denoted as W a and W a [38] 

W a = iD 2 (e~ 9V D a e 9V ) 

W a = e- gV W a e gY = e- 9V (-W a ) + e 9V (A.4.28) 
that satisfy Bianchi identities of the form 

V a W a = -V & W a . (A.4.29) 
With these fields we can construct a gauge invariant action as 

S = i-tr [ d A xd 2 9 W 2 
9 2 J 

= — Ktr l d A xd A 9 (e~ 9V D a e gV )D 2 (e- 9V D a e 9V ) . (A.4.30) 
2^ 2 J 
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A. 4. 2 Vector representation. Covariant approach 

In this approach we start defining covariant derivatives and, by means of covariant con- 
straints, express all quantities in terms of a single irreducible representation of supersymme- 
try. 

For a Lie algebra with generators T A we covariantize the derivatives with the introduction 
of connection fields Ta = F A T B as 

V A = D A - iT A ■ (A.4.31) 

Under gauge transformations these derivatives are postulated to transform with a real su- 
perfield K = K A T A as 

V A = e i9K W A e~ i9K . (A.4.32) 

By commutation, field strengths Fab are defined in terms of the connections and the flat 
superspace torsion T A ^ (where T - 77 = i5^5 - 7 is the only nonzero component) as 

OLjJ jj 

[Va,V b } = T A gv c -iF A B 

F AB = D [A r B} -i[r A ,r B }-T A grc . (A.4.33) 

Over these field strengths we can impose different constraints. 
Conventional constraints 

Due to the fact that one can always add a covariant term to the connection without changing 
the transformation of the covariant derivative, we can impose the constraint 

F a6t = , (A.4.34) 

as without imposing it we can define new connections T' A = (T a , F a & — iF a a) that iden- 
tically satisfy the constraints. Hence, the covariant derivatives take the following form 

V A = (V Q , V d , -i { V a , V A }) . (A.4.35) 
Representation-preserving constraints 

Let us define a covariantly chiral superfield $, ie. a superfield that verifies 

V d $ = , $' = e igK <$> 

V a $ = , $' = §e- igK (A.4.36) 
Hence, consistency implies that we have to impose the constraint 

F aP = F. = O, (A.4.37) 
as with the previously defined covariantly chiral superfield we have that 

= {V & ,V $ }* = -iF &$ *. (A.4.38) 
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This constraint is solved with a complex superfield f2 = Vl a Ta which allow us to write the 
covariant derivatives as 

V a = e~ 9n D a e 9n 

V A = e? Si D & e- aa (A.4.39) 
With this superfield we have two types of gauge transformations 

• K gauge transformations 

(e 9Q )' = e 9n e~ l9K . (A.4.40) 

• If we consider A to be an usual chiral superfield -D^A = 0, the derivatives defined in 
flA.4.39j) are invariant under 



(eP Q )' = e i9A e 9Q . (A.4.41) 

From the K-invariant hermitian part of Q we can define a real superfield V as 

e v = e gn e gn _ (A.4.42) 

We can also use the Q superfield to write all the quantities in a gauge chiral representation, 
where everything transforms only under A transformations [3"8] : 

Vqa = e~ 9n V A e 9n 
$ = e ~ sn $ , (A.4.43) 

where Voa are the chiral representation derivatives (IA.4.241) . $ is a covariantly chiral super- 
field and $o a chiral superfield. 

Field content. Bianchi identities 

The field content of the theory can be obtained through the Bianchi identities that are 
derived from the Jacobi identities satisfied by the covariant derivatives. With the aid of 
these identities and the different constraints imposed on the derivatives, all of the field 
strengths can be expressed in terms of a spinor superfield W a |38J that is defined as 

[Va,^] = -iCfcWp. (A.4.44) 

This superfield can be shown to be covariantly chiral (V gW a = 0) and to satisfy the following 
identity 

V a W a + V & W & = . (A.4.45) 
Finally, with these deriatives we can construct the gauge Lagrangian as 

trW 2 = -hr ([V* , {V A , V Q }]) 2 . (A.4.46) 
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A. 5 Supergraphs 

Although supersymmetric theories can be quantized at the component level with conventional 
methods, the use of superfields simplifies all the calculations. Along with the compact 
notation and automatic cancellation of graphs related by supersymmetry, in a superfield 
formalism supersymmetry is manifest. We will detail the extension of the usual functional 
methods to superspace [38], |6T] . 

Let \I/ be a generic superfield, S($f) the action and J a source of the same type (chiral,etc.) 
as The generating functional for Green functions is 



Z[J] = J [^]e 5 W+/ J * . (A.5.47) 

For connected Green functions, the generating functional is 

W[J] = \nZ[J}. (A.5.48) 
Finally, the generating functional of 1PI graphs (the effective action) is 

r[£] = W[J(if)] - J J(*)# , (A.5.49) 

where ^ is the expectation value of \I> in the presence of the source: 

SW 

* = -i--. (A.5.50) 

Let us consider two examples: a real scalar superfield (as in gauge theories) and a chiral 
superfield. With the first one, the superspace partition function can be written as 

Z[J] = j\dV]exp^j\- l -VUV + C mt {V) + JV]^ 

= exp!^J £int{jj)\ e %P |^ j J 7f J } ' (A.5.51) 
whereas in the case of a massless chiral superfield, we have 

Z\j,j] = J exp | j d 8 z [$$ + £ int ($, $)] J d G z 3 <$> + J d 6 fj$\ 

= expl^j £ int {^,-jL)jexp^-J j'^jj . (A.5.52) 

From the expansion of these expressions we can derive the propagators, vertices and sym- 
metry factors. 
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Superspace Feynman rules 

• Propagators: We present here the propagators for massless chiral(antichiral) and real 
superfields. These propagators are those that we need in our work. A detailed list with 
more propagators can be found in [38j or [61J. 

VV : -P{zi - z 2 ) = -A( Xl - x 2 )8\9 l - 9 2 ) 

<5$: P(z 1 -z 2 )=A{x 1 -x 2 )5 4 (9 1 - 9 2 ) (A.5.53) 

• Vertices: For each chiral(antichiral) line there is a D 2 (D 2 ) factor acting on the propa- 
gator. If all the lines are purely chiral or antichiral, one of the factors is omitted. 

• Apart of the usual space time integrals (or momentum space integrals), there is a d 4 9 
integration at each vertex. 

• When computing 1PI graphs (in order to obtain the effective action), for each external 
line we multiply by the corresponding superfield. In the case of a chiral (or antichiral) 
line, no D 2 (D 2 ) factor appears. 

• Some diagrams have symmetry factors. 

A superspace diagram which corresponds to a contribution to the effective action is an 
expression formed by some external fields, supercoordinate integrals J d A Xid A 9i, propagators 
(Pij) and superspace covariant derivatives acting on them [38|, [6Tj. The propagators are of 
the form Py = Aij5(9i — 9j) = Ay#y with Ay the usual spacetime propagator and <$y the 
5-function on the anticommuting coordinates. The covariant derivatives can be integrated 

by parts, obey the Leibnitz rule and a "transfer" rule of the form of 5y Dj— —Di5ij. So, we 
can choose a propagator that links two vertices and remove all the .D's from its 5-function. 
Then, if we have other propagators that link these two vertices, we can apply the following 
properties 

SijD^DfSij = 5 lJ D 2 D 2 5 ll = ^Df^Aa^y = (A.5.54) 

Now, with the free superspace 5-function, we can contract the propagator between the two 
vertices to a point in #-space. As this procedure can be repeated for other pair of vertices, 
we conclude that we can write the effective action as 

r = [ dAx i ■ ■ ■ d 4 %nd 4 9 T(x 1 , . . . , x n )$(xl, 9)... V(xi, 9)... (A.5.55) 

This expression has one important consequence. We have found that in a perturbative 
calculation a contribution to the effective action with a purely chiral (antichiral) integral 
d 2 9{d 2 9) never gets generated. Hence, if the original action had purely chiral terms (as 
mass terms $ 2 or cubic interactions <3? 3 ) they can not be modified with radiative corrections 
[731 [7i] . This is called the no-renormalization theorem for chiral superfields. 
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A. 5.1 Quantization of supersymmetric gauge theories 

As in usual gauge theories, when we quantize a supersymmetric gauge theory with functional 
methods we have to fix the gauge, as in the path integral we do not have to integrate over the 
physically equivalent gauge field configurations related by gauge transformations |47J. Hence, 
we will present here the generalization of the usual gauge fixing procedure to superspace. 
We start with the previously defined SUSY Yang-Mills [69] action 

S = — l -tr [ d A xd 4 9 (e- 9V D a e 9V )D 2 (e~ 9V D a e 9V ) . (A.5.56) 

2r J 

Then, with A (the chiral parameter of the gauge transformation), an arbitrary function / 
and a gauge-variant function F such that F = f for some value of A, we define a functional 
determinant as 



A(V) = l[dAdA]5[F(V,A,A)-f]5[F(V,A,A)-f] 

\dAdAdA'dA']J de ' zA i^ A+ ^ A )^ de ' 2A, (^ A+ & , (A.5.57) 



where the variational derivatives are evaluated at A = A = 0. 

Then, we introduce inside the partition function the unity in the form of the functional 
determinant and its inverse. With a change of variables that is a gauge transformation we 
set the A integral as a constant infinite factor reabsorbed into the normalization [38]. Also, 
in order to get a result independent of / and /, we average with a gaussian weighting factor 
of the form J[dfdf] exp(—-tr J d s zff). Hence, the partition function is written as 

Z = J [dV] (A(V))- 1 exp js - ~*r J d 8 z Ff\ . (A.5.58) 

Using as gauge fixing function F = D 2 V and replacing the parameters A, A' of (1A.5.57I) by 
anticommuting chiral ghost fields c, d (superfield extension of Faddeev-Popov ghosts [32J) 
we find [38J 

Z = J [dVdcdcdcdc}e So+s ° F+SFP , (A.5.59) 



where 



S GF = --tr I d 8 z (D 2 V)(D 2 V) 
a J 

S FP = tr j d 4 xd 4 9 {d + c')Li gV 



[c + c) + cothLi gV (c — c) 



L X Y= [X, Y] 

(A.5.60) 
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Appendix B 

Background (super)field method 



In order to quantize a gauge theory with functional methods, a gauge fixing procedure has to 
be used, so as no gauge field configurations related to a given one by gauge transformations 
are taken into account in the path integral (all of them correspond to the same physical state). 
However, as a result of this procedure, explicit gauge invariance is lost. The background field 
method was developed (see references from [I] to [EH]) to allow us to fix a gauge without 
losing explicit gauge invariance. Although originally the method was developed for the one- 
loop case, it was soon extended to include multi-loop effects [I21 [HU [851 ESI I8"7| . The basic 
idea is the splitting of the gauge field in two parts: the quantum field, which is the variable 
of integration in the functional integral, and the background field. Thus, we are allowed to 
fix the gauge for the quantum field whereas we can maintain explicit gauge invariance in 
the background one. We will first discuss the non-supersymmetric case and then obtain the 
generalization of the method to superspace [Ml EE] • 



B.l Yang- Mills theory 

We will use the conventions of [88] which are those detailed in section 13.1.11 We begin 
defining the splitting of the gauge field in two parts as 

+ (B.l.l) 

where A ^ is the quantum field and is the background field. With this splitting, let us 
consider the functional Z[B] = J[dA]e- So( - A+B ^ where S Q (A) = 1/4 / d 4 xF° u (A)F^(A) is the 
usual Yang-Mills action. After the usual gauge fixing procedure [15] this functional becomes 
(c, c are Faddeev-Popov ghost fields) 

_ 5F(B) 



Z[B] = J [dAdcdc}exp ^-S (A + B) - -Ur j d 4 x F(A, B) 2 + tr J d 4 x c 



Notice now that Sq(A + B) is invariant under two types of transformations 



c 

(B.1.2) 
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1. Quantum 



8B° = 

5AI = ^[d^w a + gr bc Blw c ]+r bc Alw c 

= ^(D^) a + f abc Alw c (B.1.3) 



where is the background covariant derivative. 
2. Background 

= -d,w a + f abc B b w c 
g 

g A a = ftc^c (B14) 

Our aim is to fix the quantum gauge invariance and at the same time maintain the back- 
ground gauge invariance. Thus, the gauge fixing function has to transform covariantly with 
respect to background gauge transformations. Hence, we choose as gauge fixing function 
F a = (D^A^) 0- , which implies that Z[.B] becomes 

Z[B] = J[dAdcdc}e- s ^ 

[dAdcdc}exp ^-S {A + B) - ^-tr J d 4 x (D M A M ) 2 + tr J d A x c[D M X> 



(B.1.5) 

with (V^w) a = d^w a + gf abc (A b ^ + B^)w c . 

As can be seen, Zi[B] is manifestly invariant under background gauge transformations. In 
order to make the connection with the usual functionals (Z,W = In Z and T[A] — J J A — W 
with A = 5W/5J), we define another functional like Z[£?] but with the quantum field coupled 
to a source 

Z[J,B] = J[dAdcdc\e^ A ^ + ^ d4xJa ^ . (B.1.6) 

We remark again that by construction, this is explicitly invariant with respect to background 
gauge transformations. Starting with Z[J, B] we can define analogous functionals as the usual 
ones like W = hxZ and t[A,B] = J J A - W[J,B], with A = 5W/5J. If we perform the 
change of variables in the partition function A a ^ — > A a ^ — B^ is straightforward to arrive to 

m 

Z[J,B] = e- lJB Z[J,B} 1 (B.1.7) 

where Z[J, B] is the usual partition function with the gauge fixing and ghost terms evaluated 
in an unusual but nevertheless valid gauge that depends on the background gauge field. So, 
we have for the other functionals 



W[J,B) = - J JB + W[J,B] (B.1.8) 
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and, with A = 5W/5J being the usual classical field 

A = -B + A 



f[A,B) = J J(A + B)-W[J,B] 
= J JA-W[J,B] 



= T[A, B] (B.1.9) 

We have found is that T and the usual effective action T are related by 

f[A,B] = T[A + B,B\. (B.1.10) 

If we restrict ourselves to diagrams with no external A we have a relevant identity: r[0, B\ = 
T[B]. This implies that the usual effective action can be obtained through the evaluation of 
r[0,-B]. This quantity is computed by summing all 1PI diagrams with B fields on external 
legs (A = implies that no A field propagators appears on external lines) and A fields inside 
loops (as the functional integral is only evaluated on A fields). 

One of the consequences of the background field method is that the renormalization of 

1 /2 

the gauge parameter [g = Z g g) and the background field (B = B) are related. As we 
have explicit background gauge invariance, the infinites appearing in r[0, B] must take the 
form of a divergent constant times {F^ u ) 2 . At the same time, is renormalized as 

(F;j = zW[drBZ-drB* + gZ g zWf*°BlB°] . (B.l.ll) 

Hence, in order to get explicit background gauge invariance, the following relation must hold 

Z g = Z B 1/2 . (B.1.12) 



B.2 Super Yang-Mills theory 

In this section we apply the conventions for superspace discussed in appendix |A] (those of 
reference [38], which we also follow in this section). As the gauge transformation is non-linear 
in the supersymmetric Yang-Mills theory, a linear splitting in the gauge field is unsuitable 
[3"8"l loT] . In order to define a splitting, we will re-examine the Yang-Mills case from a different 
point of view. The Yang-Mills action is invariant under local transformations of the form 
8A a ^ = 1/g d^w a + f abc A b ^w c . If the transformation is global we still have invariance, with 
the gauge field transforming as a matter field 5A a ^ = f abc A b ^w c . Then, considering again a 
local w, we can gauge the global transformation using the background field to covariantize 
the derivatives. This covariantization is of the form 

{D^w) a = {d^w) a + gf abc A b ^w c - {V ,w) a + g f abc A^ 

= d,w a + gf abc (Al + B b Jw c . (B.2.13) 

Hence, we have a linear splitting because the gauge field is linear in the covariant derivative. 
The procedure for the supersymmetric case is completely analogous [38]. At the end, we 
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have to covariantize the derivatives with the background field, which implies that we have 
to replace Da — > with a background covariant derivative. However, we have to take 
also into account that the covariant derivatives in a supersymmetric gauge theory can be 
formulated in two representations: 

• Chiral representation: This is more suitable for quantization. Hence, in the background 
field method, we work with a quantum V in a chiral representation. 

• Vector representation: We do not have to quantize the background fields. Hence, vector 
representation is useful for these fields as background covariance will be manifest. In 
fact, we will show that we can work with the background covariant derivatives without 
introducing explicitly the background prepotentials. 

So, we define the supersymmetric splitting writing the covariant derivatives in a quantum 
chiral but background vector representation: 

V Q = e~ 9V V a e 9V 

V aA = -i{V a ,V & } (B.2.14) 

with g the coupling constant. 

If we go to a background chiral representation (as can be seen in IA.4.21 this is achieved 
by multiplying with e~ 9n (. . -)e 9Si , where f2 is the background prepotential), we can straight- 
forwardly see that this splitting is equivalent to [38] 

e gV _^ e 9^ e9 v e gn _ (B.2.15) 

The split derivatives Va transform covariantly under two sets of transformations: 
1. Quantum: 

e 9V _^ e igA e gV e -igA ^ V a A = VqA = 

V A -> V4 (B.2.16) 
Which implies that Va transforms as 

V A -> e i9A V A e~ i9A . (B.2.17) 



2. Background: 



Which implies 



e gV _^ e tgK e gV e -igK ? K = K 

V A -> e l9K V A e- i9K (B.2.18) 



V A -> e l9K V A e- igK • (B.2.19) 
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Although Va has different transformations under background and quantum transformations, 
is not difficult to show that the transformation of the unsplit gauge field is the same in both 
cases [38] . 

If we study the background splitting of an abelian theory (SuperQED), the situation is 
simpler. In concrete, we have a linear quantum-background splitting of the form 

V^V + B, (B.2.20) 

where V and B are the quantum and background gauge fields respectively. This can be seen 
from the general supersymmetric quantum-background splitting expressed in terms of the 
background prepotential fl HB.2. 15[) 

V n V Ct 

e — > e e e 

B = fl + fl. (B.2.21) 
With this splitting the two sets of transformations are 
1. Quantum: 

V -> V + i(A-A) 

B -> B . fB.2.22) 



2. Background: 

V -> V 

B -> B + i(A-A). (B.2.23) 

Let us now consider the background field quantization. We start defining, as in the Yang- 
Mills case, a partition function Z with the gauge field split. After the gauge fixing procedure 

— 2 

with a background covariantly chiral gauge fixing function of the form of F = V V (which 
implies that the Faddeev- Popov ghosts are also background covariantly chiral) this functional 
becomes [38] 

Z = J [dVdcdc'dcdc'} 5(VV - f)5(V 2 V - f) e So+SFP . (B.2.24) 

In this case, due to the fact that we are dealing with constrained background chiral superfields 
rather than usual chiral superfields, when we gauge-average we have to consider a more 
sophisticated function. If we average with a factor like exp J fMf, with M an operator, in 
order to normalize we have to divide by detM. Hence, if M is a function of the background 
field we must average with an expression of the form 



J [dfdb]e fMf e' 



(B.2.25) 
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with b a field of opposite statistics to / that is called Nielsen-Kallosh ghost [58]. As this 
field only interacts with the background field and enters quadratically in the action, it only 
contributes at the one-loop level. In our case we gauge-average with a factor as 

J [dfdfdbdb]e- J" dSz [If+lb] , (B.2.26) 

where it is clear that b,b are background covariantly chiral ghost fields. All of this implies 
that the partition function can be written as 



[dVdcdc'dcd^dbdb}e b ^f 
Seff = S + Sgf + S fp + Jbb. (B.2.27) 

As in the usual Yang- Mills case, we can relate the background field functional with the usual 
effective action (in a special background gauge) once we set the sources to zero and consider 
diagrams with external background lines and internal quantum propagators [38J. 



B.2.1 Covariant Feynman rules 

After the supersymmetric quantum-background splitting, two approaches can be followed to 
perform the calculations. The first one is to use explicitly the background connections of the 
covariant derivatives in the calculations (V Q = D a + F a and V ft „ = d a a + ^aa)- So, in this 
situation we can apply usual D-algebra [55]. On the other approach, we do not extract the 
spinor connection of the covariant derivatives. Therefore, instead of using usual D-algebra, 
we apply the covariant D-algebra defined for these derivatives. Supergraphs obtained in this 
way give contributions that are only functions of the space-time connection T a a or the field 
strength W a . Not only these diagrams are simpler and fewer in number than those of the 
first procedure, but they are more convergent, as we do not have contributions with T a , 
which is of lower dimension that W a and r^. In this section we will detail this second 
approach. 

Let us consider a quantum-background split action of the form of 

S = — l -tr ! d 4 xd 4 9 (e- 9V V a e sV W 2 (e- 9V V a e 9V ) 

V J 

+ J d 4 xd 4 9 0e 9V > + j d A x [d 2 9 P(0) + h.c] , (B.2.28) 

where is background covariantly chiral superfield. 

After the gauge fixing procedure, we add to the action a gauge-fixing term (Sqf), 
Faddeev- Popov (Spp) and Nielsen-Kallosh ghosts terms (S^k)- All of them have the follow- 
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mg expressions 



Sgf 
Spp 

Snk 



-i t r f rfVfl(V 2 V)(vV) 
a J 

tr [ d A xd 4 9 



c'c - c'c + -(c' + c')[gV,c + c} + ... 



-tr I d A xd 4 6 bb 
a 



(B.2.29) 



The quantum gauge quadratic action can be written as [381 [6TJ (59J [60 



-tr I d 4 xd A 6V 



□ - iW a V n - iW V ( 



V 



□ 



-V aa V r 



We denote this kinetic operator as □ = □ — iW a V a — iW a ^ Also, for the background 
covariantly chiral superfields, we define operators n± by 



V 2 V 2 < 



n 



U-iW a V n - -(v a w r 



n-tW"V a --(V-W aJ 



2 

i 

2' 



(B.2.30) 



The covariant Feynman rules can be obtained from the partition function once we have 
introduced real and chiral sources J, j. This partition function can be written as 



Z = A + Aexp 



S, 



hit 



5 5 5 
6J : 8j' 8j 



exp 



(B.2.3i; 



where A, A + are one-loop contributions from real and chiral superfields (including ghosts). 
From this it is clear that covariant Feynman rules are similar to the usual ones except that 
for V-lines we have □ 1 as the propagator, and (ft fields in vertices are joined by — and 
usual D 2 factors at the vertices are replaced by V 2 factors. Making use of these covariant 
Feynman rules, in the background field approach we consider vacuum graphs with quantum 
vertices derived from S int , □ 1 and — Q+ 1 propagators and V 2 , V factors. The idea is to 
use the algebra of covariant derivatives and Bianchi identities to push, in each propagator, 
the covariant spinor derivatives to a given vertex. At this point, using the anticommutation 
relations, they can be integrated by parts or eliminated in favour of space-times derivatives. 
Finally, we have to apply the relation [59] 



5 12 V 2 V 2 5 12 = 5 



'12 



(B.2.32) 



in order to obtain free grassmanian 5-functions that allow us to evaluate 9 integrals. 
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Appendix C 



Gauge parameters and the RG 
equations 



In order to fix the gauge when quantizing the different theories that we have treated, a 
gauge fixing term that depends in an additional parameter a has been added to the action. 
Therefore, we have to take care of the running of this gauge parameter in the RG equations. 
In general, these equations applied to two-point functions are of the form 



where 7 is the anomalous dimension, (3 the beta function and 7 Q the function that takes care 
of the running of the gauge parameter. So, we need 7 a to verify the two-loop background 
RG equations. However, as in the models that we have considered the first dependence on 
a of the background two-point function can only arise at the one-loop level, we can obtain 
the relevant terms of the 7 Q (g) expansion by means of the evaluation of the one-loop RG 
equations for the quantum gauge field self-energies. 



M 



d 



+ /3^- + 7a^--2 7 r 2 = o, 



DM 



C.l Abelian examples 



C.l.l 



QED 



Figure C.l: One-loop QED diagram. Wavy lines correspond to gauge fields and solid lines 
to fermion fields. 
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The one-loop contribution to the quantum photon self-energy is shown in figure ICTTl This 
has an explicit expression of the form 



Y[ AA (! 1oo p) 



lie, 



(C.l.l) 



which is the same that we have found for the background gauge fields in (12.1.91) . Hence, with 
the same procedure as for the background fields, we find the following renormalized value 



127T 2 (47T 2 



lnx 2 M 2 e 2 
-□ = -5(x) 



x 



36tt 2 



(C.1.2) 



With this, the quantum gauge field two-point function expanded to one-loop order in a 
general gauge is 



5(x) 



9(4vr 2 



-8(x) 



e 2 lnx 2 M 2 
-□ 



3(4vr 



2\2 



1 



-d^dj(x) + 0(e 4 ) . 



As this function satisfies the usual RG equation 



-rAA _ q 



(C.1.3) 



(C.1.4) 



using the one-loop expression (1C.1.3|) we find the following one-loop values for 7^ and 7 a 

1 



lA{e) 
7« (e) 



3(4tt 2 ) 
2ft 



3(4tt 2 



e' + ... 



-e z + . . . . 



(C.1.5) 
(C.1.6) 



C.1.2 Super QED 




Figure C.2: One-loop Super QED diagram. Wavy lines correspond to gauge fields and solid 
lines to $ + or $_ propagators. 



For the supersymmetric extension of the previous model, the situation is very similar. We 
start by considering the kinetic part of the action for the quantum gauge field in a generic 
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gauge. This is of the form 



S, 



(2) 



i J d 4 xd 4 6 VD a D 2 D a V — i y d 4 xd 4 # (I)V)(5V) 

-- / d 4 xd 4 # VDIIiV - — I d A xd 4 9 vau v , 
2 ./ 2 2a ./ 



(C.1.7) 



where we have used the projection operators II i = —D a D 2 D a /n and Ho = (D 2 D 2 + 
D 2 D 2 )/D. As in QED, the one-loop renormalized contribution to the quantum gauge field 
self-energy is the same that we have evaluated for the background case (12.2.481) . So, the 
complete expansion to one-loop order of this self-energy is 



IX*) = ~nn ii( ,) - ijniww + 



s 2 Dn , n^fl + otf). (c.1.8) 



Thus, considering that this amplitude satisfies a RG equation of the form 

r = o, 



the values that we find for 7y and 7 Q are 

7v = 

7a = 



2(4tt 2 ) 



(4vr 5 



s + 



v + . . . 



(C.1.9) 



(C.1.10) 



C.2 Non-abelian examples 
C.2.1 Yang-Mills 




Figure C.3: One-loop Yang-Mills diagrams. Curvy lines correspond to gauge fields and 
dashed lines to ghosts. 

We begin writing the effective action as 

r = ~ J d 4 xd 4 y A^x)A b u (y)T^ ab (x - y) + 0(A 3 ) . (C.2.11) 

If we consider the part of the Yang-Mills lagrangian which depends only in the quantum 
fields A", in a generic gauge this is of the form 



\^ u F^ + ^l + 0(d,A;)(d u At). 



(C.2.12) 
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Notice that we have redefined the usual gauge parameter a as ^ = (1 + £). With this, the 
effective action can be written as 

r = I J d'xd'y A; [6 ab {-S^USix - y) - ^d v 8{x - y)) - U^ ab (x - y)] A a u {y) 

+0{A 3 ) . (C.2.13) 



At the one-loop level, as is shown in figure IC3|. we have contributions with gauge and 
ghost loops. We first obtain the fully expanded bare expressions (in Feynman gauge) and 
then renormalize them according to CDR rules. 



Gauge loop 



g2 jacd jbdc 



*xy 



<W^- %) + 8*M - D x p ) + 6 pa (d* -dl) 



x 



A 



9 2 C A 5 



ah 



[2(d„d v - 5 llu n)A 2 + lO^(A^A) - lOAd^A 



45, u d x (Ad x A)-25^A(nA)} . 



Ghost loop 



(C.2.14) 



-g 2 f adc f bcd A xy d* dlA xy = -g 2 C A 5 ab [d^Ad.A) - Ad^A] . (C.2.15) 
Adding the two previous results we find the total bare contribution to be 



tAA ahi 



9 2 C A 5 ab [d,d u A 2 - 5, V UA 2 + Ad,{Ad v A) - 25, u d\Ad x A) 

- AAd„d u A - VA(DA)] , (C.2.16) 



and with CDR identities it is straightforward to obtain the renormalized result as 

g 2 C A 5 ab 
g 2 C A 5 ab 



144vr 2 



- 8 pv u) 



15 \nx 2 M 2 

47T 2 X 2 



(C.2.17) 



So, Y ab y can be written as 



ri A n b (x) 



+ 



S^uSix) - ^dj(x) + 5 ab (d^ - 5 pu a) 
g 2 C A 



hg 2 C A _\nx 2 M 2 



-□- 



487T 2 (47T 2 ) X 2 



727T 2 (47T 2 



-S(x) 



+ 0(g 4 



(C.2.18) 
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Inserting this in the RG equation 



t^AA ab\ i 
1 fiv R k=0 — U 



(C.2.19) 



we can easily obtain the values for 7^ and 7^ as 



7e 

1A 



5C A 2 

-2^ 9 + 
5C A 2 

'48^ + 



(C.2.20) 



C.2.2 Super Yang-Mills 




Figure C.4: One-loop Super Yang-Mills diagrams. Wavy lines correspond to gauge fields 
and dashed lines represent ghosts. 



The kinetic term of the Super Yang-Mills action for quantum gauge fields in a generic 
gauge is 



S 



(2) 



l -tr J d A xd A yd A e V a (x, 6) [DU 1/2 + (1 + O^o] V a (y, 6)6 (x - y) , 



(C.2.21) 



where again the usual gauge parameter a is redefined as - = 1 + £. We have for these fields 
an effective action of the form of 



I> = J d A xd A yd A 9 V a (x, 6)Ty (2) (a; - y)V b (y, 6) + 



(C.2.22) 



The diagrams that contribute to the one-loop quantum gauge field self-energy are those of 
figure ICT41 As an example, we will detail the calculation of the ghost contribution (diagram 
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(a)). This is of the following form 



^l loop 

V(a) 



2 3 



d 8 Zl d 8 z 2 V a ( Zl )V a (z 2 ) {[D 2 2 D 2 2 P 12 ] [D 2 2 D 2 2 P 12 ] + [D 2 2 D 2 2 P 12 ] 

x [D 2 2 D 2 2 P 12 ] - 2 [D 2 2 D 2 2 P 12 ] [D 2 2 D 2 2 P 12 ] } 

9 2 C A 
2 3 



9 2 C A 



d 8 Zl d 8 z 2 V a ( Zl ) [D 2 D 2 V a (z 2 )} P 12 [D 2 2 D 2 2 P 12 ] 
d 8 z x d 8 z 2 V a ( Zl ) [D 2 D 2 V a (z 2 )] P l2 [D 2 2 D 2 2 P 12 ] 
d 8 z x d 8 z 2 V a ( Zl ) [D 2 D 2 V a (z 2 )] P l2 [D 2 2 D 2 2 P l2 ] 
d 8 z x d 8 z 2 V a { Zl )V a {z 2 )P 12 [UD 2 2 D 2 2 P 12 ] 

d 8 z x d 8 z 2 V a ( Zl ) [D«D a V a (z 2 )] P 12 [d 2 a& D 2 2 D 2 2 P 12 ] , (C.2.23) 



2 3 
9 2 C A 
2 2 

9 2 Ca 

2 2 
i9 2 C A 

2 2 

where in the second step we have used identity (I2.2.5ip . With this expression it is clear that 
we can apply the 5-function property (IA.5.541) that leaves a free grassmanian ^-function, 
allowing us to perform one of the 9 integrals. At this point, renormalizing and identifying 
x\ = x, x 2 = y we find for this contribution 



pl loop 
L V(a) 



9 2 C, 



d 4 xd 4 yd 4 9 V a (x,9) 



-□IIo - \au 1/2 



V a (y,9) A 2 xyR 



(C.2.24) 

Diagram (b) can be evaluated in a similar way [61]. The final result is found to be 



pl loop 
L V(b) 



9 2 C, 



J d 4 xd 4 yd 4 9 V a (x, 



--□n 1/2 + ^nn 



v a (y,9) A 2 xyR 



(C.2.25) 



The total contribution is then obtained as 



-\1 loop 
V 



pl loop i pl loop 
1 V (a) + 1 V (6) 



= "lS / <ftrfVtf ^ "^'^ 0)d H 1 -yf - (C - 2 - 26) 

From the effective action, defining Ty^ 2 \x) = 8 ab Ty\x), we have a RG equation for 
quantum gauge fields of the form of 



d d d 

M- h 0(g)— + 7 f (o)— - 27y 



o 



(C.2.27) 



and we find Fy (x) to be evaluated as 



rf (*) = -^(x)nn 1/2 - h(x){i + e)nn - ^^n^£nn lft + o( g 4 ) . 



16(4tt 



x/- 



(C.2.28) 
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Straightforward operations lead us to obtain the following values for 7y and 75 
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Appendix D 
Explicit calculations 



D.l Integrals with overlapping divergences 

Here we will show how we can obtain the different expressions for the integrals with over- 
lapping divergences presented in section I1.3.2L 



D.1.1 Notation 

Be begin by discussing some notation that we use in these calculations. As we did when we 
listed integrals of section [i.3.2l we will write the final results in terms of a variable z = x — y. 
Some intermediate local results will be found to be multiplied by a constant termed a which 
value is a = 67r 4 £(3)/(4vr 2 ) 4 - 



Integral relations 

To obtain some of the results, we use the following exact integral relations 




+ J d\d A v A xu (d^A xv )A yu (d y v dlA yv )A. 
+d x u [ d A ud A v A xu {d;A xv )A yu {d y x A yv )A 



■uv 




(D.1.1) 
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d y x J d A ud A v A xn A xv (d y x d y A yu )A yv A uv 

\dl J d A ud A v A xu A xv (aA yu )A yv A l 



-J d A ud A v A xu A xv (aA yu )A yu A uv 

+\dld\ j d\d 4 v A xu A xv (d y x A yu )A yv A uv (D.1.2) 



d\ I d A ud A v A xu A xv {dldlA uy ){dlA yv )A uv 

= \dldld\ J d A ud A v A^A^A^A^A* 



+- A dlu j d A ud A v A xu A xv {dlA uy )A yv Auv 



~U j d A ud A v A xu A xv (dldlA yu )A yv Auv (D.1.3) 

To prove these relations, one has only to perform the derivatives and expand the different 
terms. 



D.1.2 Calculations 

Evaluation of H[l, 1 ; 1,1], H[d^, 1 ; 1,1] and H[l, d x ; 1, d x ] 

These integrals are obtained by means of Gegenbauer Polynomials [7J [37] . 

Evaluation of d x H[l, ; l,d\] 
Contracting (ID.l.lj) with 5 vX we obtain 



dtH[l,d,;l,d x ] = ~cP#[l,l; □ , 1] - if [1, ^ ; 1, □] - ~n#[l, 0„ ; 1,1]. 

(D.1.4) 

In order to renormalize, we have only to write the previous expressions in terms of the integral 
form I l (x) (remember DA = —5), and use the results found in section [1.3.11 Therefore we 
find 

dZH[l, d, ■ 1, d x ] = \d,{AI l ) - {Ad,I l ) + - A {d,8) 
R 1 n \\xvz 2 M 2 

- + (D.1.5) 



32(47r 2 ) 3 " z 2 
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Evaluation of H\d^d x , d x ; 1,1] 

In this case, no integral relation is used to write H[dfj,d x , d x ; 1, 1] in terms of I 1 



H[d^d x ,d x ; 1,1] 



~d?H[d x ,d x ; 1,1] 

±d£d* x H[l,d x] 1,1]-~«9^[1,D; 1,1] 



32(4vr 



±ln 2 z 2 M 2 -lnz 2 M 2 



2^3 A* 



+ . . . 



(D.1.6) 



Evaluation of <^#[1,<9 M ; ^<9 A ,1], df#[l,d A ; 8^,1] and #[1, <9 A 5 9 A 9^, 1] 
First of all, the third integral (ll.3.54p will be evaluated with relation (ID. 1.11) 



H[l,d x ; d x d v ,l] - 

Using the previous results 
H R [l,d x ; d x d u ,l] 



±d$H[l,d x ; l,d v ] + ±0£ff[l,l; l,d x d v ] 
+^H[1, d x ; 1, d x ] + \c%d* x H[l, d x ; 1, 1] . 



±ln 2 z 2 M 2 -l\nz 2 M 2 
d n V* -i + . . . 



32(4vr 



2^3 A* 



(D.1.7) 



(D.1.8) 



However, this integral along with the other two, can be obtained with other method. 
The idea is to apply the CDR decomposition fll.2.38j) into trace part, traceless part and 
additional local terms to the divergent subdiagram (d^d^A yu )A yv A uv . Considering this in 
the general integral 



d*ud*v A xu (d x p A xv )(dyayA yu )A yv A uv 
4 -h^Ad^u-J^d^ 



16 



256tt 



(4tt 



2^5 



'peer _R 



(D.1.9) 
(D.l. 10) 



where I pea stands for the traceless part. The one-loop ambiguity fixed by CDR is reflected 
in the second term of (ID. 1.91) that at two loops has become a logarithm of the scale. In the 
renormalization of the traceless part normal differential renormalization will be used, leaving 
ambiguities (local terms) not fixed. 
The expression for I pea is 



1 pea 



pea R 



j) e-^a-^p 1 a^P - • 1 



X° 



X £ r _|_ X v r 
6 °P a + x 6 °P £ 



X 



\R 



(D.l.ll) 
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or in the form of the integrals being discussed 



I\\fi R 
dxIftXu R 

dxlx^u R 



--(An 2 ) 2 (3A-B)d,A 2 R 

1 „ „ . , 2 



-(4tt 2 ) 2 (3A- B) 
(An 2 ) 2 (3A-B) 



-d,d u A 2 R + -(47r 2 )^A| 



-^^A| + i(47r 2 )^A 



(D.1.12) 
(D.1.13) 
(D.1.14) 



The value of (3 A — B) is easily obtained using (1D.1.12|) . because this corresponds to 
integral fll.3.54j> that was obtained previously. I.e. 



d 4 ud 4 v A xu (d x x A xv )(d y x d y u A yu )A yv A uv \ R 
1 



32(4vr 



2\3 V 



Un 2 z 2 M 2 -Unz 2 M 2 
d„ □ ^ -5 



+ . . . 



2 Z 



32(4vr 



2^3 



±ln 2 z 2 M 2 + ±lnz 2 M 2 
-* 

z l 



16 



(4tt 2 ) 5 



I p£ffii , (D.1.15) 



which implies that 



(3A - B) 



3n 4 



(D.1.16) 



With this result the evaluation of ( 11. 3. 561) and ( 11.3.581) are straightforward 



^[1,^; d v d x ,l] 



256vr 2 



16 

(47T 2 ) 5 ^ ^ 



32(4vr 



2\3 



±ln 2 z 2 M 2 + ±ln;z 2 M 2 -±lnz 2 M 2 
§ -§ + 5 uv □ □ — — 



2 Z 



(D.1.17) 



d x x H R [l,d x ; 8^,1} 



-^(A^-^DAI- 



16 

(47T 2 ) 5 X X ^ uR 



1 



32(4tt 2 ) 3 
+ ... 



-±ln,2 2 M 2 |lnVM 2 + f lnz 2 M 2 



z^ 



2" 



(D.1.12 
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Evaluation of dfH[l, d\ ; 1, d^d u ] 

Using integral relation ID. 1.11 we can write this contribution in terms of others previously 
obtained. Explicitly, we find 

d x x H[l, dx ; 1, d^} = d x x H[l, d x ; d^, 1] - 8%d£H[l, d x ; 1, d v ] 

-dZd£H[l,d x] 1,1] 
R 1 [ Unz 2 M 2 Un 2 z 2 M 2 + Unz 2 M 2 

32(4tT 2 ) 3 Z = + V UU ^ 

+ ... 

(D.l. 19) 

Evaluation of H[l, ; 1, <9j,] 

Considering (I1.3.56P and applying (ID. 1.1ft we can put this as 

dtH^d,- d u d x ,l] = ±UH[l,dp; l,d v ) + ±d£d%H[l,l; 1,3„0 A ] 

+ifl?5Jjy[l,0„; l,d A ] + ^OFf[l,^; 1,1]. (D.1.20) 

Remembering previous results 

n Rr o o n i 1 ~~ ±lnVM 2 + ±lnz 2 M 2 

^[i,a M ;^ Aj i] = ^-^a^nJ 

+ ^□^[1, <9 M ; 1, 0„] + (local terms) (D.l. 21) 



So that 



i — -\nz 2 M 2 
DH^d,;!^] = ^y^ u nn^- 2 + ... (D.1.22) 



Evaluation of H[l, 1 ; d^d^, 1] 

Using fll.3.58p . fll.3.59|> and the identity 

DH[1, 1 ; d„d v , 1] = d x x H[l, d x ; d^d v , 1] + d£H[d x , 1 ; d^, 1] , (D.1.23) 

we can easily arrive to 

1 ±lnVM 2 + f ln^ 2 M 2 

= + ••• (D-1.24) 
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Evaluation of d%H [1, 1 ; d x d v , d^) 

Using (ID. 1.31) and (11 .3.60[) we can write this as 

%H[1, 1 ; d x d u , d,} = ±DH[1, 1 ; d,d u , 1] - ^c^#[l, 1 ; d x , 1] 

-^□#[1,1; ^,1]. 



(D.1.25) 



Hence, we have only to use previous results to find 

i r §: 

32(4tt 2 ) 3 M " 



o Rr o o o i 1 iln 2 z 2 M 2 + |ln^ 2 M 2 

a A a„,9 M ] = ———^□□8 ^ + ... (D.1.26) 



Evaluation of dfH[l, 1 ; d^d u ,d x 



With ( ID. 1.31) and ( 11. 3. 531) .this contribution can be evaluated with the same procedure of the 



previous one. So, we have 



1„„ „, „ 1 



<^tf[l, 1 ; 9^, 9 A ] = -d° v UH[l, 1 ; <9 M , 1] + -d£d£%H[l, 1 ; 9 A , 1] 

+^^#[1,1; d x d v ,l] 

1 ±ln 2 z 2 M 2 + fln^M 2 

d^U* -5 + . . . (D.1.27) 



32(4vr 2 ) 3 " " z 



Evaluation of H[l, d^d), ; d v d x , 1] 



In this case the CDR decomposition into trace+traceless+local terms ( 11. 2. 38ft will be used 
again, as in ( II. 3. 56ft and ( jl .3.58|) 

H R [l,d,d x ; d u d x ,l] = ~{^dJ l ) R -- A [^d v - l -8^U)I 1 ] R 

' ; (A^A),--^[A(^-k,D)A] fi 



128tt 2V m " ,n 128tt 2L_v ~^ 4' 
64 

+ J^y I » XvXR ' (D.1.28) 

where I^a^a stands for the integral with the traceless parts. This was calculated in [TJ] , and 
the result was found to be 

64 5 13 

t — ^tIuXvxr = — ^d u d v A% + — 5 UU A 3 R . (D.1.29) 
(47r 2 ) 5 M 96(47r 2 ) 48 



D.2. UV AND IR DIVERGENT INTEGRALS 
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Adding up all the terms, it is easy to arrive to 

1 



H K [l,d^ ; d v d x ,l] 



32(4tt 



2\3 



±ln 2 z 2 M 2 -^\nz 2 M 2 
d^n^ 



zs 



-±ln 2 z 2 M 2 - % In z 2 M 2 

+ V DD ^ 5"^ 



Z^ 



+ ... 



Evaluation of H[l, <9 M <9 A ; 1, d\ d„] 

In this case, applying flDTTT) . fll. 533) , fll.3.54]> . fll.3.56h and fll.3.63|) we get 



if[l,<9 M <9 A ; 



[1, d^d\ ; 0„0 A , 1] - 8£ff[l, <9 M <9 A ; 1, <9„] 
-%H[1, d„d x ; 1, 9 A ] - dZd%H[l, 8,8 X ; 1, 1] 
Mn 2 z 2 M 2 + ||ln^ 2 M 2 



32(4tt 



2^3 



-±\n 2 ^ 2 M 2 -iilnz 2 M 2 

+ <V DD — r^ 2 



+ 



(D.1.30) 



(D.1.31; 



D.2 UV and IR divergent integrals 

In this section we will detail the renormalization of three relevant integral expressions that 
appear when considering two-loop diagrams made up by the insertion of a one-loop propaga- 
tor. Examples of this are diagram (a) of QED or diagram (b) of the Yang-Mills case. These 
expressions are 



I\x-y) 



D.2.1 Renormalization of 7 C 



d 4 ud 4 v A xu A yv A 2 uv 



J d 4 ud 4 v A xu A yv (A uv c%A uv ) 
d\d 4 v A xu A yv (A uv d^d v u A uv ) 



(D.2.32) 



The renormalization of 1° is detailed in section 11.1.31 and here we only recall the final 
renormalized result found there. 



i°r(x - y) 



32(4tt 



2\2 



[In 2 x 2 M 2 R + 2 lnx 2 M 7 y 1 - lnx 2 M 2 )] + . . . 



(D.2.33) 
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APPENDIX D. EXPLICIT CALCULATIONS 



As in diagram (b) of the Yang-Mills theory we have contributions of the form An/ and 
Adf+dvl , we have to evaluate them. For the first one, it is clear that 



2 

av 



d\ A vy A 2 xv 



-/V-y) 

1 \n{x - y) 2 M 2 
~4(4tt 2 ) 2 (x-y) 2 



(D.2.34) 



where we have used the renormalized value found for I . Finally, in order to obtain Ad^d u I° 
we have to consider (ID. 2. 331) and apply usual DiffR. With this, we find 



[Ac^cUVx) 



1 



32(4tt 2 ) 3 
+ ... 



„ „ lnx 2 M 2 p 1 In 2 x 2 M 2 + \ In x 2 M 2 

2 + " 



x- 



(D.2.35) 







D.2.2 Renormalization of I\ 

The renormalization of J° is straightforward, once we recall that CDR imposes 7° R — \&^I R . 



D.2.3 Renormalization of I 



o 

[IV 



Applying CDR to the subdivergence we find 

1% = ~Jd 4 ud*vA xu A yv (d„d v -~5 fll/ n)(Al v ) R + 

J d 4 ud 4 v A xu A yv (d ^ - S^ lu a)S(u - v) 



288v 



12 



72(47r 



72(4vr 2 



With this we can evaluate the expression that appears in diagram (b) of the Yang-Mills case 
(AI® U ). We find the following result 



32(4?r 



2\3 



Mnx 2 M 2 



-±lnx 2 M 2 



dpd^ 5 + 6 



+ ... (D.2.36) 
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Renormalizacion diferencial 

Renormalizacion diferencial [7] es un metodo de renormalizacion en el espacio de posiciones 
que consiste en sustituir expresiones que son demasiado divergentes para tener una trans- 
formada de Fourier bien defmida, por derivadas de otras expresiones menos singulares. Asi, 
por ejemplo 1/x 4 no tiene una transformada de Fourier bien defmida, y renormalizacion 
diferencial proponer reemplazarla por la solution de la ecuacion diferencial 

— = uG{x 2 ) iVO, (37) 



r 1 



que es 



x 



X 



1 In x 2 M 2 



Notar que se ha introducido una constante con dimensiones de masa M que parametriza la 
ambiguedad local. Debido a que un cambio en M puede ser reabsorbido en un reescalamiento 
de la constante de acoplamiento, esto sugiere que las amplitudes renormalizadas satisfacen 
ecuaciones del grupo de renormalizacion, con M jugando el papel de escala del grupo de 
renormalizacion. 

Aunque en este trabajo se traten solo teorias sin masa, renormalizacion diferencial puede 
ser aplicada sin ningiin problema a teorias masivas, ya que las masas solo alteran el compor- 
tamiento a larga distancia de los correladores [7] [TT]. 

Renormalizacion diferencial se puede aplicar para renormalizar diagramas de orden ar- 
bitrario en teoria de perturbaciones. En concreto, en [13] se expone una implementation 
sistematica de renormalizacion diferencial a cualquier orden en teoria de perturbaciones. En 
general, cuando se aplica renormalizacion diferencial a un calculo a orden superior, apare- 
cen nuevas escalas correspondientes a la renormalizacion de los distintos subdiagramas que 
forman el diagrama completo. 

Tambien es relevante senalar que, aplicando renormalizacion diferencial en espacio de 
momentos, se pueden renormalizar divergencias IR. Asi, por ejemplo 



1 1 1 lnp 2 /M 2 R 

J R 4 P 



P 

Sin embargo, a la hora de renormalizar una teoria que tenga divergencias IR y UV, se ha 
de tener en cuenta que ambas renormalizaciones han de estar desacopladas, implicando por 
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ello que ambas escalas (IR y UV) han de ser independientes. En concreto, en este trabajo 
se discute una expresion divergente IR de la forma In p 2 j M 2 / p 4 , donde M es una escala 
UV producto de una renormalizacion previa en espacio de posiciones. En este caso, la 
independencia de las escalas se consigue en el momento que imponemos la relacion 



M- 



6_ 

6M 



\np 2 /M 2 
p 4 



M 



5 \np 2 /M 2 
6M p 1 



(40) 



R 



Esto se satisface ajustando los terminos locales en los que aparecen ambas escalas, obte- 
niendo entonces la siguiente forma para la expresion renormalizada mas general [T5] 



In p 2 /M' 
p 4 



21 



1 



R 



\n z p 2 /Mf R + 2 In p 2 /Mf R (1 + \np 2 /M 2 
p 2 



+ {a IR lnj^ + b IR )6(j)) . 

(41) 

Una de las caracteristicas mas importantes de renormalizacion diferencial es que la in- 
variancia gauge se preserva. Sin embargo, debido a las ambiguedades que se generan en el 
metodo de renormalizacion, se han de imponer en cada calculo (con una teoria gauge) las 
identidades de Ward de forma explfcita, de tal manera que se fije el esquema de renormal- 
izacion. El hecho de que se preserve invariancia gauge se refleja en que siempre es posible 
satisfacer estas identidades con las expresiones renormalizadas (excepto por supuesto las 
anomalias) . 



Renormalizacion diferencial restringida 

Para evitar la necesidad de imponer las identidades de Ward explfcitamente en cada calculo, 
se desarrollo Renormalizacion Diferencial Restringida (RDR) [21]. Este metodo consiste en 
proporcionar una serie de reglas que a priori fijan toda la ambiguedad inherente al proceso, 
de tal modo que las expresiones renormalizadas sean directamente invariantes gauge (no es 
necesario imponer las identidades de Ward). Las reglas que impone RDR son: 

1. Reduction diferencial 

• Funciones con comportamiento peor que el logaritmico se reducen a derivadas de 
(como mucho) funciones logaritmicamente divergentes sin introducir constantes 
dimensionales extra. 

• Expresiones logaritmicamente divergentes se escriben como derivadas de funciones 
regulares, introduciendo una unica constante M, que tiene dimensiones de masa 
y juega el papel de escala del grupo de renormalizacion. 

2. Integration por partes formal. No se tienen en cuenta los terminos de contorno diver- 
gentes que aparecen cuando integramos por partes. En relacion a esto, la renormal- 
izacion y la diferenciacion deben ser dos operaciones conmutativas: si F es una funcion 
arbitraria, entonces [dF] R = d[F]R. 

3. Regla de renormalizacion de la funcion delta 



[F(x, xi, . . . , x n )5(x - y)) R = [F(x, x t , . . . , x n )} R 5(x - y) (42) 
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4. Validez de la ecuacion del propagador 

[F(x,xi, . . . ,x n )(D - m 2 )A rn (x)] R = -[F(x,xi, . . . , x n )S(x)] R (43) 
donde A m es el propagador de una particula de masa m y F una funcion arbitraria. 

Aplicando estas reglas, obtenemos un conjunto basico de expresiones renormalizadas. 
Por lo tanto, el proceso de renormalizacion consta de dos partes: en un primer momento, 
se realizan todas las contracciones de indices (RDR no conmuta con esta operation) y se 
escribe la expresion desnuda en terminos de estas funciones basicas. En un segundo paso, se 
sustituyen estas funciones por sus valores renormalizados. 



Aplicacion de RDR a calculos a dos bucles 

Aunque RDR se ha desarrollado solo para calculos a un bucle, proporciona information util 
cuando tratamos calculos a dos bucles. Veremos que aplicar RDR fija umvocamente los 
coeficientes de todos los logaritmos de las escalas en la expresion a dos bucles renormalizada, 
que son los terminos que necesitamos para evaluar la ecuacion del grupo de renormalizacion. 
Es por ello que, al obtener las expresiones renormalizadas a dos bucles, no tendremos en 
cuenta los posibles terminos locales que se generen. Distinguiremos dos situaciones diferentes: 
diagramas con divergencias anidadas y diagramas con solapamiento. 

Divergencias anidadas 

En este caso, empezamos imponiendo RDR a la subdivergencia. Al hacer esto, fijamos los 
terminos locales a un bucle que tenemos en el diagrama, junto con los logaritmos de las 
escalas a un bucle (lnx 2 M 2 ). Entonces, al considerar la expresion completa del diagrama 
y aplicar renormalizacion diferencial normal, nos encontramos que todos los coeficientes 
que corresponden a logaritmos de las escalas estan umvocamente determinados, ya que los 
terminos locales a un bucle (que se promocionan a logaritmos) han sido fijados por RDR. Con 
este procedimiento, realizamos la renormalizacion de diferentes expresiones que contienen la 
siguiente integral que se utiliza a lo largo del trabajo: F = f d 4 uA xu A 



2 

yu 



[Al \ R (x) = - ——— a + (termm. locales) 



1 ln 2 x 2 M 2 + 21nx 2 M 2 

32(4^° x~ 2 

fAO rl , . , 1 ln 2 x 2 M 2 + lnx 2 M 2 . . , , . 
[Ad^I \ R (x) = - 64(47r2)3 ^ D ~2 + (termm. locales) 



[Ad,d u I l ] R (x) 



1 



96(4tt 2 ) 3 



In 2 x 2 M 2 + l\nx 2 M 2 



d,d v u 



1 ln 2 x 2 M 2 + ^lna; 2 M 2 
" I^ DD ^ 



+ (termin. locales) 



1 In x 2 M 2 

[Aai 1 ] R ( x ) = 32 ^ 47r 2)2 DD — ~2 — + (termin. locales) . (44) 
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Divergencias con solapamiento 

En el caso de divergencias con solapamiento, la situation es mas complicada, ya que muchas 
veces es dificil reconocer las expresiones a un bucle a las que hay que empezar aplicado 
RDR. Por lo tanto, lo que hemos hecho es obtener una conjunto completo de integrates 
renormalizadas con solapamiento, con como mucho cuatro derivadas actuando sobre los 
propagadores y dos indices libres. Con esta lista podemos obtener la expresion a dos puntos 
renormalizada de cualquier teoria con acoplos derivativos; lo que implica que, aplicando el 
metodo del campo de background, nos permite obtener la funcion beta. Para evaluar las 
integrates hemos empleado basicamente dos metodos: 

• Mediante igualdades integrates rescribimos las integrates en terminos de otras que ten- 
gan un d'almbertiano actuando sobre uno de los propagadores. Esto permite obtener 
la integral como suma de contribuciones de integrates con divergencias anidadas, en 
las cuales se aplica lo que hemos discutido anteriormente. 

• Utilizamos la descomposicion en parte con traza y sin traza que impone RDR (en la 
que se ahade un termino local fijo). 

Esta lista esta escrita en terminos de una expresion H que hemos definido como 

H[O u 2 ; 3 , 4 ] = J d 4 ud 4 v (0*A XV ) (0%A XV ) (O^A yu ) (0\A yv )A uv , (45) 
siendo Oi un operador diferencial. 

ifS[MiMl = W A " A (46) 

H " ld - = T^By {d - A) 55 5 a " A (47) 

1 In 7 2 M 2 

H«[l,d x ;l,d x ] = + (48) 



1 \\nz 2 M 2 

d'tH R [l,d,;l,d x ] = + (49) 

1 ±ln 2 z 2 M 2 + ^\nz 2 M 2 



^[M; dM = ^y 3 d^ - 2 4 +... (so) 

Rr o o o i 1 o l M 2 z 2 M 2 -l\nz 2 M 2 , x 

H R [l,d x ;d x d„l] = ^yA n " ^ + ••• ( 51 ) 

Rr 1 -iln 2 z 2 M 2 -In z 2 M 2 , x 

H R [d,d x ,d x ;l,l] = ^-^dp-Z + ... (52) 



d x x H R [l,d,; d u d x ,l] 



1 



2\3 



32(4vr 



iln 2 z 2 M 2 + \\nz 2 M 2 



Z 2 



1 In z 2 M 2 ^ 



4 



+ ... (53) 
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H R [l,dpi l,d v 



d%H R [l,d x ; l,d^d u 



H R [1,1- 0^,1} 



d x x H R [l,l; d x d v ,d„ 



d x x H R [l,l; d^dx 



H^dfix ; d u d x ,l] 



H R [l,d^d x ; l,d v d x 



1 c -Unz 2 M 2 



32(4tt 2 V> ;:r 
1 



+ . . . 



z c 



32(4tt 2 ) ; 



n n -±lnz 2 M 2 



±ln 2 z 2 M 2 + |lnz 2 M 2 
+<V □ □ s- 5 



+ 



32(4tt 



2^3 



klnz 2 M 2 



iln 2 z 2 M 2 + |ln^ 2 M 2 



+ . . . 



iln 2 z 2 M 2 + f lii2 2 M 2 



32(4vr 2 ) 3 



+ . 



32(4vr 2 ) 3 ^ 



Mn 2 z 2 M 2 + |lnz 2 M 2 



+ 



32(4tt 2 ) 3 " 
1 



dud„n. 



1 ln 2 2 2 M 2 + f ln2 2 M 2 



+ 



32(4vr 



2^3 



±ln 2 ^ 2 M 2 -|,lnz 2 M 2 

6 do 



2^ 



-^lnVM 2 - ffWM 2 
+<V □ □ — s- 22 



2 Z 



+ • 



32(4vr 



2^3 



iln 2 ^ 2 M 2 + §lnz 2 M 2 

C^D- 6 -36 

Z z 



-Mn 2 z 2 M 2 -%\nz 2 M 2 

+<W □ □ 5-^ 



(54) 



(55) 



(56) 

(57) 
(58) 



(59) 



(60) 



Ejemplos abelianos 



QED 

La renormalizacion diferencial de QED a dos bucles fue realizada en [13] empleando identi- 
dades de Ward para relacionar las escalas. En este trabajo reharemos este calculo empleando 
RDR a un bucle, que nos permitira obtener la funcion beta de QED sin necesidad de emplear 
dichas identidades. Utilizaremos para los calculos los convenios de [Hj, por lo que la accion 
que consideramos es 



(61) 



donde ip es el campo fermionico y se expresa en terminos del campo gauge como 
F^(x) = d^A v (x) - d v A^(x). 
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A diferencia que en |T3], realizaremos los calculos con el metodo del campo de background. 
Con este metodo, se divide el campo gauge en dos contribuciones — ► A^+B^: una cuantica 
(An), que es la variable de integration en el la funcion de partition y por lo tanto sobre la que 
se fija el gauge, y otra background (B^), en la que se mantiene invariancia gauge explfcita. 
Esto tiene multiples consecuencias, entre las que destacamos el hecho de poder obtener la 
funcion beta a partir solo de la funcion a dos puntos. 

Un bucle 

La autoenergi'a del foton a un bucle tiene una expresion desnuda de la forma 

-{lefTr [ 7m 7XA7,7 ct ^A] 
= -e 2 Tr [ 7/ ,7 A 7,7l (d x (Ad £ A) - Ad x d a A) , (62) 

a partir de la que, de acuerdo a las reglas de RDR, obtenemos el siguiente valor renormalizado 



tt(1 bucle) 

fjLV 



127T 2 (47T 2 ) 



lnx 2 M 2 



36tt 2 



S(x) 



(63) 



Obtenemos ahora la autoenergi'a del electron, ya que es una insertion dentro de uno de 
los diagramas a dos bucles. La expresion para esta contribution es 



e 2 ~f^A fMU (x)j x d x A(x)~f u 



(64) 



que utilizando el propagador del foton en un gauge general y las funciones basicas de RDR 
se renormaliza como 



e 2 7 A 



1 



4(4tt 



2^2 



_ lnx 2 M 2 . 
d x n—p— + (a-l) 



1 



In x 2 M 2 



4(4tt 



2^2 



X/ 



Autoenergia del foton background a dos bucles 



+ 



1 



167T 2 



d x S(x] 



(65) 





Figure D.l: Diagramas a dos bucles de QED. 



En primer lugar, senalar que los calculos a dos bucles los realizamos en el gauge de 
Feynman. Esto es debido a que la selection de este gauge en concreto no afecta a la ver- 
ification de las ecuaciones del grupo de renormalizacion a dos bucles, como se vera al dis- 
cutir dichas ecuaciones. Por lo tanto, el propagador desnudo del fermion en dicho gauge es 
= -2e 2 j x Ad x A(x). 
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Comenzando por el diagrama (a), este tiene la siguiente expresion desnuda 

K% a \x - y) = -(^e) 2 J d A ud% Tr [ ltll \-d^A xu )^\u - ^) 7 £ (-^A„ J/ ) 7 ^(-^A J/x )] . 

(66) 

Para simplificar la notation, definimos 7° como I® — J d A ud A v A xu A yv (A uv d ll A uv ) , que nos 
permite escribir esta contribucion como 

K% a \x) = e 4 [-32(^A)9 A ^/ A Vl6^(^A)9 A 9 CT / A Vl6(^A)n/0-8V(9 p A)n/a . 

(67) 

Por lo tanto, la renormalizacion de esta expresion pasa por estudiar la renormalizacion de 7°. 
Es facil demostrar que esto se puede escribir en terminos de la renormalizacion de la integral 
7 1 defmida previamente. Asi, d lx I^ R = — 17^ y 01® R = — \d^I R . Entonces, la contribucion 
renormalizada de este diagrama es 



r(2«)/ > e 



\n 2 x 2 M 2 - flnx 2 A7 2 ._ In 2 x 2 M 2 + f lna; 2 M 2 ' 



^ ^ + ^ 



1 V fiv / 24(4tt 2 ) 3 

+ (termin. locales) . 

(68) 

Pasamos ahora a renormalizar el diagrama (b). En este caso, lo que tenemos es un 
diagrama con divergencias de solapamiento. La expresion basica de este diagrama es 

rt^Kx-y) = -{ief J d^vTr^^d^^Yi^dlKyhv 
x (7 A ^A J/ ,) 7p (7^ CT A ra )A ra ] , 

(69) 

Mediante identidades de las matrices 7, e integrando por partes las derivadas que actuan 
sobre A xu y A yv , podemos rescribir esto en terminos de las expresiones 77 como 

n£ 6) = e 4 [ -8^D77[1, d x ■ d x , 1] + 16cF77[l, d u ; □, 1] - 8^^77[1, d x ; □ , 1] 

-16<9*77[1, d x ; d x d v , 1] + 16^77[1, d x ; d»d v , 1] - 16«^77[1, 3„ ; d x d v , 1] 
-16^77[1, □ ; d v , 1] + 8 <V^77[1, ° 5 -9a, 1] + 16^77[1, d x d„ ; d v , 1] 
-16^77[1, d^ v ; d x , 1] + 16cP77[l, d x d^ ; d x , 1] - 1677[1, □ ; d,d u , 1] 
+85^77(1, □;□,!] + 3277[1, d^d x ; d v d x , 1] - 1677[1, dfi v ; □, 1] ] . (70) 

Aquellas contribuciones que tienen un d'alamebertiano se pueden rescribir en terminos de 
la integral 7 1 , por lo que su renormalizacion es inmediata. El resto se encuentran dentro de 
la lista de integrates con divergencias de solapamiento (o pueden ser facilmente expresadas 
en terminos de esas integrates), con lo que solo hay que sustituir el valor renormalizado 
correspondiente. Por lo tanto, la contribucion renormalizada del diagrama (6) es 



4 



\n 2 x 2 M 2 + ^\nx 2 M 2 In 2 x 2 M 2 + % \nx 2 M 2 

dfivU -2 5^00 



X 2 X 2 



>'"" K ' 12(4tt 2 ) 3 

+ (termin. locales) . (71) 
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Entonces, la contribution total renormalizada a dos bucles de la autoenergfa del foton es 

n<2*(a0 = 2ng,-)(x) + ng, 6 )(x) 

p 4 1 r 2 M 2 

(d,d u -5, u U)n^f- + ... (72) 



4(4vr 2 ) 3V " " y x 2 

donde . . . corresponde a los terminos locales que no estamos teniendo en cuenta. 

Ecuacion del grupo de renormalizacion 

Estudiando la ecuacion del grupo de renormalizacion a un bucle para la funcion a dos pun- 
tos de los campos cuanticos, obtenemos que la funcion que corresponde a la variation del 
parametro gauge en dicha ecuacion , ~y a (e)d/da. Esta funcion tiene una expansion de la forma 
Ta( e ) = ~ 3(4^2) e 2 + 0(e 3 ). Esto justifica el poder realizar el calculo a dos bucles en el gauge 
de Feynman, ya que 7 a (e)<9/<9a actuando sobre cualquier diagrama a dos bucles no afecta 
a la verification de las ecuaciones del grupo de renormalizacion (es de orden superior en e, 
ya que la primera dependencia en a de la autoenergfa del foton background aparece a dos 
bucles). 

En cuanto a los campos background, si defmimos = -B' tenemos que la dimension 
anomala de este nuevo campo es nula, ya que la renormalizacion de B^ y de la constante 
de acoplamiento verifican la relation Z £ \/~Zb = 1. Por lo tanto, la ecuacion del grupo de 
renormalizacion que verifican estos campos es 



M — + 8(e)— \ T BB W = 



(73) 



con 

rjffa-y) = (d tl d v -5^n)SW(x-y)-Il^(x-y). (74) 

Con las contribuciones renormalizadas que tenemos para IT^, obtenemos que la funcion 
beta a dos bucles de QED es 

m - ^ + ^* + OW. (75) 

SuperQED 

Pasamos ahora a obtener la renormalizacion de la extension supersimetrica del caso anterior, 
SuperQED. Aquf, aplicamos los convenios del superespacio definidos en [38], con los que la 
action de SuperQED es 

S = J d 4 xd 2 9 W 2 + J d 4 xd 4 9 $ + e 9y $ + + J d 4 xd 4 9 $_e^ y <l>_ , (76) 

donde W a es un supercampo quiral, que se expresa en terminos del supercampo real gauge 
V y de superderivadas covariales D a como W a = iD 2 D a V. Q± son supercampos quirales 
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de materia. Es importante senalar tambien que se puede definir teoria de perturbaciones 
en el superespacio: en este caso tenemos diagramas defmidos en terminos de superpropa- 
gadores Pij = A^j, donde A^ es el propagador usual y es la funcion delta de variables 
grassmanianas. 

Al igual que en el caso de QED, realizamos los calculos en el gauge de Feynman (se 
justificara posteriormente su uso) y con el metodo de campo de background. El supercampo 
gauge V se divide por lo tanto en dos contribuciones V — > V + B: V es el supercampo gauge 
cuantico y B el background. 



Funcion a dos puntos del campo B a un bucle 



Figure D.2: Diagrama a un bucle en SuperQED. 



En este caso, la expresion desnuda es 



r « = »- I d 8 Zl d»z 2 B( Zl )B(z 2 ) 



DfP 12 D 2 2 



D 2 2 P 12 D\ 



que, con el algebra de derivadas covariantes, puede ser rescrita como 
T? } = 9 1 Jd i xd 4 yd 4 9B(x,9)[D 2 D 2 B(y,9)]A 2 xy 
+ ?1 J d 4 xd 4 yd 4 9 B(x,9)B(y,9)A xy (nA xy ) 
_!|! J d 4 xd 4 yd 4 9 B(x,9) [D«D a B(y,9)} A xy d y a .A xy . 

Aplicando RDR, obtenemos el siguiente valor renormalizado 

,2 



9 



R 



16(4tt 



2^2 



ln(x - y) 2 M 2 
(x - yf 



d 4 xd 4 yd 4 9 B(x,9) [D a D 2 D a B(y, 9)] □ 



(77) 



(78) 



(79) 



Funcion a dos puntos del campo B a dos bucles 

Omitiendo la expresion desnuda de cada diagrama en terminos del superpropagador P^, 
tras aplicar el algebra de superderivadas covariantes tenemos las siguientes contribuciones 
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Figure D.3: Diagramas de SuperQED a dos bucles. 



no renormalizadas 



i(2o) 



r 



(26) 



,(2c) 



^2J 



4 r 

9 — I d 4 xd 4 yd 4 9 B(x, 0)B(y, 9) [□(A/ 1 ) - 2A 3 - ^(Ad^I 1 )] (x - y) 



+ - 



d 4 xd 4 yd 4 9 B(x, 9) [D a D 2 D a B(y, 9)] [A/ 1 ] (x - y) 



9 > j We B(x ,e )BM [-□<*/>) + 2A3 + a- (a^)] <,-„) 

-^ 4 / d 4 xd 4 yd 4 9 B(x,6) [D a D 2 D a B(y, 9)] [A/ 1 ] (ar - 
d 4 xd 4 yd 4 9 B{x, 9) [D a D 2 D a B(y, 9)] [A/ 1 ] (x - y) 
d 4 xd 4 yd 4 9 B{x, 6)B{y, 9) [□(A/ 1 ) - A 3 - ^(Ad^I 1 )] (x - y) 



+ 



+ |- / d 4 xd 4 yd 4 9B(x, 9) [D l3 D 2 D a B(y, 9)} C^H[d p$ , l;l,d Q 



2 



d 4 xd 4 yd 4 9 B(x, 9)B(y, 9)A 3 xy . 



JO) 



Antes de obtener las expresiones renormalizadas, podemos sumar todas las contribuciones 
desnudas que forman la expresion a dos bucles, ya que al emplear RDR todas las estructuras 
se renormalizan siempre con las mismas escalas. En [37], donde se obtuvo la renormal- 
izacion diferencial de SuperQED a dos bucles, este calculo simplificado no se podia realizar, 
ya que habfa que renormalizar cada estructura con su escala correspondiente, para al fi- 
nal relacionarlas mediante las identidades de Ward. La expresion final renormalizada que 
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encontramos es 

( r (2a) , r (2fe) , r (2c) , r (2d 



r 2 , = 2frl 2a) + rf ) + rl 2c) + r (M A 



# 4 y d A xd A yd A e B(x, 9) [D (3 D 2 D a B(y, 6)] C^ & H R [d p$ , 1 ; 1, <9 aci ] 

^ /" j4_ ,4. j4 fl D /.MrnafS2n d/„. ml m ln ( X ~ V?^ 



16(4tt 



2\3 



/ d 4 xd 4 yd 4 9 B(x, 6) \D a D 2 D a B(y, 6)] D ln (^ V) M + 

7 J (x-y) 2 



Ecuacion del grupo de renormalizacion 

Al igual que en QED, evaluando la ecuacion del grupo de renormalizacion a un bucle de los 
campos cuanticos, obtenemos que el valor de la funcion correspondiente a la variacion del 
parametro gauge es de orden g 2 . Esto, unido al hecho de que ni el nivel arbol ni la correction 
a un bucle de la autoenenergia del campo B dependen del parametro gauge, justifica que 
hayamos podido realizar el calculo en el gauge de Feynman. Con la ecuacion del grupo de 
renormalizacion para campos background, obtenemos el valor de la funcion beta a dos bucles 
como 

P(9SQEd) = -^9SQED + 2(4^2)2 5 'SQED + G (9sQEd) , ( 82 ) 

donde empleamos la normalization usual de la constante de acoplamiento, g = V2gsQED- 
Estos valores concuerdan con resultados previos encontrados en la literatura [391 HQl E] ■ 



Ejemplos no abelianos 
Yang-Mills 

El lagrangiano de esta teoria es 

C = ^ v F^ + ^(d,A,) a (d u A v r + (d,mD,v) a , (83) 

donde A a ^ es el campo gauge, 77 y 77 son los fantasmas de Fadeev-Popov, F" = d^A^ — 
d u A^ + gf abc A b Al y f abc las constantes de estructura del algebra de Lie asociada al grupo 
de simetria. Al igual que en los ejemplos previos abelianos, los calculos los realizaremos en 
el gauge de Feynman y con el metodo del campo de background, obteniendo la funcion de 
dos puntos renormalizada del campo background 5" 

Un bucle 

Si estudiamos la autoenergfa del campo background en el gauge de Feynman, obtenemos la 
siguiente expresion desnuda (suma de las contribuciones del bucle de gluones y fantasmas) 

< B;(x)B b u (y) > = g 2 C A 6 ab [4d,d v A 2 -4S^nA 2 + 2d,(Ad v A)-4Ad,d u A] , 

(84) 
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Figure D.4: Diagramas de Yang-Mills a un bucle. 



que se renormaliza como 

<B;(x)B b u (0)> R = g 2 C A 5 a \d,d u -5, u U) 



11 \nx 2 M 2 
-□- 



487r 2 (47r 2 ) x 2 72tt 2 



S(x) 



(85) 



Tambien consideramos la funcion a dos puntos del campo gauge cuantico, ya que sera 
una insertion en uno de los diagramas a dos bucles. El valor total no renormalizado de esta 
funcion es 



< A;(x)Al(y) > 

y aplicando RDR 

< Al(x)A»(0) > R 



g 2 C A 5 ab [d,d u A 2 - 5, v uA 2 + 4d^Ad u A) - 25, u d\Ad x A) 

- AAd^A - VA(DA)] , (86) 



g 2 C A 5 ab 
g 2 C A 5 a 

1447T 2 



\{d,d v - S, u D)A 2 R (x) - ^(d,d„ - S, u n)S(x) 



15 lna; 2 M 2 

47T 2 X 2 



(87) 



Accion efectiva en un gauge generico 



A diferencia que en los ejemplos abelianos, el restringirnos al gauge de Feynman no es inocuo 
en este caso. Para incluir las variaciones del parametro gauge en la ecuacion del grupo de 
renormalizacion, obtendremos la dependencia lineal en £ (el parametro gauge defmido a 
partir del usual como - = 1 + £) de la accion efectiva background expandida a segundo 
orden en los campos B^. Para obtener dicha accion efectiva, consideramos el generador de 
funciones de Green conectadas W 



1 



W = - trln[5 flu n ab -2gf cab B; u + an»V u ) ab } , 



(88) 



con D^ c = d^5 ac + gf abc B b y O ab = (D^D^) ab . A primer orden en £ y segundo orden en B^, 
esto se rescribe como 



W = iC A g 2 tr 



-AB° AR" 

2 [IV 



2ABLAB a uX Ad x d, 



Renormalizando esta expresion, se obtiene facilmente la accion efectiva como 
tC A g 2 



4(4tt 2 ) 



dVy B;(x)B a u (y)(d^ - 5^D)(DA(x - y)) . 



(89) 



(90) 
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Renormalizacion del propagador background a dos bucles 



4 

(a) 





+ 

(b) 






Figure D.5: Diagramas de Yang-Mills a dos bucles (a)-(e). 

Los diagramas (a)-(h) tienen divergencias anidadas, mientras que (i), (j) y (k) correspon- 
den a expresiones con divergencias de solapamiento. En concreto, si evaluamos el diagrama 
(a), tenemos la siguiente expresion desnuda 



< B;(x)B b u (y) > a = ^r c f cd f 9dS P h j d 4 ud 4 v A xy (d; -d*)(dv- ay)A 

x(d v x A uv )A uv (d u x A xu ) , 



yv 



(91) 



que puede ser rescrita en terminos de la integral I 1 previamente definida como 

< B;(x)Bl(y) > a = -g 4 C 2 A 5 ab [Ad^Ad.I 1 ) - d ll d u (AI 1 ) - AAd.dJ 1 } . (92) 

Con los resultados mostrados para I 1 , esto se puede renormalizar de forma inmediata y 
escribir como 

< B a (x)B b u (0) > aR 



9 A C\5 ab 
32(4vr 2 ) 3 



-\\n 2 x 2 M 2 -l\nx 2 M 2 
d„d v D-^ 



\\n 2 x 2 M 2 + %\nx 2 M 2 
+<W □ □ 3 ^ 



r 



+ (termin. locales) 



(93) 
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Figure D.6: Diagramas de Yang-Mills a dos bucles (f)-(k). 

Tomando ahora (i) como un ejemplo de integrales de solapamiento, esta contribution se 
puede escribir en terminos de las integrales H como 



< B*{x)Bl(y) > t = --9 4 C 2 A 5 ab 



dffiHil, d x ; d x , 1] - 2d?H[l, d x ; d x d„, 1] 
-2dy u H[l,d x d,; d x ,l]+4H[l,d,d x ; d u d x ,l]] . 



(94) 



Una vez que tenemos esto, con la lista de expresiones H renormalizadas, llegamos inmedi- 
atamente al siguiente resultado renormalizado 



< b;( x )bI(o) > lR 



g 4 C\5 ab 



32(4vr 



2\3 



-±\n 2 x 2 M 2 -mnx 2 M 2 



±\n 2 x 2 M 2 + fAnx 2 M 2 

+ <W DD - 2^ 

X 1 



+ (termin. locales) . 



(95) 



Procediendo de forma similar en el resto de los diagramas, obtenemos la renormalizacion 
de todas las contribuciones. Sumando todos los resultados, obtenemos el valor renormalizado 
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de la funcion a dos puntos del campo background como 

< Bfa)B b v (0) > R = ~ 2(4 ^2)3 (^A ~ V D ) D x 2 + (termin. locales) 



Ecuacion del grupo de renormalizacion 



(96) 



Evaluando la ecuacion del grupo de renormalizacion a un bucle para la funcion a dos puntos 
del campo cuantico, obtenemos el valor del coeficiente que se ocupa de las variaciones del 
parametro gauge. Dicho valor es 

Entonces, empleando 7^, la accion efectiva background a un bucle en un gauge generico y el 
resultado para la autoenergfa renormalizada del campo B^ a uno y dos bucles, obtenemos 
a partir de la ecuacion del grupo de renormalizacion background el valor de la funcion beta 
como 

P(J9) = i3 ig 3 + P 2 9 5 + 0(g 7 ) 
11CU 



0-. 



48vr 2 



A = "aS- < 98) 



Super Yang-Mills 



Pasamos ahora a estudiar la version supersimetrica del modelo anterior, Super Yang-Mills. 
Al igual que con SuperQED, aplicamos los convenios de [38]. En este caso, la division del 
campo gauge en parte cuantica y parte background es no lineal e 9V <- 3plit ) = e n e gV e n , con 

V el campo cuantico gauge y f2 el prepotencial background. Por lo tanto, las derivadas 
covariantes gauge se escriben en una represent acion quiral cuantica y vectorial background, 
por lo que la accion dividida tiene la forma de 

S = — l -tr ! (e- 9V V a e 9V )V 2 (e- 9V V a e 9V ) , (99) 

2r J 

donde V a es la derivada covariante background. Esto implica que el la parte cuadratica en 

V de la accion con el gauge fijado (de la cual se deriva el propagador cuantico) depende de 
los campos background como 



□ - iW a V a - iW a Va 



donde denotamos el operador cinetico como □ = □ — i\V ai V a — iW ai Va- Otra consecuencia 
del metodo de campo de background en el superespacio es la aparicion de los fantasmas 
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de Nielsen-Kallosh, que corresponden a la normalization de la funcion que se utiliza para 
realizar el promedio sobre los parametros gauge en el procedimiento funcional estandar de 
cuantizacion. Hay que destacar que, como estos fantasmas entran cuadraticamente en la 
action y solo interaccionan con el campo background B (ft = ft = \B), solo contribuyen a 
primer orden en teoria de perturbaciones. 

Debemos notar tambien que para realizar los calculos a dos bucles empleamos supergraficos 
covariantes. Basicamente, realizamos los calculos sin extraer de la derivada covariante la 
conexion espinorial, mediante el algebra de las derivadas covariantes. Por lo tanto, tenemos 
menos graficos y estos son mas convergentes. 

Un bucle 

A la autoenergia background a un bucle solo contribuyen los fantasmas (tanto los de Fadeev- 
Popov como los de Nielsen-Kallosh), siendo la expresion desnuda 



r« = J d 4 xd 4 yd 4 9 B a (x,9)[D 2 D 2 B a (y,9)]A 

M 1 t d 4 xd 4 yd 4 9 B a (x,0)B a (y,e)A xy aA xy 



2 

xy 



+ t3CA j d 4 xd 4 yd 4 9 B a { ^ Q) [ D * D * B a^ ^ A xy d^A xy , (100) 

que se renormaliza de acuerdo con las reglas de RDR como 

= ^^Jd 4 xd 4 yd 4 eB%x,9)[D-D 2 D a B%y,9)]n l ^ X - y)2M2 



(x - yf 



(101) 



En cuanto a la funcion a dos puntos de los campos gauge cuanticos, las contribuciones 
corresponden a los siguientes diagramas: 



WV: AAAv 
(a) 

Figure D.7: Contribuciones a la funcion de dos puntos cuantica a un bucle. 

La contribution final renormalizada en este caso es 
r (i) 




*' = -||/ Aiy,1, ' (l ' ,|DV * W)D w ! ' (102) 
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Accion en un gauge generico 

Al igual que en caso de Yang-Mills, tenemos que evaluar la accion efectiva background 
en un gauge generico para tener en cuenta el termino de variation del parametro gauge £ 
(redefinido a partir del usual como £ + 1 = -) en la ecuacion del grupo de renormalizacion. 
Obtendremos el termino lineal en £ correspondiente a la contribucion a segundo orden en 
campos background de la expansion de la accion efectiva. Para ello, a partir de un calculo 
funcional, escribimos dicha accion como 



eff 



--tr In 
2 



□ + £(v 2 v 2 + V 2 V 2 ) 



— tr\nn + tr\nD_ + T € , 
2 



+ tr In 



□_ + £V 2 V 2 



(103) 



donde 



□_ 



□ -iW a V n - -(V a W r 



(104) 



Considerando la expresiones inversas de los operadores, y quedandonos a segundo orden 
en campos background, obtenemos el valor renormalizado como 



J^y tr J d^yd 20 W a (x,6)W a (y,6)n 



lm> - yfM? R 
(x - y) 2 



(105) 



Dos bucles 

Para realizar la renormalizacion a dos bucles con supergraficos covariantes, simplemente 
tenemos que considerar el siguiente diagrama de vacio: 



Figure D.8: Contribucion a dos bucles a la accion efectiva background. 

En dicho diagrama, los propagadores gauge □ 1 dependen de los campos background, 
por lo que hemos de expandirlos y quedarnos con el segundo orden en B. De esta expansion 
tenemos dos tipos de contribuciones. Unas tienen campos W a expHcitos, mientras que otras 
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se construyen en base a las conexiones espacio-temporales background T a a. Asi, tenemos los 
siguients resultados renormalizados para las contribuciones con campos W a 



(2) | 



\R 



2>%g 2 C 2 



i=i 



/ 



tr I d 4 xd 4 yd 4 6 



W a (x, 9)81^ + W\x, 6)d y a& W a (y, 9) 



x 



R 
2ni 



3ig 2 C 



16(4tt 2 ) j 



x \n(x-y) 2 M 2 



tr / d 4 xd 4 yd 2 6W a (x,6)W a (y,6)d x aA 



(x - y)< 



+ 0(B 3 ) , 
(106) 



la cual, mediante una identidad de Bianchi se puede escribir en una forma invariante gauge 
como 

2 

E^l* = tfCltrj d 4 xd 4 yd 2 9W a (x,9)W a (y : 9)n[AI°} R (x-y) 

ln(x - y) 2 M 2 



i=i 



3 ^ r tr / d 4 xd 4 yd 2 Q W a {x,9)W a {y,9)W- 



16(4vr 



2\3 



(x - yf 



+ 0(B 3 ) . 



(107) 



Por otro lado, la suma de los diagramas con conexiones espacio-temporales es 

£lf>| fl = -Zg 2 C\tr I ^^V^r-(x^)r^(y^)(^^-(2C a/3 C' d/ j)n) 



-[A/%(x-y)- 



1 ln(x-y) 2 M 2 



+ 0(B 3 ) , 



(108) 



4 L— j«v- 32 ( 47r 2)3 ( X - 2/ )2 

que, debido a ser una expresion transversa se puede escribir en terminos de W a como 
tr J d 4 xd 4 yd 4 9 T™(x, 9)T^(y, 6) (d^ - 2C Q/3 C^n) /(x - y) 
= -3tr J d 4 xd 4 yd 4 9 [D a B(x,9)] [D 2 D a B(y, 9)] af(x-y) + 0(B 3 ) 

^rjWeW^W^U^-y)^). (109) 

Por lo tanto, la expresion total renormalizada a dos bucles y segundo orden en los campos 
background es 

I^rf } = tr J d 4 xd 4 yd 2 6W a {x,6)W a {y,6)T i2 \x-y) , (110) 



con 



V m {x) = 3fl 2 Cj D |ln 2 x 2 M| fi + |lnx 2 M f yi -lnx 2 M 2 ) + lnx 2 M 2 



64(4tt 2 ) 3 
+ (termin. locales) 



(111) 
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Ecuacion del grupo de renormalizacion 

A la hora de evaluar la ecuacion del grupo de renormalizacion, los pasos que hemos de dar 
son identicos al caso de Yang-Mills. En primer lugar, con la ecuacion correspondiente a la 
funcion a dos puntos de los campos cuanticos a un bucle, obtenemos el valor del coeficiente 
de variation del parametro gauge, •y^d/dC,. En concreto, tenemos el siguiente resultado 

-» = ~wk/ +0 ^- (112) 

Entonces, con 7^, la action efectiva en una gauge generico y la contribution a uno y dos 
bucles a la autoenergfa background, de la ecuacion del grupo de renormalizacion satisfecha 
por los campos B 



M m + ^ 9) h +ld9) k 



T(x) 



= (113) 

C=o 



obtenemos la expansion de la funcion beta como 

(3(9sym) = -(3/2)[C A /(8n 2 )]gl YM - (3/2)[C A /(8n 2 )] 2 gl YM + 0(g 7 SYM ) , (114) 

donde hemos utilizado, igual que en el caso de SuperQED, la constante de acoplamiento 
usual, que difiere en un factor a/2 de la que se emplea en [38]. Notar que, al igual que 
en el resto de las teorfas consideradas, no tenemos dimension anomala (7b), ya que con la 
normalization que tenemos del campo background y la relation entre la renormalizacion de 
dicho campo y la constate de acoplamiento, jb se anula. 

Este calculo nos permite dar una nueva vision sobre un punto controvertido: el origen de 
las correcciones mas alia de un bucle a la funcion beta de Super Yang-Mills. En un princi- 
pio, Novikov, Shifman, Vainshtein and Zakharov (NSVZ) en [8] obtuvieron la "funcion beta 
exacta" (Pnsvz) empleando un calculo de instantones, siendo posteriormente reobtenido 
el coeficiente a dos bucles de esta funcion mediante calculos perturbativos (empleando re- 
duction dimensional) [551 EE E3] . Aunque algunos calculos paretian indicar un origen IR a 
las correcciones de orden superior de Pnsvz [HI HH E3j, esto fue cuestionado en [9], [10], donde 
empleando un formalismo wilsoniano (por lo tanto, en principio solo dependiente del compor- 
tamiento UV de la teorfa) y diferenciando entre una constante de acoplamiento holomorfica 
y una canonica, se obtuvo un ffujo NSVZ para la ultima. 

Nuestro calculo tiene la virtud de evitar uno de los puntos conflictivos de la aplicacion 
de metodos dimensionales a Super Yang-Mills: la regularizacion de tanto las divergencias 
UV como IR con el mismo parametro infinitesimal, que implica que se mezclen ambas con- 
tribuciones en los resultados renormalizados. Nosotros en cambio tenemos las divergencias 
claramente diferenciadas al tener asociadas dos escalas independientes. Lo que hemos en- 
contrado es que la escala correspondiente a la renormalizacion a un bucle es la que genera la 
contribution a dos bucles de Pnsvz- No existe una escala a dos bucles UV (lo cual coincide 
con la conclusion obtenida en [56J, conforme a la cual, en un esquema de regularizacion en 
cuatro dimensiones no hay divergencias superficiales), aunque esto no implica que el coefi- 
ciente a dos bucles de la funcion beta sea nulo. Como se ve en este caso, la escala UV a un 
bucle sobrevive a dos bucles al tener en cuenta los efectos IR. 
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